Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,068,407 A→G 5.0% pseudogene (13/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,4070AG 5.0% 349.3 / 5.6 122pseudogene (13/159 nt)yoeHpseudogene, CP4‑44 putative prophage remnant;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base A (52/61);  new base G (3/3);  total (58/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

Only 100 of 122 total aligned reads displayed.

GACCTGATACTGATATGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAAT  >  NC_000913/2068259‑2068554
                                                                                                                                                    |                                                                                                                                                   
gACCTGATACTGATATGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGAACTCTTAt                                                                                                                                                    >  1:27004/1‑150 (MQ=255)
 aCCTGATACTGATATGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATg                                                                                                                                                   <  1:370137/150‑1 (MQ=255)
 aCCTGATACTGATATGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATg                                                                                                                                                   <  1:243199/150‑1 (MQ=255)
  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)
           tatatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGc                                                                                                                                         <  2:972299/149‑1 (MQ=255)
           gataTGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGCCCGCACTCTTATGTCAAGAGTGc                                                                                                                                         >  2:408650/1‑150 (MQ=255)
           gataTGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGc                                                                                                                                         >  2:500928/1‑150 (MQ=255)
            atatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCg                                                                                                                                        <  2:769230/150‑1 (MQ=255)
             tatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCgg                                                                                                                                       <  2:783846/150‑1 (MQ=255)
             tatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCCCACTCTTATGTCAAGAGTGCgg                                                                                                                                       >  1:1179591/1‑150 (MQ=255)
             tatGAGTCTTACCGCATATACGGTCAATTTCAGCAAATAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTGACATAAGAGTGCgg                                                                                                                                       <  2:1870412/150‑1 (MQ=255)
              atGAGTCTTACCGCATATACTGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTGACATAAGAGTGCGGt                                                                                                                                      <  1:192430/150‑1 (MQ=255)
               tGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTc                                                                                                                                     >  1:220586/1‑150 (MQ=255)
               tGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTc                                                                                                                                     >  1:1292517/1‑150 (MQ=255)
               tGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTc                                                                                                                                     >  1:323497/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)
                      tACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTTCGGTCCGTGGAt                                                                                                                              >  2:1382589/1‑150 (MQ=255)
                       aCCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCCCGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATa                                                                                                                             <  1:1742701/150‑1 (MQ=255)
                          gCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGATTGCGGTCCGTGGATACac                                                                                                                          >  2:104046/1‑149 (MQ=255)
                          gCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACaa                                                                                                                          <  1:1603143/150‑1 (MQ=255)
                             tataCGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCa                                                                                                                       >  2:328788/1‑150 (MQ=255)
                             tataCGGTCAATTTCAGCAATTAATTACAATATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCa                                                                                                                       <  1:1865658/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)
                                    tCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTACGTGGATACAACCAGAGACCg                                                                                                                >  1:473764/1‑150 (MQ=255)
                                       aTTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGAca                                                                                                             >  1:1375686/1‑149 (MQ=255)
                                       aTTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACt                                                                                                             >  2:841920/1‑150 (MQ=255)
                                        tttCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTg                                                                                                            <  2:266871/150‑1 (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)
                                                                 cGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGat                                                                                   >  1:892333/1‑150 (MQ=255)
                                                                   ccAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCACCTCCGTTTTCCCCGGACCGCACGCTTTTGTCAAGTGTGCGGTCCGTGGATACACCCAGAGACCGGCTGACACGTGTGAGAGGAAACGCCGGatat                                                                                 >  2:633277/1‑150 (MQ=255)
                                                                   ccAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGatat                                                                                 <  2:9494/150‑1 (MQ=255)
                                                                    cAATGTATTTTTCATCAATATTATGTTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACAGCACTCTTATGTCAAGAGTGCGTTCCGTTGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATg                                                                                <  1:1615695/150‑1 (MQ=255)
                                                                      aaGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACGCTTATGTCAATAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGtt                                                                              >  1:1597687/1‑150 (MQ=255)
                                                                         tATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGGCAGAGGAAACGACGGATATGTTCAg                                                                           <  1:392059/150‑1 (MQ=255)
                                                                          aTTTGTCATCACAATGATGGTACCTTCTTTGAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCTGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGt                                                                          <  2:1107763/150‑1 (MQ=255)
                                                                          aTTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTGTTTTTTCATGATTGCGGTCCGTGGTTTCAACAATAGTGCGACTTACACGATTCAGAGGAAAGGACGGAGTTGTTCGGt                                                                          >  1:1166100/1‑150 (MQ=255)
                                                                            tagTCCTCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCGGTTTTCCACGGACCGCACTCTTATGTCAAGACGGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCATAGGAAACGACGGATATGTTCAGTCg                                                                        <  2:310736/148‑1 (MQ=255)
                                                                                  tcaCAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACAAGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAat                                                                  >  1:817955/1‑150 (MQ=255)
                                                                                     caATGATGGTACCTTCTTTCAGCCACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAAtatc                                                               >  2:963965/1‑150 (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)
                                                                                              tACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCaaaaa                                                      <  1:1475035/150‑1 (MQ=255)
                                                                                                 cttcttTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACTGATATGTTCAGTCGTAAAATATCTATCAAAAAACa                                                   <  2:1622624/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)
                                                                                                               cATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGGAAAATATCTATCAAAAAACATGATTAAGGTCaaa                                     >  1:734289/1‑150 (MQ=255)
                                                                                                                      tCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGttt                              <  1:521601/150‑1 (MQ=255)
                                                                                                                       cAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTg                             <  1:497668/150‑1 (MQ=255)
                                                                                                                       cAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTg                             <  2:191925/150‑1 (MQ=255)
                                                                                                                       cAACTCCGTTTTCCACGGACCGCACGCTTATGTCAAGAGTGCGGTCCGTGGATCCAACCAGAGCCCGACAGCCACGAGTCAGAGGAAACGACGGCTATGGTCAGTCGTAAAATATCTATCAAAACACCTGCTTAAGGGCAAACATGTTTg                             >  1:900611/1‑150 (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)
                                                                                                                                tttCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACa                    <  2:1159683/150‑1 (MQ=255)
                                                                                                                                  tCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAAt                  >  1:1462321/1‑150 (MQ=255)
                                                                                                                                   ccACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAAtt                 >  1:533009/1‑150 (MQ=255)
                                                                                                                                   ccACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAAtt                 >  1:399639/1‑150 (MQ=255)
                                                                                                                                     acagACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTa               <  1:1322759/147‑1 (MQ=255)
                                                                                                                                      cGGACCGCACTCTTGACATAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTAt              >  2:1773069/1‑150 (MQ=39)
                                                                                                                                      cGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTAt              <  1:618675/150‑1 (MQ=255)
                                                                                                                                      cGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTAt              >  2:608934/1‑150 (MQ=255)
                                                                                                                                       ggACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATg             >  2:1812096/1‑150 (MQ=255)
                                                                                                                                        gACCGCACTCTTGACATAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGa            >  2:978685/1‑150 (MQ=255)
                                                                                                                                        gACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGa            >  2:1552415/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)
                                                                                                                                         aCCGCACTCTTGACATAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGAAATCTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGaa           >  1:481902/1‑150 (MQ=255)
                                                                                                                                          ccGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAATAATGTTTGATATTTACAATTTATGAAg          >  1:847296/1‑150 (MQ=255)
                                                                                                                                            gCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGAt        >  1:1293608/1‑150 (MQ=255)
                                                                                                                                             cACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATg       >  1:477358/1‑150 (MQ=255)
                                                                                                                                               ctctTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGAc     <  1:1086764/150‑1 (MQ=255)
                                                                                                                                               ctctTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGAc     <  1:583113/150‑1 (MQ=255)
                                                                                                                                                  ttATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGAGTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACaat  >  2:115822/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GACCTGATACTGATATGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAAT  >  NC_000913/2068259‑2068554

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: