Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,068,408 T→A 7.6% pseudogene (12/159 nt) yoeH ← pseudogene, CP4‑44 putative prophage remnant;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,068,4080TA 7.6% 500.9 / 24.8 153pseudogene (12/159 nt)yoeHpseudogene, CP4‑44 putative prophage remnant;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base T (62/78);  new base A (6/6);  total (68/85)
Fisher's exact test for biased strand distribution p-value = 7.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

Only 100 of 154 total aligned reads displayed.

GACCTGATACTGATATGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAA  >  NC_000913/2068259‑2068556
                                                                                                                                                     |                                                                                                                                                    
  ccTGATACTGATATGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGt                                                                                                                                                    >  2:1340223/1‑150 (MQ=255)
   ctgatactgataTGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTc                                                                                                                                                   <  1:2493450/150‑1 (MQ=255)
     gatactgataTGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCaa                                                                                                                                                 >  1:2224532/1‑150 (MQ=255)
       tactgataTGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAaga                                                                                                                                               <  1:3750454/150‑1 (MQ=255)
       tactgataTGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAaga                                                                                                                                               <  1:3314384/150‑1 (MQ=255)
         ctgataTGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGt                                                                                                                                             <  1:3933900/150‑1 (MQ=255)
          tgataTGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTg                                                                                                                                            <  2:842809/150‑1 (MQ=255)
          tgataTGAGTCTTACCGCATAGACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTg                                                                                                                                            <  2:2952165/150‑1 (MQ=255)
            atatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCTAGAGTGCg                                                                                                                                          >  1:1133256/1‑150 (MQ=255)
            atatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCg                                                                                                                                          <  1:3977284/150‑1 (MQ=255)
            atatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCg                                                                                                                                          <  1:3977268/150‑1 (MQ=255)
            atatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCg                                                                                                                                          <  1:3830014/150‑1 (MQ=255)
            atatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCg                                                                                                                                          <  1:2911157/150‑1 (MQ=255)
             tatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTGACATAAGAGTGCgg                                                                                                                                         <  1:624078/150‑1 (MQ=255)
               tGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTGACATAAGAGTGCGGTc                                                                                                                                       <  2:221906/150‑1 (MQ=255)
               tGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTGACATAAGAGTGCGGTc                                                                                                                                       <  2:3777463/150‑1 (MQ=255)
                  gTCTTACCGCATATACGGTCAATTTCAGCAATTAATTATATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTGACATAAGAGTGCGGTCCGt                                                                                                                                    <  2:2251866/150‑1 (MQ=39)
                  gTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTGACATAAGAGTGCGGTCCGt                                                                                                                                    <  2:817092/150‑1 (MQ=39)
                  gTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGt                                                                                                                                    <  2:180934/150‑1 (MQ=255)
                    cgtcccGCATCTACGGTCAATTTCAGCAATTAATTACAGTATCCACGCCAAAGTATTTGTCATCACACTGATGGTACCTTCTTGCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTgg                                                                                                                                  <  2:533301/146‑1 (MQ=255)
                    cTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTgg                                                                                                                                  <  2:3793252/150‑1 (MQ=255)
                     ttACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGa                                                                                                                                 >  2:253129/1‑150 (MQ=255)
                     ttACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGa                                                                                                                                 >  2:4000008/1‑150 (MQ=255)
                       aCCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATa                                                                                                                               >  1:2424731/1‑150 (MQ=255)
                          gCATATACGGTCACTTTCAGCAAGTAATTACATTATCCCCGCCAACGTATTAGTCACCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCCTATGTCAAGAGTGCGGTCCGTGGATACaa                                                                                                                            <  2:2616759/150‑1 (MQ=255)
                          gCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCGCTCTTATGTCAAGAGTGCGGTCCGTGGATACaa                                                                                                                            <  2:1196214/150‑1 (MQ=255)
                              ataCGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGATTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCag                                                                                                                        >  2:1315124/1‑150 (MQ=255)
                                 cGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAgaga                                                                                                                     >  2:3608163/1‑150 (MQ=255)
                                 cGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAgaga                                                                                                                     <  1:4083905/150‑1 (MQ=255)
                                    tCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAGGAGTGCGGTCCGTGGATACAACCAGAGACCg                                                                                                                  >  1:2217556/1‑150 (MQ=255)
                                    tCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCg                                                                                                                  >  2:1927071/1‑150 (MQ=255)
                                     cAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGa                                                                                                                 <  2:2858064/150‑1 (MQ=255)
                                       aTTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACt                                                                                                               <  2:575783/150‑1 (MQ=255)
                                         ttCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGa                                                                                                             <  1:967223/150‑1 (MQ=255)
                                         ttCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGa                                                                                                             <  1:3266417/150‑1 (MQ=255)
                                         ttCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGa                                                                                                             <  2:3368943/150‑1 (MQ=255)
                                          tCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGac                                                                                                            <  1:3822911/150‑1 (MQ=255)
                                          tCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAACGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGac                                                                                                            >  2:2582125/1‑150 (MQ=255)
                                           cAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGaca                                                                                                           <  1:185338/150‑1 (MQ=255)
                                            aGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGacac                                                                                                          >  2:4009817/1‑150 (MQ=255)
                                                attaattaCATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGt                                                                                                      <  1:3220493/150‑1 (MQ=255)
                                                attaattaCATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGt                                                                                                      <  1:755946/150‑1 (MQ=255)
                                                    attaCATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCTCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCaga                                                                                                  <  1:3222386/150‑1 (MQ=255)
                                                    attaCATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCaga                                                                                                  <  1:3771136/150‑1 (MQ=255)
                                                    attaCATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCaga                                                                                                  <  1:2552909/150‑1 (MQ=255)
                                                          ttATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAac                                                                                            >  1:1924909/1‑150 (MQ=255)
                                                                  gCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGata                                                                                    >  2:2135487/1‑150 (MQ=255)
                                                                         tATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAg                                                                             <  2:354287/150‑1 (MQ=255)
                                                                         tATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAg                                                                             >  2:362160/1‑150 (MQ=255)
                                                                         gatgtGTCATCACAATGATGGTACCTTCTTTCCGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAg                                                                             <  2:354276/146‑1 (MQ=255)
                                                                               tcatcaCAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTaa                                                                       >  2:1962551/1‑150 (MQ=255)
                                                                                 atcaCAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTaaaa                                                                     >  1:4021258/1‑150 (MQ=255)
                                                                                      aaTGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACACTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGGTCAGTCGTAAAAtatct                                                                >  1:1515984/1‑150 (MQ=255)
                                                                                            ggTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGTATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCaaa                                                          <  1:2200016/150‑1 (MQ=255)
                                                                                              tACCTTCTTTCAGACACAATTTTTTCAACTCCGTTGTCCACGGACCGCACACTTAGGTCAAGAGTGCGGTCCGTGTATACAACCAGAGACCGACTGACACGAGTCAGTGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCaaaaa                                                        <  1:1744080/150‑1 (MQ=255)
                                                                                               aCCTTCTTTCAGACACCATTTTTTAAACTCCTTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACTAGTCAGTGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCaaaaaa                                                       <  1:1357162/150‑1 (MQ=255)
                                                                                               aCCTTCTTTAAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCaaaaaa                                                       <  1:1357164/150‑1 (MQ=255)
                                                                                                 cttcttTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACa                                                     >  2:3194831/1‑150 (MQ=255)
                                                                                                  ttcttTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACAt                                                    >  2:2311568/1‑150 (MQ=255)
                                                                                                  ttcttTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCACAGACCCACTGAAACGAGGCAGAGGAAACGACCGATATGTTAAGTAGTAAAATATCTATCAAAAAACAt                                                    >  1:3688284/1‑150 (MQ=255)
                                                                                                     tttCAGACACCATTTTTTCAACTCCGTTTTCCACGGGCCCCACACTTAAGTCCAGAGTGGGGGCCGGGGGTACAACCAGAGACCCACAGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGAt                                                 >  1:53590/1‑150 (MQ=255)
                                                                                                     tttCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCCGAGGAAACGACGGATATTTTCAGTCGTAAAATATCTATAAAAAAACAAGAt                                                 >  1:2729179/1‑150 (MQ=255)
                                                                                                     tttCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACGCATATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTAGTAAAAAATCAAAAAAAAAACATGAt                                                 >  1:954043/1‑150 (MQ=255)
                                                                                                     tttCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACGCATATGTCAAGAGTGCGGTCCGTGGATAAAACCAGAGACCGACAGAAACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATATATAAAAAAAAATGAt                                                 >  1:53681/1‑150 (MQ=255)
                                                                                                        cAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGTGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTaa                                              <  1:2815992/150‑1 (MQ=255)
                                                                                                         aGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAg                                             <  1:1224201/150‑1 (MQ=255)
                                                                                                         aGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAg                                             <  1:180882/150‑1 (MQ=255)
                                                                                                         aGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAg                                             <  1:371504/150‑1 (MQ=255)
                                                                                                         aGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAg                                             <  2:697652/150‑1 (MQ=255)
                                                                                                         aGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAg                                             <  2:673950/150‑1 (MQ=255)
                                                                                                               cATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCaaa                                       <  1:2307321/150‑1 (MQ=255)
                                                                                                               cATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCaaa                                       <  1:2307108/150‑1 (MQ=255)
                                                                                                               cATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCaaa                                       >  2:1140616/1‑150 (MQ=255)
                                                                                                               cATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCaaa                                       >  2:1865438/1‑150 (MQ=255)
                                                                                                                     ttCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGtt                                 <  2:3415092/150‑1 (MQ=255)
                                                                                                                     ttCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGtt                                 <  1:2377286/150‑1 (MQ=255)
                                                                                                                      tCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGttt                                <  2:3962666/150‑1 (MQ=255)
                                                                                                                       cAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTg                               <  2:3678615/150‑1 (MQ=255)
                                                                                                                       cAACTCCGATTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGACAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTg                               <  2:3678432/150‑1 (MQ=255)
                                                                                                                           tCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGatat                           <  1:3608163/150‑1 (MQ=255)
                                                                                                                            gttgtatCCACGGACCGCACTCTTGACATAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATAtt                          >  2:3900636/7‑150 (MQ=38)
                                                                                                                            ccGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATAtt                          >  2:2602284/1‑150 (MQ=255)
                                                                                                                            ccGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATAtt                          <  1:3612158/150‑1 (MQ=255)
                                                                                                                            ccGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATAAAGGTCAAAAATGTTTGATAtt                          >  2:2738823/1‑150 (MQ=255)
                                                                                                                             cGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATAttt                         >  2:330166/1‑150 (MQ=255)
                                                                                                                                    cACGGACCGCATTCGTATGTAAAGAGTGCGGTCCGTGGATACAACCAGTGCCCGACTGACAATAGTCAGAGGAAACGACGGATATGTTCTGTCGTAAAATATCTATCCAAAAACATGATTAAGGTCAAAAATGTTTGATATTTAGTAttt                  <  1:612908/150‑1 (MQ=255)
                                                                                                                                     aCGGACCGCACTCTTGACATAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTa                 >  1:1214525/1‑150 (MQ=39)
                                                                                                                                     aCGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTa                 <  1:1109755/150‑1 (MQ=255)
                                                                                                                                       ggACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATAGGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATg               <  2:1756598/150‑1 (MQ=255)
                                                                                                                                           cGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGa           >  1:1657321/1‑150 (MQ=255)
                                                                                                                                           cGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGa           >  1:829291/1‑150 (MQ=255)
                                                                                                                                           cGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGa           >  2:916959/1‑150 (MQ=255)
                                                                                                                                            gCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATTTTTGATATTTACAATTTATGAAGAt          <  2:2372008/150‑1 (MQ=255)
                                                                                                                                            gCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGAt          >  1:2727055/1‑150 (MQ=255)
                                                                                                                                            gCACTATTATTTCAAGAGAGCGGTCCGTGTATACAACCAGAGACCTACTGACACGAGTCAGATGAAACGACGGATATGATAAGACGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAAAAAATGATATATACAATATATGAAGAt          >  1:2867875/1‑150 (MQ=255)
                                                                                                                                             cACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATg         >  1:756830/1‑150 (MQ=255)
                                                                                                                                             cACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATg         >  2:879888/1‑150 (MQ=255)
                                                                                                                                                tctTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACa      >  2:811073/1‑150 (MQ=255)
                                                                                                                                                tctTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACa      >  2:2439380/1‑150 (MQ=255)
                                                                                                                                                    aTGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACaataa  <  2:4021258/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GACCTGATACTGATATGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAA  >  NC_000913/2068259‑2068556

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