Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 3,178,277 G→T E27* (GAA→TAA)  yqiC → BMFP family putative fusogenic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,178,2770GT92.8% 263.8 / NA 99E27* (GAA→TAA) yqiCBMFP family putative fusogenic protein
Reads supporting (aligned to +/- strand):  ref base G (3/4);  new base T (36/54);  total (40/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGC  >  CP009273/3178131‑3178420
                                                                                                                                                  |                                                                                                                                               
tCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTaa                                                                                                                                               >  1:434514/1‑149 (MQ=255)
tCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTaa                                                                                                                                               <  2:1779678/149‑1 (MQ=255)
tCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTaa                                                                                                                                               >  1:1571706/1‑149 (MQ=255)
  tGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAAGGCCTAAAGGAATCAGGGAGTTCGGGTAAGa                                                                                                                                             >  2:1827581/1‑149 (MQ=255)
    gTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGAtg                                                                                                                                           <  2:1481314/149‑1 (MQ=255)
     tttAAGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGAtgt                                                                                                                                          <  2:507971/149‑1 (MQ=255)
      ttATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGGTTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGCTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGAtgtg                                                                                                                                         <  1:1337717/149‑1 (MQ=255)
       tATGACTTACCCCTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTgg                                                                                                                                        <  1:1005425/149‑1 (MQ=255)
       tATGACTACCCCTTATCGCACTTCAATGGCAATCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGAAAATCGCTCGCCAGGTTCACGAATCAATGCCAAAAGGAATCAGGGAGTTAGGGTAAGATGTgg                                                                                                                                        >  1:355290/1‑149 (MQ=255)
        aTGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCGAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGATCGGAGAAGATGTggg                                                                                                                                       >  2:912719/1‑148 (MQ=255)
        aTGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTCCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGga                                                                                                                                       <  2:416852/149‑1 (MQ=255)
        aTGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTAAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAATGGAATCAGGGAGTTCGGGTAAGATGTGga                                                                                                                                       <  2:1189132/149‑1 (MQ=255)
          gACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCCCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGgaga                                                                                                                                     >  2:1923540/1‑149 (MQ=255)
             ttATCCTTATCGCACACCAATGGCACTCAACACCAATCAGGACAGGGAAACTACAAGGAAAGAACCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGaaaa                                                                                                                                  <  1:166955/149‑1 (MQ=255)
             ttACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGGCCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGGATCAGGGAGTTCGGGTAAGATGTGGAGaaaa                                                                                                                                  >  2:1302430/1‑149 (MQ=255)
             ttACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGTCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGaaaa                                                                                                                                  >  1:904703/1‑149 (MQ=255)
             ttACCCTTATCGCACTACAATGGCACTCAACACATATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCAGGTAAGATGTGGAGaaaa                                                                                                                                  >  2:1242823/1‑149 (MQ=255)
                  cTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCATAAAGGAATCAGGGAGTTCGGAGAAGATGTGGAGAAAAAAAta                                                                                                                             >  1:1157204/1‑148 (MQ=255)
                   ttAGCGCACGACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCCAATCGCTCGCCAGGTTCACGAAAAAATGCCTAAAAGAATCAGGGCGTTCGGGAAAGGTGTGGAGAAAAAAATcc                                                                                                                            >  1:257117/1‑149 (MQ=255)
                   ttACCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATcc                                                                                                                            <  1:1548835/149‑1 (MQ=255)
                           cTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGAAGTTCGGGTAAGATGTGGAGAACAGAATCCGCCAAAcc                                                                                                                    >  2:822093/1‑149 (MQ=255)
                              cAACGGCACTCAACATCTATCAGCACATGGAAACTACAATGATTGGCCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTa                                                                                                                 <  1:5701/149‑1 (MQ=255)
                                 tGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACaa                                                                                                              <  2:368248/149‑1 (MQ=255)
                                  ggCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAg                                                                                                             <  2:1859812/149‑1 (MQ=255)
                                    cACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGAGGTGGAGAAAAAAATCCGCCAAACCCTACAAgcg                                                                                                           >  2:890060/1‑149 (MQ=255)
                                     aCTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAgcgc                                                                                                          >  2:540860/1‑149 (MQ=255)
                                       tCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAg                                                                                                        <  1:147979/149‑1 (MQ=255)
                                       tCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGGGGAGAAAAAAATCCGCCAACCCCTACAAGCGCAg                                                                                                        >  2:1715374/1‑149 (MQ=255)
                                       tCAACACCTATCAGAACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAACCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAg                                                                                                        <  2:54157/149‑1 (MQ=255)
                                       tCAACACCTACCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAg                                                                                                        <  2:412598/149‑1 (MQ=255)
                                          acacCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAACAAAATCCGCCAAACCCTACAAGCGCAGCTg                                                                                                     >  2:321765/1‑149 (MQ=255)
                                          acacCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTg                                                                                                     >  2:1275750/1‑149 (MQ=255)
                                              cTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGAcgc                                                                                                 >  2:1741627/1‑149 (MQ=255)
                                              cTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGAcgc                                                                                                 >  2:1817756/1‑149 (MQ=255)
                                                  cAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCt                                                                                             <  2:1910869/149‑1 (MQ=255)
                                                     tACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGa                                                                                          <  2:1043377/149‑1 (MQ=255)
                                                     caCAGGGAAACTACAATGCTGGACCCGAAAAAAATTGGGCAATGCGCCCGCCAGGTTCACGAATCAATGCCTAACGGAAGCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAACCCCTGCAAGCGCAGCTGACGCGCCTCGa                                                                                          <  2:1694901/148‑1 (MQ=255)
                                                        aGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTAACGAATAAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATACGCAAAACCCTACAAGCGAAGCTGAAGCGCTTCGACCt                                                                                       >  1:168738/1‑149 (MQ=255)
                                                            aaaCTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTa                                                                                   <  1:116753/149‑1 (MQ=255)
                                                               cTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGc                                                                                <  2:1158884/149‑1 (MQ=255)
                                                                tACAATGATTGACCCGAAAAACATTGAACAAATCGCTCGCCAGGTACCCGAATCAATGCCAAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCGTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGcc                                                                               >  1:677124/1‑149 (MQ=255)
                                                                 aCAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTCAGCCg                                                                              >  1:367212/1‑149 (MQ=255)
                                                                  ctaTGAGTGACCCGAAAAAAATTGAGCCAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGt                                                                             <  1:1146319/147‑1 (MQ=255)
                                                                  cAATGATTGACCCGCAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGGTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGACTCGGCCTTGTAAGCCgg                                                                             >  1:524443/1‑148 (MQ=255)
                                                                   aagggttgcgccGAAAAAAAGTGAGCAGATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTg                                                                            <  1:395350/139‑1 (MQ=255)
                                                                   aaTGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTg                                                                            <  2:1925994/149‑1 (MQ=255)
                                                                   aaTGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTg                                                                            >  2:1490874/1‑149 (MQ=255)
                                                                        tttaCCCGAAAAAAAGTGAGCAAATCGCTCGCGAGGTTCACGAATCAATGCCGAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAagag                                                                       <  1:1399130/146‑1 (MQ=255)
                                                                           aCCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTc                                                                    >  2:705062/1‑149 (MQ=255)
                                                                           aCCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTc                                                                    <  1:795943/149‑1 (MQ=255)
                                                                            cccGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCg                                                                   >  2:819010/1‑149 (MQ=255)
                                                                              cGAAAAAAATTGAGCAAGTCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGTAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAACCCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTCAGCCGTGAAGAGTTCGAc                                                                 >  1:1323222/1‑149 (MQ=255)
                                                                              cGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTTGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGCGTTCAAc                                                                 >  1:272608/1‑149 (MQ=255)
                                                                               gAAATAAATTGATCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTGAGCCGGGAAGAGTTCGACg                                                                >  1:1359389/1‑149 (MQ=255)
                                                                                aaaaaaaTTGCGCAAAACGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGt                                                               <  1:1623324/149‑1 (MQ=255)
                                                                                aaaaaaaTTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGAAGTTCGGGTACGACGTGGAGAAAAAAATCCACCAAACCCTACAAGCGCACCTGACGCGCCTCGCCCTTGTAAACCGCGAagctaccgccga                                                               >  1:1789872/1‑139 (MQ=255)
                                                                                 aaaaaaTTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTc                                                              <  2:776895/149‑1 (MQ=255)
                                                                                    aaaTTGAGCAAATCGCCCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCaa                                                           <  2:1284658/149‑1 (MQ=255)
                                                                                     aaTGGAGCGAATCGCTCGCCAGGTGCACGAATCAATGGCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCaaa                                                          <  2:1596152/149‑1 (MQ=255)
                                                                                      aTTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGGAAGCCGTGAAGAGTTCGACGTCCAAAc                                                         <  2:389047/149‑1 (MQ=255)
                                                                                       ttGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACg                                                        >  2:965118/1‑149 (MQ=255)
                                                                                           gCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATGAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAg                                                    <  1:855580/149‑1 (MQ=255)
                                                                                               aTCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACACGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCt                                                <  2:492708/149‑1 (MQ=255)
                                                                                               aTCGCTCGCCAGGTTCACGAAGCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCt                                                <  1:1625090/149‑1 (MQ=255)
                                                                                                  gcccgcCACGTTCCCGAATCAATGCCTAAAGGAATCAGGAACTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGtt                                             <  2:1617180/146‑1 (MQ=255)
                                                                                                   cacgcCAGGTCCACGAATCAATGCCTAAAGTAATCAGGGCGTTGGGGGAAGATGGGGAGAAAAAAATCCGCCAAACCCTATAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTa                                            <  1:145496/147‑1 (MQ=255)
                                                                                                    tcgcCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAAGGCAAGTCCTGTTAc                                           >  2:820398/1‑149 (MQ=255)
                                                                                                       ccAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGTCGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGGccgtaacagga                                        >  1:599530/1‑140 (MQ=255)
                                                                                                            ttCACGAATCAATGCCCAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGt                                   <  2:830039/149‑1 (MQ=255)
                                                                                                               gcGACTCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAAGCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGaa                                <  1:1656603/148‑1 (MQ=255)
                                                                                                               aCGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGaa                                >  2:872837/1‑149 (MQ=255)
                                                                                                                  aaTCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGaaaaa                             >  1:270002/1‑149 (MQ=255)
                                                                                                                  aaTCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGaaaaa                             <  2:368832/149‑1 (MQ=255)
                                                                                                                  aaTCAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGaaaaa                             >  1:828536/1‑149 (MQ=255)
                                                                                                                     cAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTg                          <  2:256344/149‑1 (MQ=255)
                                                                                                                     cAATGCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTg                          <  1:1791851/149‑1 (MQ=255)
                                                                                                                      aaggCCTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGGGAAAAAAATCCGCCAAACCCTACAGGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTgg                         <  1:874583/146‑1 (MQ=255)
                                                                                                                      aaTGCCTAAAGGAATCAGGGGGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTgg                         <  2:1689599/149‑1 (MQ=255)
                                                                                                                       aTGCCTAAAGGAATCAGGGAGTTTGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGgc                        <  2:295857/149‑1 (MQ=255)
                                                                                                                           ccaaaGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGTAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCgctg                    <  1:1138805/147‑1 (MQ=255)
                                                                                                                           cTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCgctg                    <  2:1385444/149‑1 (MQ=255)
                                                                                                                           cTAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCgctg                    <  1:1195209/149‑1 (MQ=255)
                                                                                                                            tAAAGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCgctgc                   >  2:708332/1‑149 (MQ=255)
                                                                                                                            tAAAGGAAGCAGGGGGTTCGGGTAAGAGGTGGAGAAAAAAACCCGCCCAACCCTACCAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCgctgc                   <  1:406989/149‑1 (MQ=255)
                                                                                                                             aaaGGAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCgctgct                  >  2:442020/1‑149 (MQ=255)
                                                                                                                              gagcaaTGAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCtt                 <  1:837877/145‑1 (MQ=255)
                                                                                                                              aaggacaccggcaGTTCGGGCAAGACGCGTAGAAAAAAATCCGCCAAACCCGACAAGCGAAGCTGACGCGCCCCGACCTGGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCtt                 <  2:255365/137‑1 (MQ=255)
                                                                                                                                 gAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGGCTCGTGAAAAACTGGCGCTGCTGGag              >  1:275284/1‑148 (MQ=255)
                                                                                                                                 gAATCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGaa              <  1:1690042/149‑1 (MQ=255)
                                                                                                                                 gAATAAGGGAGTTCGGGTAAGAAGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGTTGACGCGCCTCGACCTTGTAAGACGTCAAGATTTCGACGTCCAAACGCAAGTCCTGTTACGAACTCGTGAAAAACTGGCGCTGCTTGaa              <  1:202100/149‑1 (MQ=255)
                                                                                                                                  aaTCAGGGAGTTCGGTGAAGAGGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGCCCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCCTGacc             <  1:747401/149‑3 (MQ=255)
                                                                                                                                    tCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACaa           <  1:807905/149‑1 (MQ=255)
                                                                                                                                    tCAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGAGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACaa           <  2:870749/149‑1 (MQ=255)
                                                                                                                                     cAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCGAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACaaa          <  2:1575652/149‑1 (MQ=255)
                                                                                                                                     cAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACaaa          <  2:619454/149‑1 (MQ=255)
                                                                                                                                     cAGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCGACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACaaa          <  1:1711806/149‑1 (MQ=255)
                                                                                                                                      aGGGAGTTCGGGTAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGGACGAAg         >  2:1639775/1‑149 (MQ=255)
                                                                                                                                       gggAGGTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCGGACGCGCCTCGACCTTGTACGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGa        <  1:1186886/149‑1 (MQ=255)
                                                                                                                                             tCGGGTAAGATGTGGAGAAAAAAATCCGCCCAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGc  <  2:91746/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TCTGGTTTATGACTTACCCTTATCGCACTACAATGGCACTCAACACCTATCAGTACAGGGAAACTACAATGATTGACCCGAAAAAAATTGAGCAAATCGCTCGCCAGGTTCACGAATCAATGCCTAAAGGAATCAGGGAGTTCGGGGAAGATGTGGAGAAAAAAATCCGCCAAACCCTACAAGCGCAGCTGACGCGCCTCGACCTTGTAAGCCGTGAAGAGTTCGACGTCCAAACGCAAGTCCTGTTACGGACTCGTGAAAAACTGGCGCTGCTTGAACAAAGAATCAGC  >  CP009273/3178131‑3178420

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: