Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,730,144 A→T 8.8% A170A (GCA→GCT sodB → superoxide dismutase, Fe

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,730,1440AT 8.8% 733.4 / 47.9 239A170A (GCA→GCTsodBsuperoxide dismutase, Fe
Reads supporting (aligned to +/- strand):  ref base A (113/105);  new base T (14/7);  total (127/112)
Fisher's exact test for biased strand distribution p-value = 2.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 246 total aligned reads displayed.

CTGGCTGGACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGAT  >  CP009273/1729996‑1730292
                                                                                                                                                    |                                                                                                                                                    
ctggctggACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCa                                                                                                                                                      <  2:1606962/149‑1 (MQ=255)
 tggctggACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCAc                                                                                                                                                     <  2:1917253/149‑1 (MQ=255)
  ggctggACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTCCTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACAACCACTATCGCAATGCaag                                                                                                                                                    >  2:750752/1‑147 (MQ=255)
  ggctggACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCTCg                                                                                                                                                    >  2:3818246/1‑149 (MQ=255)
     tggACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTcc                                                                                                                                                 >  1:598698/1‑149 (MQ=255)
     tggACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTcc                                                                                                                                                 <  1:322371/149‑1 (MQ=255)
     tggACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTcc                                                                                                                                                 <  1:226471/149‑1 (MQ=255)
      ggACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCt                                                                                                                                                >  1:3458648/1‑149 (MQ=255)
        aCCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTgg                                                                                                                                              >  2:1156676/1‑149 (MQ=255)
          ctggctggTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCt                                                                                                                                            <  2:844747/149‑1 (MQ=255)
              ctggTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCt                                                                                                                                        <  2:1320762/149‑1 (MQ=255)
               tggTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTg                                                                                                                                       <  2:176360/149‑1 (MQ=255)
                ggTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTgg                                                                                                                                      >  1:2779064/1‑149 (MQ=255)
                 gTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCTCGTCCTGGCTATCTGGa                                                                                                                                     >  1:963967/1‑149 (MQ=255)
                  tGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAg                                                                                                                                    >  1:1194217/1‑149 (MQ=255)
                     aaaaCAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCAc                                                                                                                                 >  1:1397/1‑149 (MQ=255)
                     aaaaCAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCAc                                                                                                                                 <  1:1911641/149‑1 (MQ=255)
                      aaaCAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCTCGTCCTGGCTATCTGGAGCACt                                                                                                                                >  2:2966897/1‑149 (MQ=255)
                      aaaCAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACt                                                                                                                                >  1:1704761/1‑149 (MQ=255)
                          aGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCTCAATGCTCGTCCTGGCTATCTGGAGCACTTCTg                                                                                                                            >  1:346748/1‑149 (MQ=255)
                          aGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCTCGTCCTGGCTATCTGGAGCACTTCTg                                                                                                                            >  1:2693884/1‑149 (MQ=255)
                           gCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTgg                                                                                                                           >  2:259752/1‑149 (MQ=255)
                           gCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTgg                                                                                                                           >  2:2001389/1‑149 (MQ=255)
                               tGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGgcgc                                                                                                                       <  2:911978/149‑1 (MQ=255)
                                ggCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCt                                                                                                                      <  1:124356/149‑1 (MQ=255)
                                   aaaCTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGt                                                                                                                   >  2:3573654/1‑149 (MQ=255)
                                   aaaCTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGt                                                                                                                   >  1:3798041/1‑149 (MQ=255)
                                    aaCTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTg                                                                                                                  <  1:1300936/149‑1 (MQ=255)
                                     aCTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGa                                                                                                                 >  2:918086/1‑149 (MQ=255)
                                     aCTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGa                                                                                                                 <  2:576518/149‑1 (MQ=255)
                                       tGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAAc                                                                                                               <  2:1276607/149‑1 (MQ=255)
                                       tGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAAc                                                                                                               >  1:3722103/1‑149 (MQ=255)
                                        ggCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACt                                                                                                              <  1:1008904/149‑1 (MQ=255)
                                            aTCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGa                                                                                                          >  1:4085927/1‑149 (MQ=255)
                                             tCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGaa                                                                                                         <  2:3054118/149‑1 (MQ=255)
                                               ntttCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAAtt                                                                                                       >  1:3119716/2‑149 (MQ=255)
                                               gTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAAtt                                                                                                       <  1:1656127/149‑1 (MQ=255)
                                               gTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAAtt                                                                                                       <  1:374584/149‑1 (MQ=255)
                                                 ttCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCg                                                                                                     <  2:3336086/149‑1 (MQ=255)
                                                  tCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGt                                                                                                    <  2:4008243/149‑1 (MQ=255)
                                                   cAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTa                                                                                                   <  1:2840971/149‑1 (MQ=255)
                                                    aaCCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAg                                                                                                  >  2:2312303/1‑149 (MQ=255)
                                                     aCCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGc                                                                                                 >  1:3449129/1‑149 (MQ=255)
                                                           aaCGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGaaaaa                                                                                           >  2:3442562/1‑149 (MQ=255)
                                                            aCGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAAt                                                                                          >  2:2241258/1‑149 (MQ=255)
                                                                 ggTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGACTTCGTAGCGAAAAATCTCGc                                                                                     >  1:1934774/1‑149 (MQ=255)
                                                                 ggTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGc                                                                                     >  1:2841559/1‑149 (MQ=255)
                                                                 ggTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGc                                                                                     <  2:3500404/149‑1 (MQ=255)
                                                                   tACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTg                                                                                   >  2:1701654/1‑149 (MQ=255)
                                                                         gCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCataat                                                                             >  2:3476284/1‑149 (MQ=255)
                                                                           tGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCTCGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCAtaataa                                                                           >  1:771308/1‑149 (MQ=255)
                                                                            gACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAAc                                                                          <  1:814422/149‑1 (MQ=255)
                                                                              ccaccGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTg                                                                        <  2:869186/149‑1 (MQ=255)
                                                                                  cGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATgg                                                                    <  1:3011015/149‑1 (MQ=255)
                                                                                   gATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGc                                                                   >  2:3206659/1‑149 (MQ=255)
                                                                                   gATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGc                                                                   <  2:3634036/149‑1 (MQ=255)
                                                                                       cGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAAt                                                               <  1:3204268/149‑1 (MQ=255)
                                                                                         aCTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATgc                                                             <  2:672583/149‑1 (MQ=255)
                                                                                            ccGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATgcagc                                                          <  1:2504610/149‑1 (MQ=255)
                                                                                             cGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATgcagca                                                         <  1:3404321/149‑1 (MQ=255)
                                                                                              gctgctGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCAt                                                        <  1:3587982/149‑1 (MQ=255)
                                                                                              gctgctGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCAt                                                        >  2:2044476/1‑149 (MQ=255)
                                                                                               ctgctgACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCAtt                                                       <  1:686751/149‑1 (MQ=255)
                                                                                                 gctgACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGc                                                     <  1:1513012/149‑1 (MQ=255)
                                                                                                 gctcACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGACCTGGCTATCTGGAGCACTTCTGGGCGCTGGGGAACTGGCAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATCCAGCATTGc                                                     >  2:3955500/1‑149 (MQ=255)
                                                                                                  ctgACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGcc                                                    <  2:1808355/149‑1 (MQ=255)
                                                                                                  ctgACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGcc                                                    <  2:3194515/149‑1 (MQ=255)
                                                                                                  ctgACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGcc                                                    <  2:3850244/149‑1 (MQ=255)
                                                                                                   tgACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCt                                                   <  1:470864/149‑1 (MQ=255)
                                                                                                     aCCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGa                                                 <  1:2756247/149‑1 (MQ=255)
                                                                                                      ccGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGaa                                                >  2:1200771/1‑149 (MQ=255)
                                                                                                      ccGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGaa                                                >  1:361321/1‑149 (MQ=255)
                                                                                                       cGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAg                                               <  1:2250345/149‑1 (MQ=255)
                                                                                                         ttGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCTCGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAgcg                                             >  2:1839377/1‑149 (MQ=255)
                                                                                                             tGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTAc                                         >  2:496555/1‑149 (MQ=255)
                                                                                                                  gggAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCtaa                                    >  2:3755379/1‑147 (MQ=255)
                                                                                                                   ggAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTAt                                   <  1:2588129/149‑1 (MQ=255)
                                                                                                                   ggAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTAt                                   >  1:2308149/1‑149 (MQ=255)
                                                                                                                      acacGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAg                                >  1:1343615/1‑149 (MQ=255)
                                                                                                                      acacGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGAGCTACGCTTATCAg                                >  1:560406/1‑149 (MQ=255)
                                                                                                                          gCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCt                            <  1:2975693/149‑1 (MQ=255)
                                                                                                                            ttattaCATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTAc                          >  2:2417124/1‑149 (MQ=255)
                                                                                                                            ttattaCATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTAc                          >  1:229843/1‑149 (MQ=255)
                                                                                                                               ttaCATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTAcgcg                       >  2:3895093/1‑149 (MQ=255)
                                                                                                                                taCATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCgg                      <  2:352931/149‑1 (MQ=255)
                                                                                                                                  cATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGat                    >  1:2449818/1‑149 (MQ=255)
                                                                                                                                   aTCGACTATCGCAATGCTCGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGatc                   >  2:3318967/1‑149 (MQ=255)
                                                                                                                                    tCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGatca                  <  1:3886474/149‑1 (MQ=255)
                                                                                                                                     cGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGatcat                 <  1:3709480/149‑1 (MQ=255)
                                                                                                                                       aCTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCg               >  2:1076235/1‑149 (MQ=255)
                                                                                                                                       aCTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGAGCATCg               >  2:1244680/1‑149 (MQ=255)
                                                                                                                                          aTCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATg            <  2:2841307/149‑1 (MQ=255)
                                                                                                                                           tCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGt           >  2:2379116/1‑149 (MQ=255)
                                                                                                                                           tCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGt           >  2:932392/1‑149 (MQ=255)
                                                                                                                                             gCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAg         <  2:1582845/149‑1 (MQ=255)
                                                                                                                                               aaTGCTCGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGt       >  1:2078427/1‑149 (MQ=255)
                                                                                                                                                aTGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTc      <  2:4074171/149‑1 (MQ=255)
                                                                                                                                                aTGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTc      <  1:839564/149‑1 (MQ=255)
                                                                                                                                                 tGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCg     <  2:2779056/149‑1 (MQ=255)
                                                                                                                                                 tGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCg     >  2:266109/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGGCTGGACCTGGCTGGTGAAAAACAGCGATGGCAAACTGGCTATCGTTTCAACCTCTAACGCGGGTACTCCGCTGACCACCGATGCGACTCCGCTGCTGACCGTTGATGTCTGGGAACACGCTTATTACATCGACTATCGCAATGCACGTCCTGGCTATCTGGAGCACTTCTGGGCGCTGGTGAACTGGGAATTCGTAGCGAAAAATCTCGCTGCATAATAACTGATGGCAAATGCAGCATTGCCTGAAGCGCTACGCTTATCAGGCCTACGCGGATCATCGATGTAGGTCGGAT  >  CP009273/1729996‑1730292

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