Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,068,411 C→T 5.0% pseudogene (9/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,4110CT 5.0% 335.4 / 4.6 122pseudogene (9/159 nt)yoeHpseudogene, CP4‑44 putative prophage remnant;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (56/57);  new base T (3/3);  total (61/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.92e-01

Only 100 of 128 total aligned reads displayed.

CCTGATACTGATATGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAATTAT  >  NC_000913/2068261‑2068560
                                                                                                                                                      |                                                                                                                                                     
ccTGATACTGATATGAGTCATACCGCATATACGGTCAATTTCAGCAATTAATAACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTc                                                                                                                                                       <  1:128804/150‑1 (MQ=255)
 ctgatactgataTGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTc                                                                                                                                                       <  1:1760966/150‑1 (MQ=255)
 ctgatactgataTGAGTATTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTc                                                                                                                                                       <  1:1426709/150‑1 (MQ=255)
  tgatactgataTGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCa                                                                                                                                                      >  2:1823568/1‑150 (MQ=255)
       ctgataTGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGt                                                                                                                                                 >  1:17920/1‑150 (MQ=255)
       ctgataTGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGt                                                                                                                                                 >  1:1635297/1‑150 (MQ=255)
        tgataTGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTg                                                                                                                                                <  1:1688044/150‑1 (MQ=255)
         gataTGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGCCCGCACTCTTATGTCAAGAGTGc                                                                                                                                               >  2:408650/1‑150 (MQ=255)
           tatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCgg                                                                                                                                             <  1:1224896/150‑1 (MQ=255)
           tatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCgg                                                                                                                                             <  2:783846/150‑1 (MQ=255)
           tatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCCCACTCTTATGTCAAGAGTGCgg                                                                                                                                             >  1:1179591/1‑150 (MQ=255)
            atGAGTCTTACCGCATATACTGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTGACATAAGAGTGCGGt                                                                                                                                            <  1:192430/150‑1 (MQ=255)
            atGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTGACATAAGAGTGCGGt                                                                                                                                            <  2:799845/150‑1 (MQ=255)
             tGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTc                                                                                                                                           >  1:323497/1‑150 (MQ=255)
             tGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTc                                                                                                                                           >  1:220586/1‑150 (MQ=255)
             tGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTc                                                                                                                                           >  1:992564/1‑150 (MQ=255)
             tGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTACGTTTTCCTCGGACAGCACTCTTTTGTCAAGAGTGGGGTc                                                                                                                                           >  1:705758/1‑150 (MQ=255)
                gTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTACGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGt                                                                                                                                        >  1:738644/1‑150 (MQ=255)
                 tCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTg                                                                                                                                       <  2:544040/150‑1 (MQ=255)
                   ttACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACTCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGa                                                                                                                                     <  1:779152/150‑1 (MQ=255)
                   ttACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCACGAGTGCGGTCCGTGGa                                                                                                                                     <  1:624717/150‑1 (MQ=255)
                        gCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACaa                                                                                                                                <  1:1603143/150‑1 (MQ=255)
                           tataCGGTCAATTTCAGCAATTAATTACAATATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCa                                                                                                                             <  1:1865658/150‑1 (MQ=255)
                             taCGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCaga                                                                                                                           <  1:328788/150‑1 (MQ=255)
                              aCGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGAGGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTATTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCagag                                                                                                                          <  1:991064/150‑1 (MQ=255)
                              aCGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAAAGATGGTACCTTCATACAGACACCATTTTTTCAACTCCGTTTTCCCCGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCagag                                                                                                                          <  2:924081/150‑1 (MQ=255)
                                ggTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGAc                                                                                                                        >  2:497668/1‑150 (MQ=255)
                                ggTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGAc                                                                                                                        >  2:52434/1‑150 (MQ=255)
                                ggTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGAc                                                                                                                        >  2:1050173/1‑150 (MQ=255)
                                  tCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTACGTGGATACAACCAGAGACCg                                                                                                                      >  1:473764/1‑150 (MQ=255)
                                     aTTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGAca                                                                                                                   >  1:1375686/1‑149 (MQ=255)
                                      tttCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTg                                                                                                                  <  2:266871/150‑1 (MQ=255)
                                         cAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGaca                                                                                                               >  2:1228580/1‑150 (MQ=255)
                                          aGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGacac                                                                                                              >  2:1009546/1‑150 (MQ=255)
                                           gCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACg                                                                                                             <  1:793834/150‑1 (MQ=255)
                                           gCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTAAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACg                                                                                                             >  2:1788976/1‑150 (MQ=255)
                                             aattaattACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAg                                                                                                           <  2:806588/150‑1 (MQ=255)
                                               ttaattaCATTATCCAAGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTc                                                                                                         <  2:1823827/150‑1 (MQ=255)
                                                    aaCATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTATCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAgg                                                                                                    <  1:397334/149‑1 (MQ=255)
                                                      cATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATTCAACCAGAGACCGACTGACACGAGTCAGAGGaa                                                                                                  >  2:1657287/1‑150 (MQ=255)
                                                        ttATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCGTATGTGAATAGTGCGGTCCGGGGGTACAACCCGAGACCGACGGACACGAGGCAGAGGAAac                                                                                                >  1:1536304/1‑150 (MQ=255)
                                                         tATCCACGCCAAAGGATATGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAacg                                                                                               <  2:966622/150‑1 (MQ=255)
                                                               cGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGat                                                                                         >  1:892333/1‑150 (MQ=255)
                                                                 ccAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCACCTCCGTTTTCCCCGGACCGCACGCTTTTGTCAAGTGTGCGGTCCGTGGATACACCCAGAGACCGGCTGACACGTGTGAGAGGAAACGCCGGatat                                                                                       >  2:633277/1‑150 (MQ=255)
                                                                 ccAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGatat                                                                                       <  2:9494/150‑1 (MQ=255)
                                                                    aaGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACGCTTATGTCAATAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGtt                                                                                    >  1:1597687/1‑150 (MQ=255)
                                                                     aGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGTATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTc                                                                                   <  2:487332/150‑1 (MQ=255)
                                                                        aTTTGTCATCACAATGATGGTACCTTCTTTGAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCTGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGt                                                                                <  2:1107763/150‑1 (MQ=255)
                                                                          tagTCCTCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCGGTTTTCCACGGACCGCACTCTTATGTCAAGACGGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCATAGGAAACGACGGATATGTTCAGTCg                                                                              <  2:310736/148‑1 (MQ=255)
                                                                                   caATGATGGTACCTTCTTTCAGCCACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAAtatc                                                                     >  2:963965/1‑150 (MQ=255)
                                                                                    aaTGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTAAGTCGTAAAAtatct                                                                    >  1:399937/1‑150 (MQ=255)
                                                                                     atgatgGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACAGCACTCTTATGTCAAGAGTTCGGTCCGTGTATACAACCATAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAAtatcta                                                                   <  1:80641/150‑1 (MQ=255)
                                                                                        atgGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCa                                                                <  1:1668914/150‑1 (MQ=255)
                                                                                         tgGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCaa                                                               <  2:98021/150‑1 (MQ=255)
                                                                                           gTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAATCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGTTATGTTCAGTCGTAAAATATCTATCaaaa                                                             >  2:424736/1‑150 (MQ=255)
                                                                                           gTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCaaaa                                                             <  2:1491097/150‑1 (MQ=255)
                                                                                                        tacacCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAgg                                                <  1:996583/149‑1 (MQ=255)
                                                                                                        gACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAgg                                                <  2:547759/150‑1 (MQ=255)
                                                                                                            ccATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCaa                                            >  1:588119/1‑150 (MQ=255)
                                                                                                             cATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTGAGTCGTAAAATTTAGATCAAAAAACATGATTAAGGTCaaa                                           >  2:1669223/1‑150 (MQ=255)
                                                                                                                     cAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTg                                   <  2:191925/150‑1 (MQ=255)
                                                                                                                     cAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTg                                   <  1:497668/150‑1 (MQ=255)
                                                                                                                     cAACTCCGTTTTCCACGGACCGCACGCTTATGTCAAGAGTGCGGTCCGTGGATCCAACCAGAGCCCGACAGCCACGAGTCAGAGGAAACGACGGCTATGGTCAGTCGTAAAATATCTATCAAAACACCTGCTTAAGGGCAAACATGTTTg                                   >  1:900611/1‑150 (MQ=255)
                                                                                                                     cAACTCCGTTTTCCACGGACAGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTg                                   <  1:52434/150‑1 (MQ=255)
                                                                                                                      aaCTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGa                                  >  2:105840/1‑150 (MQ=255)
                                                                                                                          ccGTTTTCCACGGACCGCACTCTTATGTATAGAGTGCGGTCCTTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATAtt                              <  2:1179591/150‑1 (MQ=255)
                                                                                                                           cGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATAttt                             <  2:1635297/150‑1 (MQ=255)
                                                                                                                            gTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTa                            <  2:1292517/150‑1 (MQ=255)
                                                                                                                             ttttCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTAc                           <  1:1050173/150‑1 (MQ=255)
                                                                                                                             ttttCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTAc                           <  2:1414802/150‑1 (MQ=255)
                                                                                                                              tttCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACa                          <  2:986444/150‑1 (MQ=255)
                                                                                                                                tCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAAt                        >  1:1462321/1‑150 (MQ=255)
                                                                                                                                 ccACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAAtt                       >  1:399639/1‑150 (MQ=255)
                                                                                                                                 ccACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAAtt                       >  1:533009/1‑150 (MQ=255)
                                                                                                                                   acagACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTa                     <  1:1322759/147‑1 (MQ=255)
                                                                                                                                    cGGACCGCACTCTTGACATAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTAt                    >  2:1773069/1‑150 (MQ=39)
                                                                                                                                    cGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTAt                    >  2:608934/1‑150 (MQ=255)
                                                                                                                                     ggACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATg                   >  2:1812096/1‑150 (MQ=255)
                                                                                                                                      gACCGCACTCTTGACATAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGa                  >  2:978685/1‑150 (MQ=255)
                                                                                                                                      gACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGa                  <  2:17920/150‑1 (MQ=255)
                                                                                                                                      gACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGa                  >  1:10377/1‑150 (MQ=255)
                                                                                                                                      gACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGa                  >  1:1427366/1‑150 (MQ=255)
                                                                                                                                      gACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGa                  >  2:1552415/1‑150 (MQ=255)
                                                                                                                                       aCCGCACTCTTGACATAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGAAATCTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGaa                 >  1:481902/1‑150 (MQ=255)
                                                                                                                                        ccGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAATAATGTTTGATATTTACAATTTATGAAg                >  1:847296/1‑150 (MQ=255)
                                                                                                                                          gCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGAt              >  1:1293608/1‑150 (MQ=255)
                                                                                                                                          gCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGAt              <  2:1536304/150‑1 (MQ=255)
                                                                                                                                           cACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATg             <  2:27004/150‑1 (MQ=255)
                                                                                                                                             ctctTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGAc           <  1:583113/150‑1 (MQ=255)
                                                                                                                                             ctctTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGAc           <  1:1086764/150‑1 (MQ=255)
                                                                                                                                                ttATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGAGTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACaat        >  2:115822/1‑150 (MQ=255)
                                                                                                                                                   tGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACaataat     >  1:634946/1‑150 (MQ=255)
                                                                                                                                                   tGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACaataat     >  2:185877/1‑150 (MQ=255)
                                                                                                                                                    gTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAAtt    <  1:500928/150‑1 (MQ=255)
                                                                                                                                                    gTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAAtt    >  1:1189862/1‑150 (MQ=255)
                                                                                                                                                     tCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAATta   <  1:1669223/150‑1 (MQ=255)
                                                                                                                                                     tCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAATta   >  1:1062941/1‑150 (MQ=255)
                                                                                                                                                      cAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAATtat  >  1:143566/1‑150 (MQ=255)
                                                                                                                                                      cAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGACAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAATtat  >  2:678157/1‑150 (MQ=255)
                                                                                                                                                      cAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGAGATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAATtat  >  2:672590/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CCTGATACTGATATGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAATTAT  >  NC_000913/2068261‑2068560

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