Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,068,411 C→T 9.4% 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 9.4% 301.6 / 17.9 121pseudogene (9/159 nt)yoeHpseudogene, CP4‑44 putative prophage remnant;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (61/43);  new base T (4/7);  total (71/50)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-03

Only 100 of 123 total aligned reads displayed.

GATACTGATATGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAATTAT  >  NC_000913/2068264‑2068560
                                                                                                                                                   |                                                                                                                                                     
gatactgataTGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCaa                                                                                                                                                     >  2:1435112/1‑150 (MQ=255)
      gataTGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGTTACCTTCTATCAGACACCATTTTTTCAACTCAGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGc                                                                                                                                               >  1:97320/1‑150 (MQ=255)
      gataTGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACTGACCGCACTCTTATGTCAAGAGTGc                                                                                                                                               >  1:1723394/1‑150 (MQ=255)
       atatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCg                                                                                                                                              >  2:327070/1‑150 (MQ=255)
        tatGAGTCTTACCTCATATACGGTCAATTTCAGAAATTACTTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTGACATAAGAGTGCgg                                                                                                                                             <  2:120401/150‑1 (MQ=255)
        tatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTGACATAAGAGTGCgg                                                                                                                                             <  1:802229/150‑1 (MQ=255)
        tatGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTGACATAAGAGTGCgg                                                                                                                                             <  1:561024/150‑1 (MQ=255)
         atGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTGACATAAGAGTGCGGt                                                                                                                                            <  1:1355367/150‑1 (MQ=255)
          tGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATAATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTc                                                                                                                                           >  2:703487/1‑150 (MQ=255)
           gAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCATTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTcc                                                                                                                                          >  2:1194827/1‑150 (MQ=255)
              tCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTg                                                                                                                                       <  1:993100/150‑1 (MQ=255)
              tCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTAGACATAAGAGTGCGGTCCGTg                                                                                                                                       <  1:1489797/150‑1 (MQ=39)
               cTTACCGCATATACGGTCAATTTCAGCAATTAATAACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCAATTTTACAAATCCGTTTTCCACGGACCGCACTCTAATGTCAAGAGTGCGGTCCGTgg                                                                                                                                      >  1:1132158/1‑150 (MQ=255)
                  aCCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGCCCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATa                                                                                                                                   >  2:942227/1‑150 (MQ=255)
                    cGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACa                                                                                                                                 <  2:561640/150‑1 (MQ=255)
                     gCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCTAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACaa                                                                                                                                >  1:288813/1‑150 (MQ=255)
                        tataCGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCa                                                                                                                             <  2:1237148/150‑1 (MQ=255)
                        tataCGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGAAACCATTTTTTCAACTCAGTTTTCCACGGACCGCAATATTATTTAAAGAGTTCCGTACGTGGAaacaacaa                                                                                                                             >  1:603022/1‑148 (MQ=255)
                         ataCGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCag                                                                                                                            >  1:227798/1‑150 (MQ=255)
                         ataCGGTCAATTTCAGCAATTAATTACATTAACCACGCCAAAGTATTTGTCATCACAATGATGTTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCag                                                                                                                            <  1:1524220/150‑1 (MQ=255)
                           aCGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTAAAGAGTGCGGTCCGTGGATACAACCagag                                                                                                                          >  2:596657/1‑150 (MQ=255)
                             ggTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGAc                                                                                                                        <  2:457514/150‑1 (MQ=255)
                             ggTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACAGCACTCTTATGTCAAGAGTCCGGTCCCTGGAAACAACCAGAAAc                                                                                                                        >  1:478125/1‑150 (MQ=255)
                              gTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTAATCACAATGATGGTACCTTCTTTAATACACCATTTTTTAAACTAAGTATTCAACGGACAGCACTATTAAGACAAGAGTGATCTCCGTGTATACAACCagataac                                                                                                                       >  1:1216507/1‑146 (MQ=255)
                                 aaTTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTACGTGGATACAACCAGAGACCGaa                                                                                                                    >  1:1372579/1‑149 (MQ=255)
                                   tttCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTg                                                                                                                  <  1:913779/150‑1 (MQ=255)
                                   tttCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTg                                                                                                                  <  2:756802/150‑1 (MQ=255)
                                      cAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGaca                                                                                                               <  1:942227/150‑1 (MQ=255)
                                       aGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACCGACCGCACTCTTATGTAAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGacac                                                                                                              >  1:367275/1‑150 (MQ=255)
                                        gCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGacacc                                                                                                             >  1:624283/1‑149 (MQ=255)
                                        gCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACg                                                                                                             >  1:55259/1‑150 (MQ=255)
                                        gCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACACTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACg                                                                                                             >  1:1002756/1‑150 (MQ=255)
                                         cAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGa                                                                                                            >  1:1441753/1‑150 (MQ=255)
                                          aattaattACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACACTTATGTCAAGAGTGCGGTCCGTGGATTCAACCAGAGGCCCACTCACACGaa                                                                                                           >  1:592031/1‑149 (MQ=255)
                                            ttaattaCATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTc                                                                                                         <  1:1722493/150‑1 (MQ=255)
                                               attaCATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCaga                                                                                                      <  2:1354622/150‑1 (MQ=255)
                                                           aCGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTAAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGa                                                                                          >  1:1113478/1‑150 (MQ=255)
                                                            cGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTACGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAAAGACGGat                                                                                         >  2:195120/1‑150 (MQ=255)
                                                             gCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGata                                                                                        >  1:1189087/1‑150 (MQ=255)
                                                               cAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATg                                                                                      >  1:358104/1‑150 (MQ=255)
                                                                aaaGTATTTGTCATCACAATGATGGTACCTCCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGt                                                                                     >  2:768198/1‑150 (MQ=255)
                                                                  aGTATTTGTCATCACAATGATGGTTCCTTTTTTCAGACACCATTTTTTCAACTCCGATTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTc                                                                                   <  2:1688265/150‑1 (MQ=255)
                                                                    tATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCGTATGTCAAGAGTGCGGTCCGTGGCTAAAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAg                                                                                 >  1:1362067/1‑150 (MQ=255)
                                                                     aTTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGt                                                                                <  2:1442710/150‑1 (MQ=255)
                                                                     aTTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGt                                                                                >  1:1392990/1‑150 (MQ=255)
                                                                        tGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCCCGCTTATGTCAAGCGTGCGGTCCGTGGATTCAACCCGAGACCGCCTGCCACGAGTCAGAGGAAACGACGGATATGTTCAGTCgg                                                                             >  1:1398604/1‑149 (MQ=255)
                                                                        tGTCATCACAATGATGGTACCTTATTTCAGACACCATTTTTTCAACTCCGCTTTCCACGGACCGCACTCTTATGTCAAGAGGGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGt                                                                             <  2:760446/150‑1 (MQ=255)
                                                                         gTCAGCACACTGATGGTACCTTCTTTCAGACACCATTTTTTCACCTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTa                                                                            <  1:116116/150‑1 (MQ=255)
                                                                           catcaCAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTaaa                                                                          >  2:615616/1‑150 (MQ=255)
                                                                           catcaCAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTaaa                                                                          >  2:895330/1‑150 (MQ=255)
                                                                            atcaCAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTaaaa                                                                         <  2:227798/150‑1 (MQ=255)
                                                                             tcaCAATGATGGTACATTCTTTCAGTCACGATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGTTTGCGGTCCGTGGATACAACCATAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAat                                                                        <  1:759375/150‑1 (MQ=255)
                                                                                 aaTGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACTGACCGCACTCTTATGTCAAGAGTGCGGTCCGAGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAAtatct                                                                    <  1:327070/150‑1 (MQ=255)
                                                                                 aaTGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAAtatct                                                                    <  2:785131/150‑1 (MQ=255)
                                                                                  atgatgGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAAtatcta                                                                   >  2:130867/1‑150 (MQ=255)
                                                                                      tgGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTTTGTCAAGAGTGCGGGCCGTGGATACCACCAGAGACCGACAGACACGAGTCAGGGGAAACGACGGATAAGTTCAGTCGGAAAATATATATCaa                                                               >  2:1095798/1‑150 (MQ=255)
                                                                                               cttTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGa                                                      >  1:312119/1‑150 (MQ=255)
                                                                                                  tCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGCCACGAGTCAGAGGAAACGCCGGATATGTTCAATAGTAAAATATCAATCAAAAAACATGATTa                                                   >  1:861899/1‑150 (MQ=255)
                                                                                                    aGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAg                                                 >  2:631443/1‑150 (MQ=255)
                                                                                                     gACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAgg                                                <  2:1098258/150‑1 (MQ=255)
                                                                                                     gACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAgg                                                >  1:1048035/1‑150 (MQ=255)
                                                                                                     gACACCATTTTTTCAACTCCGTTATCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAgg                                                <  1:1435112/150‑1 (MQ=255)
                                                                                                      acacCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGt                                               <  2:546053/150‑1 (MQ=255)
                                                                                                       cagcATTTTTTCAACTCCGTTTTCCACTGACCGCACTCTTATGTCAAGAGTGCGGTCTGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTc                                              <  1:962643/147‑1 (MQ=255)
                                                                                                       cacCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTc                                              >  2:585566/1‑150 (MQ=255)
                                                                                                       cacCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTAAGTCGTAAAATATCTATCAAAAAACATGATTAATGTc                                              >  2:325425/1‑150 (MQ=255)
                                                                                                        acCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGCCCGACTGCCACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCa                                             >  1:1382921/1‑150 (MQ=255)
                                                                                                         ccATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCaa                                            >  1:955957/1‑150 (MQ=255)
                                                                                                             tttttCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAAt                                        <  1:540479/150‑1 (MQ=255)
                                                                                                                 tCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACAAGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGttt                                    <  1:1140938/150‑1 (MQ=255)
                                                                                                                      tCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGatat                               >  1:234992/1‑150 (MQ=255)
                                                                                                                         tttttCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGAAAGACACGAGTCAGAGGAAACGACGGATATGTTCAGTATTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTa                            <  1:1124781/149‑1 (MQ=255)
                                                                                                                          ttttCCACGGACCGCGCTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTAc                           >  1:898742/1‑150 (MQ=255)
                                                                                                                          ttttCCACGGACCGCACTCTTATGTCAATAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCATAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTAc                           <  2:1534298/150‑1 (MQ=255)
                                                                                                                          ttttCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTAc                           <  2:1716620/150‑1 (MQ=255)
                                                                                                                              ccACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAAtt                       <  2:601988/150‑1 (MQ=255)
                                                                                                                              ccACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGATTGATATTTAAAAtt                       >  2:277800/1‑150 (MQ=255)
                                                                                                                               cACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAAttt                      >  2:932839/1‑150 (MQ=255)
                                                                                                                                 cGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTAt                    <  2:358104/150‑1 (MQ=255)
                                                                                                                                 cGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTAt                    <  1:1095798/150‑1 (MQ=255)
                                                                                                                                  ggACCGCACTCTTGACATAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGCCACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACATTTTATg                   >  2:1023580/1‑150 (MQ=39)
                                                                                                                                   gACCGCACTCTTGACATAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGa                  >  2:1542993/1‑150 (MQ=255)
                                                                                                                                    aCCGCACTCTTGACATAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGaa                 >  1:1468504/1‑150 (MQ=255)
                                                                                                                                    aCCGCACTCTTGACATAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGaa                 >  1:1034194/1‑150 (MQ=255)
                                                                                                                                     ccGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAg                <  1:1104120/150‑1 (MQ=255)
                                                                                                                                       gCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGac              >  2:488417/1‑149 (MQ=255)
                                                                                                                                       gCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGAt              <  1:1210464/150‑1 (MQ=255)
                                                                                                                                       gCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGAt              >  1:1030207/1‑150 (MQ=255)
                                                                                                                                       gCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGAt              >  1:1274355/1‑150 (MQ=255)
                                                                                                                                       gCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGAt              >  2:1553235/1‑150 (MQ=255)
                                                                                                                                       gCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGATTGATATTTACAATTTATGAAGAt              >  2:1516090/1‑150 (MQ=255)
                                                                                                                                       gCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGACAAAAATGTTTGATATTTAAAATTTATGAAGAt              >  1:1274364/1‑150 (MQ=255)
                                                                                                                                           tctTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACa          <  1:596657/150‑1 (MQ=255)
                                                                                                                                             ttATGTCAAGAGTGCGGTGCGTGGATACAACCAGAGACCGAATGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACaat        <  1:1643527/150‑1 (MQ=255)
                                                                                                                                                 gTCAAGAGTGCGGTCCGTGTATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAAtt    <  1:1623082/150‑1 (MQ=255)
                                                                                                                                                 gTCAAGAGTGCGGTCCGTGGATACGACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAATACATGATTAAGGTCAAAATTGTTTGATATTTACAATTTTTGAATATGACAattatt    >  1:1603342/1‑150 (MQ=255)
                                                                                                                                                 gTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGATACGACGGATATGTTCAGTCGTAAACTATCTATCACAAACCAGGAGTAAGGGCAAAAATGTTTGAGATTTACAATTTATGAAGATGACAATAAtt    >  1:1596910/1‑150 (MQ=255)
                                                                                                                                                 gTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATACCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAAtt    >  2:1463051/1‑150 (MQ=255)
                                                                                                                                                   cAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAATtat  >  2:401911/1‑150 (MQ=255)
                                                                                                                                                   cAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAATtat  <  2:1441753/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GATACTGATATGAGTCTTACCGCATATACGGTCAATTTCAGCAATTAATTACATTATCCACGCCAAAGTATTTGTCATCACAATGATGGTACCTTCTTTCAGACACCATTTTTTCAACTCCGTTTTCCACGGACCGCACTCTTATGTCAAGAGTGCGGTCCGTGGATACAACCAGAGACCGACTGACACGAGTCAGAGGAAACGACGGATATGTTCAGTCGTAAAATATCTATCAAAAAACATGATTAAGGTCAAAAATGTTTGATATTTACAATTTATGAAGATGACAATAATTAT  >  NC_000913/2068264‑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: