Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,641,686 G→A 8.0% pseudogene (208/279 nt) ydfX ← pseudogene, Qin prophage;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,641,6860GA 8.0% 1010.6 / 53.9 299pseudogene (208/279 nt)ydfXpseudogene, Qin prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (129/145);  new base A (8/16);  total (137/162)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 302 total aligned reads displayed.

TTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCTGTGTCATTC  >  CP009273/1641537‑1641835
                                                                                                                                                     |                                                                                                                                                     
  aaCATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGtt                                                                                                                                                     <  1:2380094/150‑1 (MQ=255)
    cATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCa                                                                                                                                                   <  2:2785800/150‑1 (MQ=255)
    cATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCa                                                                                                                                                   <  1:4118427/150‑1 (MQ=255)
     aTAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCaa                                                                                                                                                  <  2:4878202/150‑1 (MQ=255)
       aGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAat                                                                                                                                                >  2:133808/1‑150 (MQ=255)
        gTATTTTGAGCGTAAATCAATACGGGCAATATCCACTCGTACATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAata                                                                                                                                               >  2:2216110/1‑150 (MQ=255)
             ttGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTc                                                                                                                                          <  1:1156365/150‑1 (MQ=255)
              tGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCATTCAAATATCTTCa                                                                                                                                         <  1:4461190/150‑1 (MQ=255)
              tGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCAGAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCa                                                                                                                                         >  1:2235891/1‑150 (MQ=255)
              cgaccGTAAATCACTACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCa                                                                                                                                         <  1:4864777/146‑1 (MQ=255)
                     aaaTCAATACGGGCATTATCCACTCGGTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGAtt                                                                                                                                  >  2:1015094/1‑150 (MQ=255)
                      aaTCAATACGGGCATTATCCACTCGTTCATGGACTCTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTa                                                                                                                                 <  2:3899415/150‑1 (MQ=255)
                      aaTCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTa                                                                                                                                 <  1:1851267/150‑1 (MQ=255)
                      aaTCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCAGCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTa                                                                                                                                 <  2:585655/150‑1 (MQ=255)
                       aTCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTaa                                                                                                                                <  2:1861301/150‑1 (MQ=255)
                         cAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACg                                                                                                                              <  1:3684320/150‑1 (MQ=255)
                         cAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACg                                                                                                                              >  1:3514149/1‑150 (MQ=255)
                         cAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACg                                                                                                                              >  2:2429335/1‑150 (MQ=255)
                          aaTACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCATTCACATATCTTCACTCGATTAACGc                                                                                                                             <  2:1972864/150‑1 (MQ=255)
                          aaTACGGGCATTATCCACTCGTTCATGGACACTCATTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCATTCAAATATCTTCACTCGATTAACGc                                                                                                                             <  2:730344/150‑1 (MQ=255)
                             aCGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCaa                                                                                                                          >  1:267378/1‑150 (MQ=255)
                              cGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAg                                                                                                                         <  1:4161670/150‑1 (MQ=255)
                                  cATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGcgct                                                                                                                     <  2:4652433/150‑1 (MQ=255)
                                  cATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGcgct                                                                                                                     >  1:4664538/1‑150 (MQ=255)
                                    ttATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGcgctcg                                                                                                                   <  1:4019194/150‑1 (MQ=255)
                                    ttATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACACGATTAACGCCAAGcgctcg                                                                                                                   >  2:3593723/1‑150 (MQ=255)
                                     tATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGcgctcgc                                                                                                                  >  1:1686819/1‑150 (MQ=255)
                                      aTCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGcgctcgct                                                                                                                 >  2:4395585/1‑150 (MQ=255)
                                      aTCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGcgctcgct                                                                                                                 >  1:1632943/1‑150 (MQ=255)
                                          aCTCGGTCATGGACACTCTTATCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATc                                                                                                             >  2:3130866/1‑150 (MQ=255)
                                             cGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAAt                                                                                                          >  2:340592/1‑150 (MQ=255)
                                             cGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAAt                                                                                                          >  1:4787883/1‑150 (MQ=255)
                                               ttCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATcg                                                                                                        <  1:2681124/150‑1 (MQ=255)
                                               ttCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATcg                                                                                                        <  1:355728/150‑1 (MQ=255)
                                                 cATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATcgcg                                                                                                      >  2:4895132/1‑150 (MQ=255)
                                                 cATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATcgcg                                                                                                      >  1:748365/1‑150 (MQ=255)
                                                      acacTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCAcc                                                                                                 >  1:4616431/1‑150 (MQ=255)
                                                       cacTCTTTTCATCCTCCAGTGGGAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCt                                                                                                <  2:3189935/150‑1 (MQ=255)
                                                       cacTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCAATCGATTAACGACAAGCGCTCGCTGATCAATCGCGCCACCt                                                                                                >  2:3155025/1‑150 (MQ=255)
                                                         ctctTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTc                                                                                              <  1:4387953/150‑1 (MQ=255)
                                                          tctTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCATTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCa                                                                                             <  1:2122481/150‑1 (MQ=255)
                                                           ctTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCat                                                                                            <  2:644255/150‑1 (MQ=255)
                                                               tcatcaTCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATAtgt                                                                                        >  1:1489341/1‑150 (MQ=255)
                                                                catcatCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATAtgtg                                                                                       >  2:1146023/1‑150 (MQ=255)
                                                                 atcatcCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGa                                                                                      >  1:3611271/1‑150 (MQ=255)
                                                                 atcatcCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGa                                                                                      <  1:5072148/150‑1 (MQ=255)
                                                                  tcatcCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGac                                                                                     >  1:961471/1‑150 (MQ=255)
                                                                   catcCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGaca                                                                                    <  1:4132174/150‑1 (MQ=255)
                                                                   catcCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGaca                                                                                    >  2:1263995/1‑150 (MQ=255)
                                                                    atcCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGacac                                                                                   >  1:3496263/1‑150 (MQ=255)
                                                                     tcCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACg                                                                                  <  1:3868220/150‑1 (MQ=255)
                                                                      ccAGGGGGAGCCTCCAGCAGTGCAGGCACTAGCGCCGTGCACTTTTCACGCGTCTCCCTGGGGTCGATAGCCGGCCAGCGTTCAAATAGCGTCACGCGAGGACCGCCAAGCGCGCGCGGATCAATCGCGCCACCGTCAGAGGGGACACGc                                                                                 >  2:5169575/1‑150 (MQ=255)
                                                                       cAGTGGTAGCCTCAATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCt                                                                                <  2:4381489/150‑1 (MQ=255)
                                                                         gTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGa                                                                              <  1:225551/150‑1 (MQ=255)
                                                                         gTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGa                                                                              >  1:553428/1‑150 (MQ=255)
                                                                         atGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTGCATATGTGACACGCTGa                                                                              >  2:4027203/2‑150 (MQ=255)
                                                                          tGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGaa                                                                             <  2:5145666/150‑1 (MQ=255)
                                                                          tGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGaa                                                                             <  2:2555449/150‑1 (MQ=255)
                                                                          tGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCGTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGaa                                                                             <  2:4809969/150‑1 (MQ=255)
                                                                          tGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCATTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGaa                                                                             <  1:804974/150‑1 (MQ=255)
                                                                              aGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCATTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATc                                                                         <  2:4072007/150‑1 (MQ=255)
                                                                                    cATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGtt                                                                   <  1:3221600/150‑1 (MQ=255)
                                                                                     aTAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTc                                                                  <  2:1176017/150‑1 (MQ=255)
                                                                                     aTAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTc                                                                  >  1:228153/1‑150 (MQ=255)
                                                                                     aTAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTc                                                                  <  2:2015554/150‑1 (MQ=255)
                                                                                       aaTTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCATTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCAGATGTGACACGCTGAACATCGATGTTCgg                                                                >  2:3883534/1‑150 (MQ=255)
                                                                                         ttGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAACTATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGgc                                                              <  1:2602695/150‑1 (MQ=255)
                                                                                          tGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACCTCGCTGTTCGGgcg                                                             >  2:4505329/1‑150 (MQ=255)
                                                                                             aGGCACTAGCGCCGTGAACTTTTCACGATTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGctc                                                          <  2:2065560/150‑1 (MQ=255)
                                                                                                cACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCttt                                                       <  1:100648/150‑1 (MQ=255)
                                                                                                 aCTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTc                                                      >  1:620588/1‑150 (MQ=255)
                                                                                                   tAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCaa                                                    >  2:2487361/1‑150 (MQ=255)
                                                                                                    aGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCaaa                                                   >  1:3082508/1‑150 (MQ=255)
                                                                                                       gcCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGcc                                                >  1:198366/1‑150 (MQ=255)
                                                                                                        ccGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGccc                                               <  1:2764816/150‑1 (MQ=255)
                                                                                                        ccGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGccc                                               <  1:1352998/150‑1 (MQ=255)
                                                                                                         cGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCa                                              <  1:4144305/150‑1 (MQ=255)
                                                                                                           tGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGa                                            >  2:1201687/1‑150 (MQ=255)
                                                                                                               cTTGTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACTCGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGc                                        <  2:2657835/150‑1 (MQ=255)
                                                                                                                   tCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCa                                    <  2:862197/150‑1 (MQ=255)
                                                                                                                     aCGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGt                                  <  2:2891161/150‑1 (MQ=255)
                                                                                                                     aCGCTTATCCCTGGTGTCGATAGCCTTCCAGCATTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGt                                  <  2:4292085/150‑1 (MQ=255)
                                                                                                                           aTCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTaa                            >  2:4258380/1‑150 (MQ=255)
                                                                                                                              ccTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAatat                         <  1:3751626/150‑1 (MQ=255)
                                                                                                                                tGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATAtcg                       <  1:4165500/150‑1 (MQ=255)
                                                                                                                                tGGTGTCGATAGCCTTCCAGCATTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATAtcg                       >  2:5243504/1‑150 (MQ=255)
                                                                                                                                   tgtCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATAtcgtcg                    >  2:2495303/1‑150 (MQ=255)
                                                                                                                                    gtCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGc                   >  1:2216483/1‑150 (MQ=255)
                                                                                                                                    gtCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGc                   >  1:452605/1‑150 (MQ=255)
                                                                                                                                     tCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGcc                  <  2:2628733/150‑1 (MQ=255)
                                                                                                                                     tCGATAGCCTTCCAGCATTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGcc                  >  1:3350944/1‑150 (MQ=255)
                                                                                                                                      cGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCt                 >  2:4203938/1‑150 (MQ=255)
                                                                                                                                      cGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCt                 >  1:2177605/1‑150 (MQ=255)
                                                                                                                                      cGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCt                 <  1:3009456/150‑1 (MQ=255)
                                                                                                                                        aTAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGc               >  1:4301609/1‑150 (MQ=255)
                                                                                                                                          agccgaccAGCGATCAAATATCTTCACTCGATTAACGCCACGCGCTCGCTGATCAATCGCGCCACCTTCAAATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGAGATAAATATCGTCGCCAGCTc             >  2:5118238/7‑150 (MQ=255)
                                                                                                                                            ccTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCt           <  2:328634/150‑1 (MQ=255)
                                                                                                                                              ttCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCtgt         >  1:2866032/1‑150 (MQ=255)
                                                                                                                                                 cAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCTGTGTc      <  1:5244570/150‑1 (MQ=255)
                                                                                                                                                    cGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCTGTGTCAtt   <  1:2643355/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTAACATAGTATTTTGAGCGTAAATCAATACGGGCATTATCCACTCGTTCATGGACACTCTTTTCATCATCCAGTGGTAGCCTCCATAATTGCAGGCACTAGCGCCGTGAACTTTTCACGCTTATCCCTGGTGTCGATAGCCTTCCAGCGTTCAAATATCTTCACTCGATTAACGCCAAGCGCTCGCTGATCAATCGCGCCACCTTCATATGTGACACGCTGAACATCGATGTTCGGGCGCTCTTTCAAAGCCCAGAATGCTTCAGTGATTAATATCGTCGCCTGCTCCTGTGTCATTC  >  CP009273/1641537‑1641835

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