Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,959,915 C→A 89.9% P178P (CCG→CCT)  flhB ← flagellin export apparatus, substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,959,9150CA89.9% 374.0 / 39.0 149P178P (CCG→CCT) flhBflagellin export apparatus, substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base C (8/7);  new base A (81/53);  total (89/60)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

Only 100 of 149 total aligned reads displayed.

AATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGG  >  CP009273/1959771‑1960058
                                                                                                                                                |                                                                                                                                               
aaTATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGga                                                                                                                                              <  1:2369916/148‑1 (MQ=255)
  tatCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGgaga                                                                                                                                            >  1:148679/1‑148 (MQ=255)
  tatCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGgaga                                                                                                                                            >  1:2059638/1‑148 (MQ=255)
     ccTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTc                                                                                                                                         >  1:3123285/1‑148 (MQ=255)
     ccTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTc                                                                                                                                         >  2:178390/1‑148 (MQ=255)
       tGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCgg                                                                                                                                       >  2:2440435/1‑148 (MQ=255)
            tGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATc                                                                                                                                  >  2:2207163/1‑148 (MQ=255)
             gACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCa                                                                                                                                 <  2:2333186/148‑1 (MQ=255)
              aCATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCaa                                                                                                                                <  2:1501624/148‑1 (MQ=255)
                 tACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAgcg                                                                                                                             >  1:2303989/1‑148 (MQ=255)
                  aCGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAgcgc                                                                                                                            <  1:1598832/148‑1 (MQ=255)
                  aCGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAgcgc                                                                                                                            <  1:2220464/148‑1 (MQ=255)
                        cTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCatcatc                                                                                                                      >  2:1342335/1‑148 (MQ=255)
                          ttttCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTg                                                                                                                    <  1:1492158/148‑1 (MQ=255)
                          ttttCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTg                                                                                                                    <  1:2310061/148‑1 (MQ=255)
                                   ggCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTg                                                                                                           <  2:1689418/148‑1 (MQ=255)
                                     cTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGAt                                                                                                         >  1:1012442/1‑148 (MQ=255)
                                     cTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGAt                                                                                                         >  2:1062127/1‑148 (MQ=255)
                                       gAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAACTCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGc                                                                                                       <  2:2740794/148‑1 (MQ=255)
                                       gAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGc                                                                                                       >  2:2657198/1‑148 (MQ=255)
                                       gAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGc                                                                                                       >  2:478499/1‑148 (MQ=255)
                                           atTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCaga                                                                                                   >  1:405505/1‑148 (MQ=255)
                                           atTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCaga                                                                                                   >  1:1153205/1‑148 (MQ=255)
                                           atTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCaga                                                                                                   >  2:3047349/1‑148 (MQ=255)
                                                gAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGaaaa                                                                                              >  1:1208751/1‑148 (MQ=255)
                                                gAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGaaaa                                                                                              >  2:2337538/1‑148 (MQ=255)
                                                gAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGaaaa                                                                                              <  2:1747808/148‑1 (MQ=255)
                                                  aaaaGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGaaaaaa                                                                                            >  2:1794827/1‑148 (MQ=255)
                                                   aaaGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAAc                                                                                           >  1:2231387/1‑148 (MQ=255)
                                                    aaGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCCTGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAAcc                                                                                          <  1:1287616/148‑1 (MQ=255)
                                                       aCGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCg                                                                                       >  1:1707061/1‑148 (MQ=255)
                                                          tCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTc                                                                                     <  1:3310952/147‑1 (MQ=255)
                                                               tCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTg                                                                               >  2:1794245/1‑148 (MQ=255)
                                                               tCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTg                                                                               <  2:1883532/148‑1 (MQ=255)
                                                                cccACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGc                                                                              >  2:325522/1‑148 (MQ=255)
                                                                 ccaccaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGcc                                                                             <  1:856270/148‑1 (MQ=255)
                                                                 ccaccaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGcc                                                                             >  2:2639035/1‑148 (MQ=255)
                                                                 ccaccaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGcc                                                                             <  1:856198/148‑1 (MQ=255)
                                                                  caccaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCa                                                                            >  2:298600/1‑148 (MQ=255)
                                                                   accaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCaa                                                                           >  1:3181892/1‑148 (MQ=255)
                                                                    ccaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAAc                                                                          <  1:437055/148‑1 (MQ=255)
                                                                    ccaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAAc                                                                          <  2:1667680/148‑1 (MQ=255)
                                                                    ccaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAAc                                                                          >  2:3252462/1‑148 (MQ=255)
                                                                    ccaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCCCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAAc                                                                          <  1:1695813/148‑1 (MQ=255)
                                                                     caTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAAcc                                                                         <  1:201627/148‑1 (MQ=255)
                                                                     caTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAAcc                                                                         <  1:969791/148‑1 (MQ=255)
                                                                        tGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGACAAAACCCCGTCACGCTGCCAACCAgg                                                                      >  2:1432333/1‑148 (MQ=255)
                                                                        tGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAgg                                                                      >  2:768569/1‑148 (MQ=255)
                                                                        tGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAgg                                                                      >  2:1386911/1‑148 (MQ=255)
                                                                        tGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAgg                                                                      >  1:902231/1‑148 (MQ=255)
                                                                        tGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAgg                                                                      >  2:2920463/1‑148 (MQ=255)
                                                                         ggAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGa                                                                     <  2:2325648/148‑1 (MQ=255)
                                                                            aTGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATgg                                                                  <  1:2282903/148‑1 (MQ=255)
                                                                             tGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGANTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGt                                                                 >  1:1680479/1‑148 (MQ=255)
                                                                               acacCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGttt                                                               >  2:1086465/1‑148 (MQ=255)
                                                                               acacCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGttt                                                               <  1:2390717/148‑1 (MQ=255)
                                                                                 acCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTc                                                             <  1:55479/148‑1 (MQ=255)
                                                                                 acCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTc                                                             <  1:3083132/148‑1 (MQ=255)
                                                                                  ccAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCa                                                            >  2:2180140/1‑148 (MQ=255)
                                                                                    aaGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCTCGCTGCCAACCAGGATGGTTTTCaaa                                                          <  1:249315/148‑1 (MQ=255)
                                                                                    aaGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCaaa                                                          >  1:1512549/1‑148 (MQ=255)
                                                                                      gCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAAt                                                        >  2:336969/1‑148 (MQ=255)
                                                                                      gCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAAt                                                        >  2:3189294/1‑148 (MQ=255)
                                                                                       caccacCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAAtt                                                       <  1:2809257/148‑1 (MQ=255)
                                                                                        accaccAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTg                                                      <  2:2167302/148‑1 (MQ=255)
                                                                                        accaccAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGACATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTg                                                      <  2:503669/148‑1 (MQ=255)
                                                                                          caccaGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCt                                                    >  2:2553283/1‑148 (MQ=255)
                                                                                          caccaGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCt                                                    <  2:752011/148‑1 (MQ=255)
                                                                                             cagcagTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTa                                                 <  1:2414246/148‑1 (MQ=255)
                                                                                               gcagTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAg                                               <  1:1609653/148‑1 (MQ=255)
                                                                                                cagTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGc                                              >  2:1472203/1‑148 (MQ=255)
                                                                                                 agTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCa                                             <  2:2231387/148‑1 (MQ=255)
                                                                                                 agTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCa                                             <  1:1213724/148‑1 (MQ=255)
                                                                                                   tGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCACCCAGGATGGTTTTCAAAATTGCTTTAAGCAAc                                           <  1:3111469/148‑1 (MQ=255)
                                                                                                   tGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAAc                                           >  1:2497703/1‑148 (MQ=255)
                                                                                                    gcgcATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACt                                          >  2:2916300/1‑148 (MQ=255)
                                                                                                          aGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCgcgc                                    <  2:1208751/148‑1 (MQ=255)
                                                                                                            cccTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCa                                  >  1:1420206/1‑148 (MQ=255)
                                                                                                             ccTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAg                                 >  2:2361858/1‑148 (MQ=255)
                                                                                                              cTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGt                                >  1:3222217/1‑148 (MQ=255)
                                                                                                                aCCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCt                              >  2:2366670/1‑148 (MQ=255)
                                                                                                                 ccAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTg                             <  1:1062127/148‑1 (MQ=255)
                                                                                                                  cAAATCCATCGCATTCCCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGa                            <  1:1545740/148‑1 (MQ=255)
                                                                                                                     aTCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGcc                         >  1:42705/1‑148 (MQ=255)
                                                                                                                      tCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCCGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCg                        >  2:287887/1‑148 (MQ=255)
                                                                                                                      tCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCg                        >  2:1955456/1‑148 (MQ=255)
                                                                                                                       ccATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCga                       >  2:595186/1‑148 (MQ=255)
                                                                                                                           cGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAAc                   >  1:1313941/1‑148 (MQ=255)
                                                                                                                              aTTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATc                >  1:557286/1‑148 (MQ=255)
                                                                                                                              aTTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATc                <  1:178390/148‑1 (MQ=255)
                                                                                                                                tACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCg              >  2:65377/1‑148 (MQ=255)
                                                                                                                                tACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCg              >  2:239912/1‑148 (MQ=255)
                                                                                                                                tACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCg              >  2:130497/1‑148 (MQ=255)
                                                                                                                                    cATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTa          >  1:2985191/1‑148 (MQ=255)
                                                                                                                                    cATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTa          <  1:543671/148‑1 (MQ=255)
                                                                                                                                     aTGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTaa         <  2:352217/148‑1 (MQ=255)
                                                                                                                                      tGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAAt        >  2:2837602/1‑148 (MQ=255)
                                                                                                                                           gTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCg   >  2:1426820/1‑148 (MQ=255)
                                                                                                                                            tAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCgg  >  1:1338969/1‑148 (MQ=255)
                                                                                                                                            tAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCgg  <  2:1153205/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
AATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGG  >  CP009273/1959771‑1960058

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: ‑