Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,959,915 C→A 47.1% 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,9150CA47.1% 170.6 / inf 299P178P (CCG→CCTflhBflagellin export apparatus, substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base C (86/72);  new base A (65/76);  total (151/148)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 301 total aligned reads displayed.

ACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGG  >  CP009273/1959768‑1960064
                                                                                                                                                   |                                                                                                                                                     
aCGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCgg                                                                                                                                                     >  2:1660419/1‑150 (MQ=255)
aCGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCgg                                                                                                                                                     <  2:4715647/150‑1 (MQ=255)
    atatCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGAc                                                                                                                                                 >  2:3075349/1‑150 (MQ=255)
    atatCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGAc                                                                                                                                                 >  1:716813/1‑150 (MQ=255)
    atatCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGAc                                                                                                                                                 <  2:119393/150‑1 (MQ=255)
    atatCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGAc                                                                                                                                                 <  2:118299/150‑1 (MQ=255)
     tatCCTGCCGCGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACt                                                                                                                                                >  2:4641614/1‑150 (MQ=255)
       tCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCg                                                                                                                                              <  2:2621128/150‑1 (MQ=255)
          tGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGcc                                                                                                                                           >  2:4121377/1‑150 (MQ=255)
           gCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCa                                                                                                                                          >  2:3013305/1‑150 (MQ=255)
           gCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCa                                                                                                                                          <  1:2754703/150‑1 (MQ=255)
            ccGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCAt                                                                                                                                         <  2:2403222/150‑1 (MQ=255)
              gTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCa                                                                                                                                       >  2:1130767/1‑150 (MQ=255)
                gACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAg                                                                                                                                     <  1:1514505/150‑1 (MQ=255)
                gACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAg                                                                                                                                     <  2:4667779/150‑1 (MQ=255)
                 aCATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAgc                                                                                                                                    >  1:2302333/1‑150 (MQ=255)
                   aTACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAgcgc                                                                                                                                  <  2:1855840/150‑1 (MQ=255)
                     aCGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGcat                                                                                                                                <  2:1236964/150‑1 (MQ=255)
                      cGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGcatc                                                                                                                               <  1:4601344/150‑1 (MQ=255)
                      cGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGcatc                                                                                                                               >  2:2282118/1‑150 (MQ=255)
                             ttttCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCg                                                                                                                        >  2:1927315/1‑150 (MQ=255)
                                tCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGcc                                                                                                                     <  1:794952/150‑1 (MQ=255)
                                tCAGGTGGCTGAATATTTGGAAAAAGAAGTCAAATCCCACCATTGGAATGTCACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGcc                                                                                                                     <  2:3456685/150‑1 (MQ=255)
                                 cAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCa                                                                                                                    >  2:452510/1‑150 (MQ=255)
                                      ggCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGAt                                                                                                               >  1:1150288/1‑150 (MQ=255)
                                         tGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGcc                                                                                                            >  2:3189809/1‑150 (MQ=255)
                                          gAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCa                                                                                                           >  2:3599671/1‑150 (MQ=255)
                                          gAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCa                                                                                                           <  1:4140108/150‑1 (MQ=255)
                                          gAATATTTGGAAAAAGACGTCAAATCACACCATTGGAATGACACCAAGCACCAACAGCAGTGCTAATAGAAATACCAAATCAATCGCATTAAACATGGCGGTAATATGAGAATCGGCAATCAAGAGCATCATCTGCTGCCAGTGAttcaa                                                                                                           >  1:2585258/1‑146 (MQ=255)
                                           aaTATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCag                                                                                                          <  1:281493/150‑1 (MQ=255)
                                            atatTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCaga                                                                                                         >  1:1124213/1‑150 (MQ=255)
                                            atatTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCaga                                                                                                         <  2:1568676/150‑1 (MQ=255)
                                             tatTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCagag                                                                                                        <  1:1005550/150‑1 (MQ=255)
                                              atTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAgaga                                                                                                       <  1:1806808/150‑1 (MQ=255)
                                                ttGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGaaa                                                                                                     >  2:2103070/1‑150 (MQ=255)
                                                 tGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGaaaa                                                                                                    >  2:3258685/1‑150 (MQ=255)
                                                  ggAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGaaaaa                                                                                                   >  1:2328304/1‑150 (MQ=255)
                                                        aGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCg                                                                                             >  1:266616/1‑150 (MQ=255)
                                                         gACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGt                                                                                            <  1:54072/150‑1 (MQ=255)
                                                          aCGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTc                                                                                           <  2:4409426/150‑1 (MQ=255)
                                                           cGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCa                                                                                          >  1:2206766/1‑150 (MQ=255)
                                                               aaaTCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCt                                                                                      >  2:1398725/1‑150 (MQ=255)
                                                               aaaTCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCt                                                                                      >  2:3128907/1‑150 (MQ=255)
                                                               aaaTCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCt                                                                                      >  1:3629852/1‑150 (MQ=255)
                                                               aaaTCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCt                                                                                      >  2:3756296/1‑150 (MQ=255)
                                                                aaTCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTg                                                                                     <  1:3605923/150‑1 (MQ=255)
                                                                  tCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGcc                                                                                   >  2:421716/1‑150 (MQ=255)
                                                                  tCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGcc                                                                                   >  1:1639761/1‑150 (MQ=255)
                                                                   cccACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCa                                                                                  <  2:497391/150‑1 (MQ=255)
                                                                   cccACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCa                                                                                  <  1:3095480/150‑1 (MQ=255)
                                                                     caccaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAAc                                                                                <  2:1430169/150‑1 (MQ=255)
                                                                      accaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAAcc                                                                               >  2:778236/1‑150 (MQ=255)
                                                                       ccaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCa                                                                              <  1:838239/150‑1 (MQ=255)
                                                                        caTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAg                                                                             <  1:1769082/150‑1 (MQ=255)
                                                                         aTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAgg                                                                            >  1:919608/1‑150 (MQ=255)
                                                                         aTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAgg                                                                            >  1:3313521/1‑150 (MQ=255)
                                                                         aTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAgg                                                                            >  1:4323642/1‑150 (MQ=255)
                                                                          ttGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGa                                                                           >  1:3354323/1‑150 (MQ=255)
                                                                            ggAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATg                                                                         <  1:4547608/150‑1 (MQ=255)
                                                                            ggAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATg                                                                         <  2:4688423/150‑1 (MQ=255)
                                                                             gAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATgg                                                                        >  1:302672/1‑150 (MQ=255)
                                                                             gAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATgg                                                                        <  1:408049/150‑1 (MQ=255)
                                                                               aTGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGtt                                                                      <  2:60584/150‑1 (MQ=255)
                                                                                 gACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGtttt                                                                    <  1:2109369/150‑1 (MQ=255)
                                                                                            ccaccaGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCtt                                                         <  2:252091/150‑1 (MQ=255)
                                                                                                cagcagTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAg                                                     <  1:3089866/150‑1 (MQ=255)
                                                                                                cagcagTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAg                                                     <  1:2840771/150‑1 (MQ=255)
                                                                                                cagcagTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAg                                                     <  1:642751/150‑1 (MQ=255)
                                                                                                 agcagTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGc                                                    >  2:880114/1‑150 (MQ=255)
                                                                                                 agcagTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGc                                                    <  1:2219461/150‑1 (MQ=255)
                                                                                                 agcagTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGc                                                    >  2:2861223/1‑150 (MQ=255)
                                                                                                   cagTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCaa                                                  <  1:1635822/150‑1 (MQ=255)
                                                                                                         gcATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCcgc                                            >  2:678501/1‑150 (MQ=255)
                                                                                                             aGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCa                                        >  2:1003035/1‑150 (MQ=255)
                                                                                                                ccTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTc                                     >  1:4653327/1‑150 (MQ=255)
                                                                                                                    ccAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAg                                 >  1:376193/1‑150 (MQ=255)
                                                                                                                         tCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCgag                            >  2:4542149/1‑150 (MQ=255)
                                                                                                                           cATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGaa                          >  1:3956545/1‑150 (MQ=255)
                                                                                                                           cATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGaa                          <  1:1953995/150‑1 (MQ=255)
                                                                                                                             tCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACa                        <  1:2948240/150‑1 (MQ=255)
                                                                                                                                cATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATcc                     <  2:485133/150‑1 (MQ=255)
                                                                                                                                cATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATcc                     >  1:3465452/1‑150 (MQ=255)
                                                                                                                                cATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGATAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATcc                     <  1:2456689/150‑1 (MQ=255)
                                                                                                                                 aTTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCg                    <  2:1340996/150‑1 (MQ=255)
                                                                                                                                 aTTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCg                    <  2:3387776/150‑1 (MQ=255)
                                                                                                                                  ttACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGt                   >  2:4626358/1‑150 (MQ=255)
                                                                                                                                  ttACCCAGGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGt                   <  1:1230910/150‑1 (MQ=255)
                                                                                                                                   tACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCTCGCCAGTCTGAGCCGAGAACATCCGtt                  <  2:692987/150‑1 (MQ=255)
                                                                                                                                   tACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGtt                  >  1:3062125/1‑150 (MQ=255)
                                                                                                                                       cATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAAt              <  1:3387486/150‑1 (MQ=255)
                                                                                                                                        aTGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATg             >  1:1447475/1‑150 (MQ=255)
                                                                                                                                        aTGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATg             <  2:2772335/150‑1 (MQ=255)
                                                                                                                                          ggCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGcc           >  2:666835/1‑150 (MQ=255)
                                                                                                                                           gCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGccc          <  2:2206766/150‑1 (MQ=255)
                                                                                                                                              gTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGt       <  1:3189809/150‑1 (MQ=255)
                                                                                                                                              gTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGt       >  2:2332621/1‑150 (MQ=255)
                                                                                                                                               tAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTa      >  2:3677979/1‑150 (MQ=255)
                                                                                                                                                aaTCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAg     >  1:1082192/1‑150 (MQ=255)
                                                                                                                                                  tCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCg   >  1:4569625/1‑150 (MQ=255)
                                                                                                                                                   cGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCgg  >  1:3698240/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGG  >  CP009273/1959768‑1960064

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: