Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,959,915 C→A 57.4% 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,9150CA57.4% 31.2 / inf 235P178P (CCG→CCT) flhBflagellin export apparatus, substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base C (61/39);  new base A (75/60);  total (136/99)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

Only 100 of 236 total aligned reads displayed.

CACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGG  >  CP009273/1959767‑1960064
                                                                                                                                                    |                                                                                                                                                     
   gAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGgag                                                                                                                                                   >  1:2184845/1‑150 (MQ=255)
   gAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGgag                                                                                                                                                   <  1:2036391/150‑1 (MQ=255)
    aaTATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGgaga                                                                                                                                                  >  1:4138155/1‑150 (MQ=255)
     atatCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGATAc                                                                                                                                                 >  2:4175561/1‑150 (MQ=255)
        tCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCg                                                                                                                                              >  1:740348/1‑150 (MQ=255)
           tGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGcc                                                                                                                                           >  1:2664898/1‑150 (MQ=255)
           tGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGcc                                                                                                                                           <  1:2743129/150‑1 (MQ=255)
               gTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCa                                                                                                                                       >  2:1307261/1‑150 (MQ=255)
                tGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCaa                                                                                                                                      >  1:1173956/1‑150 (MQ=255)
                tGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCaa                                                                                                                                      <  1:1447601/150‑1 (MQ=255)
                    aTACGCAGCTTTTTCGGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAgcgc                                                                                                                                  >  2:2127034/1‑150 (MQ=255)
                           gCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCt                                                                                                                           >  1:1724632/1‑150 (MQ=255)
                           gCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCt                                                                                                                           <  1:3599485/150‑1 (MQ=255)
                            cTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTg                                                                                                                          >  2:331034/1‑150 (MQ=255)
                             tttttCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGc                                                                                                                         <  2:1173843/150‑1 (MQ=255)
                              ttttCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCg                                                                                                                        >  2:2536117/1‑150 (MQ=255)
                               tttCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCgg                                                                                                                       >  1:1555344/1‑150 (MQ=255)
                                  cAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCa                                                                                                                    >  2:3975531/1‑150 (MQ=255)
                                  cAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCa                                                                                                                    >  2:4098155/1‑150 (MQ=255)
                                       ggCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGAt                                                                                                               >  2:932472/1‑150 (MQ=255)
                                       ggCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGAt                                                                                                               >  2:588634/1‑150 (MQ=255)
                                       ggCTCAAAATTTGGAAAAAGACGTCAAATCCCACCAATTGAAAGACAACAAGCACCACCAGCAGTGAGCATAGCCCTACCAAAACCATCGCATTACCCATGGAGGGAATAGGAGACTCGGCCATCAATCGCATCATCTGCGGCAAGTTAt                                                                                                               >  2:588713/1‑150 (MQ=255)
                                         cTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGc                                                                                                             >  1:3974579/1‑150 (MQ=255)
                                           gAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCa                                                                                                           >  2:2011194/1‑150 (MQ=255)
                                            aaTATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCag                                                                                                          >  2:2332758/1‑150 (MQ=255)
                                               atTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAgaga                                                                                                       <  2:2925636/150‑1 (MQ=255)
                                                 ttGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGaaa                                                                                                     <  2:2902405/150‑1 (MQ=255)
                                                  tGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGaaaa                                                                                                    <  2:3390734/150‑1 (MQ=255)
                                                    gAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGaaaaaa                                                                                                  >  1:4045071/1‑150 (MQ=255)
                                                     aaaaaGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAAc                                                                                                 <  1:2226298/150‑1 (MQ=255)
                                                     aaaaaGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAAc                                                                                                 >  2:3183918/1‑150 (MQ=255)
                                                        aaGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAAcccc                                                                                              >  2:1279544/1‑150 (MQ=255)
                                                         aGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCg                                                                                             >  1:314687/1‑150 (MQ=255)
                                                          gACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGt                                                                                            >  1:2773455/1‑150 (MQ=255)
                                                          gACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGt                                                                                            >  2:107528/1‑150 (MQ=255)
                                                           aCGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTc                                                                                           <  1:3652010/150‑1 (MQ=255)
                                                           aCGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAGAAACCCCGTc                                                                                           >  1:1738968/1‑150 (MQ=255)
                                                               cAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGAATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGTAGACTCGGCCATCAATCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGc                                                                                       >  1:3706565/1‑150 (MQ=255)
                                                                aaaTCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCt                                                                                      <  2:1601933/150‑1 (MQ=255)
                                                                aaaTCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCt                                                                                      <  2:4054696/150‑1 (MQ=255)
                                                                aaaTCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCt                                                                                      >  2:2396800/1‑150 (MQ=255)
                                                                aaaTCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCt                                                                                      >  2:1825829/1‑150 (MQ=255)
                                                                  aTCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGc                                                                                    >  1:3450634/1‑150 (MQ=255)
                                                                   tCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGcc                                                                                   >  1:2757178/1‑150 (MQ=255)
                                                                    cccACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCa                                                                                  <  2:1907535/150‑1 (MQ=255)
                                                                     ccaccaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCaa                                                                                 >  1:3512287/1‑150 (MQ=255)
                                                                      caccaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCCGAGAAAAAACCCCGTCACGCTGCCAAc                                                                                >  2:135313/1‑150 (MQ=255)
                                                                      caccaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAAc                                                                                >  2:1571080/1‑150 (MQ=255)
                                                                       accaTTGTAATGACACAAAGCACCACCAGCAGTGCGAATAGCCCTACCATATCCATCGCATTTCCCATGTCGGTAATCGGTGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAAcc                                                                               <  1:80718/150‑1 (MQ=255)
                                                                       accaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAAcc                                                                               <  2:665990/150‑1 (MQ=255)
                                                                       accaTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAAcc                                                                               <  2:3918625/150‑1 (MQ=255)
                                                                         caTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAg                                                                             >  1:2421105/1‑150 (MQ=255)
                                                                          aTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAgg                                                                            >  2:3848809/1‑150 (MQ=255)
                                                                          aTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAgg                                                                            <  1:1127433/150‑1 (MQ=255)
                                                                          aTTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAgg                                                                            >  2:1285785/1‑150 (MQ=255)
                                                                               aaTGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGt                                                                       <  1:328240/150‑1 (MQ=255)
                                                                               aaTGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGt                                                                       >  1:2876601/1‑150 (MQ=255)
                                                                                aTGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGtt                                                                      <  2:3203711/150‑1 (MQ=255)
                                                                                aTGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGtt                                                                      >  1:4018464/1‑150 (MQ=255)
                                                                                  gACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGtttt                                                                    <  2:3948484/150‑1 (MQ=255)
                                                                                    cacCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCa                                                                  >  1:2333715/1‑150 (MQ=255)
                                                                                     acCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCaa                                                                 <  2:268876/150‑1 (MQ=255)
                                                                                        aaGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAAt                                                              <  1:263117/150‑1 (MQ=255)
                                                                                         aGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTNACGCTGCCAACCAGGATGGTTTTCAAAAtt                                                             <  1:3851773/150‑1 (MQ=255)
                                                                                         aGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAAtt                                                             <  2:1180152/150‑1 (MQ=255)
                                                                                            accaccAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCt                                                          <  2:2573944/150‑1 (MQ=255)
                                                                                                ccaGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTaa                                                      <  1:1896779/150‑1 (MQ=255)
                                                                                                    cagTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCaa                                                  <  2:3240482/150‑1 (MQ=255)
                                                                                                    cagTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCaa                                                  >  2:83629/1‑150 (MQ=255)
                                                                                                     agTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAAc                                                 >  2:2768484/1‑150 (MQ=255)
                                                                                                     agTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAAc                                                 >  1:58612/1‑150 (MQ=255)
                                                                                                      gTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACt                                                <  1:1552276/150‑1 (MQ=255)
                                                                                                       tGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTc                                               >  2:735593/1‑150 (MQ=255)
                                                                                                       tGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTc                                               <  1:4025453/150‑1 (MQ=255)
                                                                                                         cgcATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCcg                                             >  2:2659409/1‑150 (MQ=255)
                                                                                                         cgcATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCcg                                             <  1:2642174/150‑1 (MQ=255)
                                                                                                         cgcATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCcg                                             >  2:185172/1‑150 (MQ=255)
                                                                                                           cATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCcgcg                                           >  2:4104309/1‑150 (MQ=255)
                                                                                                             tAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGcc                                         >  1:3823034/1‑150 (MQ=255)
                                                                                                               gCCCTACCAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAg                                       <  1:3186830/150‑1 (MQ=255)
                                                                                                                     ccAAATCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAg                                 >  1:239359/1‑150 (MQ=255)
                                                                                                                       aaaTCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGcc                               >  2:2590454/1‑150 (MQ=255)
                                                                                                                        aaTCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCg                              >  1:2422637/1‑150 (MQ=255)
                                                                                                                        aaTCCATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCg                              >  2:2924438/1‑150 (MQ=255)
                                                                                                                            cATCGCATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGaa                          <  2:869568/150‑1 (MQ=255)
                                                                                                                               cGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACAt                       >  2:3835790/1‑150 (MQ=255)
                                                                                                                                gCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATc                      <  2:2890129/150‑1 (MQ=255)
                                                                                                                                 cATTACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATcc                     <  1:3054507/150‑1 (MQ=255)
                                                                                                                                   ttACCCATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGt                   >  1:4107470/1‑150 (MQ=255)
                                                                                                                                        cATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAAt              >  1:210523/1‑150 (MQ=255)
                                                                                                                                        cATGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAAt              <  1:3092324/150‑1 (MQ=255)
                                                                                                                                         aTGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATg             <  2:2299827/150‑1 (MQ=255)
                                                                                                                                         aTGGCGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAAta             <  1:931256/150‑2 (MQ=255)
                                                                                                                                             cGGTAATAGGAGACTCGGCCATCAATCGAATCATCTGCGGCCAGTGATGCAAGAGAAAAAACCCCGTCACGCTGCCAACCAGGTTGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCg         <  1:745681/150‑1 (MQ=255)
                                                                                                                                             cGGTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCg         <  1:992868/150‑1 (MQ=255)
                                                                                                                                               gTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGt       <  1:131995/150‑1 (MQ=255)
                                                                                                                                               gTAATAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGt       <  1:3568106/150‑1 (MQ=255)
                                                                                                                                                 aaTCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAg     >  2:3639411/1‑150 (MQ=255)
                                                                                                                                                 aaTAGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAg     >  1:531223/1‑150 (MQ=255)
                                                                                                                                                    cGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCgg  <  2:1381764/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CACGAATATCCTGCCGTGACATACGCAGCTTTTTCAGGTGGCTGAATATTTGGAAAAAGACGTCAAATCCCACCATTGGAATGACACCAAGCACCACCAGCAGTGCGCATAGCCCTACCAAATCCATCGCATTACCCATGGCGGTAATCGGAGACTCGGCCATCAAGCGCATCATCTGCGGCCAGTGATGCCAGAGAAAAAACCCCGTCACGCTGCCAACCAGGATGGTTTTCAAAATTGCTTTAAGCAACTCCGCGCCAGTCTGAGCCGAGAACATCCGTTTAATGCCCGGTAGCGG  >  CP009273/1959767‑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: ‑