Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,959,915 C→A 64.2% 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,9150CA64.2% 163.2 / 310.9 246P178P (CCG→CCT) flhBflagellin export apparatus, substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base C (57/31);  new base A (95/63);  total (152/94)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

Only 100 of 250 total aligned reads displayed.

TCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGG  >  CP009273/1959766‑1960064
                                                                                                                                                     |                                                                                                                                                     
tcaCGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATc                                                                                                                                                       <  2:463787/150‑1 (MQ=255)
tcaCGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAAGCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATc                                                                                                                                                       <  2:4337623/150‑1 (MQ=255)
  aCGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCgg                                                                                                                                                     >  2:606111/1‑150 (MQ=255)
   cGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGga                                                                                                                                                    <  1:3255190/150‑1 (MQ=255)
    gAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGgag                                                                                                                                                   >  1:1916479/1‑150 (MQ=255)
    gAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGgag                                                                                                                                                   >  1:60207/1‑150 (MQ=255)
      atatCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCCGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGAc                                                                                                                                                 >  2:3065879/1‑150 (MQ=255)
        atCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTc                                                                                                                                               >  2:4116008/1‑150 (MQ=255)
         tCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCg                                                                                                                                              <  1:2051631/150‑1 (MQ=255)
         tCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCg                                                                                                                                              <  1:4300738/150‑1 (MQ=255)
          ccTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCgg                                                                                                                                             >  1:4521480/1‑150 (MQ=255)
          ccTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCgg                                                                                                                                             >  1:1956957/1‑150 (MQ=255)
          ccTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCgg                                                                                                                                             <  1:3773014/150‑1 (MQ=255)
             gCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCa                                                                                                                                          >  2:1352298/1‑150 (MQ=255)
              ccGTGACATACGCAGCTTTTTCAGGTGTCTGAATATTTTTAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATATCCCTACCAAATCCATCGCATTAACCATGGCGGTAATAGGAGACTCTTCCAt                                                                                                                                         <  2:2780497/150‑1 (MQ=255)
              ccGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCAt                                                                                                                                         >  1:3391094/1‑150 (MQ=255)
                gTGACCTACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCAGTACCCATGGCGGTAATAGGAGACTCGGCCATCa                                                                                                                                       <  2:3159391/150‑1 (MQ=255)
                 tGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCaa                                                                                                                                      >  1:546955/1‑150 (MQ=255)
                   aCATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAgc                                                                                                                                    >  2:4093058/1‑150 (MQ=255)
                   aCATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAgc                                                                                                                                    >  2:3589526/1‑150 (MQ=255)
                   aCATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAgc                                                                                                                                    >  1:1029906/1‑150 (MQ=255)
                     aTACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAgcgc                                                                                                                                  <  1:96670/150‑1 (MQ=255)
                     aTACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAgcgc                                                                                                                                  <  1:94298/150‑1 (MQ=255)
                        cGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGcatc                                                                                                                               >  2:2685031/1‑150 (MQ=255)
                            gCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCt                                                                                                                           <  1:1363947/150‑1 (MQ=255)
                             cTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTACTAGGAGACTCGGCCATCAAGCGCATCATCTg                                                                                                                          >  2:2346239/1‑150 (MQ=255)
                             cTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTg                                                                                                                          >  2:766907/1‑150 (MQ=255)
                             cTTTTTCACGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAATCACCAACAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCAAGGAGGTAATAGGAGACTCGGCCATCTAGCTCATCATCTg                                                                                                                          >  2:768139/1‑150 (MQ=255)
                              tttttCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGc                                                                                                                         >  1:373926/1‑150 (MQ=255)
                              tttttCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGc                                                                                                                         >  1:2233757/1‑150 (MQ=255)
                                tttCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCgg                                                                                                                       >  1:816701/1‑150 (MQ=255)
                                  tCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGcc                                                                                                                     >  1:2487388/1‑150 (MQ=255)
                                   cAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCa                                                                                                                    >  2:3841825/1‑150 (MQ=255)
                                      gTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTg                                                                                                                 >  1:1056074/1‑150 (MQ=255)
                                        ggCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGAt                                                                                                               <  2:3090952/150‑1 (MQ=255)
                                        ggCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGAt                                                                                                               <  1:4031888/150‑1 (MQ=255)
                                           tGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGCTGcc                                                                                                            >  2:4354985/1‑150 (MQ=255)
                                           tGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGcc                                                                                                            >  2:2860620/1‑150 (MQ=255)
                                                atTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAgaga                                                                                                       <  1:4224567/150‑1 (MQ=255)
                                                 tttGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGaa                                                                                                      >  1:1447619/1‑150 (MQ=255)
                                                   tGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGaaaa                                                                                                    >  2:2025978/1‑150 (MQ=255)
                                                   tGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGaaaa                                                                                                    >  1:1009859/1‑150 (MQ=255)
                                                       aaaaGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAAcc                                                                                                >  1:4483074/1‑150 (MQ=255)
                                                        aaaGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAAccc                                                                                               >  2:3524155/1‑150 (MQ=255)
                                                        aaaGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAAccc                                                                                               >  2:4144405/1‑150 (MQ=255)
                                                        aaaGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAAccc                                                                                               >  2:972481/1‑150 (MQ=255)
                                                        aaaGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAAccc                                                                                               >  2:1282463/1‑150 (MQ=255)
                                                             cGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCa                                                                                          <  1:1266948/150‑1 (MQ=255)
                                                             cGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCa                                                                                          <  2:1919666/150‑1 (MQ=255)
                                                              gTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCAc                                                                                         >  1:2765009/1‑150 (MQ=255)
                                                              gTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCAc                                                                                         >  1:3167653/1‑150 (MQ=255)
                                                              gTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCAc                                                                                         >  1:3166627/1‑150 (MQ=255)
                                                                cAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGc                                                                                       >  1:1891426/1‑150 (MQ=255)
                                                                cAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGc                                                                                       >  2:186826/1‑150 (MQ=255)
                                                                 aaaTCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCt                                                                                      >  1:592542/1‑150 (MQ=255)
                                                                        accaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAAcc                                                                               >  1:1320980/1‑150 (MQ=255)
                                                                         ccaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCa                                                                              <  2:4222419/150‑1 (MQ=255)
                                                                          caTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAg                                                                             <  2:2277314/150‑1 (MQ=255)
                                                                           aTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAgg                                                                            >  1:2951963/1‑150 (MQ=255)
                                                                             tGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGAt                                                                          >  2:3312326/1‑150 (MQ=255)
                                                                              ggAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATg                                                                         <  1:4139445/150‑1 (MQ=255)
                                                                               gAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATgg                                                                        >  2:3124988/1‑150 (MQ=255)
                                                                                aaTGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGt                                                                       >  2:4512863/1‑150 (MQ=255)
                                                                                    acacCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTc                                                                   >  1:1709076/1‑150 (MQ=255)
                                                                                    acacCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTc                                                                   <  2:4149582/150‑1 (MQ=255)
                                                                                        cAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCaaaa                                                               >  1:1584454/1‑150 (MQ=255)
                                                                                         aaGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAAt                                                              <  2:3398598/150‑1 (MQ=255)
                                                                                          aGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAAtt                                                             <  1:2407966/150‑1 (MQ=255)
                                                                                           gCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTg                                                            >  2:665107/1‑150 (MQ=255)
                                                                                              ccaccaGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCtt                                                         <  2:927814/150‑1 (MQ=255)
                                                                                               caccaGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCttt                                                        <  1:2854823/150‑1 (MQ=255)
                                                                                                accaGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTa                                                       >  2:2823583/1‑150 (MQ=255)
                                                                                                 ccaGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTaa                                                      <  1:3281019/150‑1 (MQ=255)
                                                                                                  cagcagTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAg                                                     <  2:1770125/150‑1 (MQ=255)
                                                                                                     cagTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCaa                                                  >  2:101823/1‑150 (MQ=255)
                                                                                                     cagTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCaa                                                  >  2:327043/1‑150 (MQ=255)
                                                                                                       gTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACt                                                >  2:1461001/1‑150 (MQ=255)
                                                                                                         gcgcATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTcc                                              <  2:2599292/150‑1 (MQ=255)
                                                                                                              tAGCCCTACCAAATCCATCGCATTACCCATGGCGGTACTCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGcc                                         <  1:1099143/150‑1 (MQ=255)
                                                                                                                 cccTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGt                                      >  2:1524663/1‑150 (MQ=255)
                                                                                                                  ccTACCAAATCCATCTCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTc                                     <  1:4262906/150‑1 (MQ=255)
                                                                                                                   cTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTcg                                    >  1:3672275/1‑149 (MQ=255)
                                                                                                                   cTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCt                                    <  1:4010508/150‑1 (MQ=255)
                                                                                                                        aaaTCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGcc                               >  1:3940097/1‑150 (MQ=255)
                                                                                                                        aaaTCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGcc                               >  1:3940272/1‑150 (MQ=255)
                                                                                                                         aaTCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCg                              <  1:2214495/150‑1 (MQ=255)
                                                                                                                          aTCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCga                             <  2:4264826/150‑1 (MQ=255)
                                                                                                                          aTCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCga                             <  1:3589526/150‑1 (MQ=255)
                                                                                                                          aTCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCga                             >  1:3509637/1‑150 (MQ=255)
                                                                                                                          aTCCATCGCATGACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCga                             <  1:433387/150‑1 (MQ=255)
                                                                                                                           tCCCTCGCATGACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCgag                            <  2:2934411/150‑1 (MQ=255)
                                                                                                                            ccATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCgaga                           <  1:3065879/150‑1 (MQ=255)
                                                                                                                              aTCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAAc                         <  1:3364910/150‑1 (MQ=255)
                                                                                                                              aTCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAAc                         <  1:67539/150‑1 (MQ=255)
                                                                                                                                   aTTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCg                    >  2:2005521/1‑150 (MQ=255)
                                                                                                                                      aCCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGttt                 <  2:3253307/150‑1 (MQ=255)
                                                                                                                                              cGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCg         >  1:4199585/1‑150 (MQ=255)
                                                                                                                                                 tAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTa      >  1:3608901/1‑150 (MQ=255)
                                                                                                                                                  aaTAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAg     >  2:2124507/1‑150 (MQ=255)
                                                                                                                                                    tCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCg   <  2:3822089/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGG  >  CP009273/1959766‑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: ‑