Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,212,088 (C)8→9 intergenic (‑85/+615) iraM ← / ← ycgX RpoS stabilzer during Mg starvation, anti‑RssB factor/DUF1398 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,212,0801.C82.1% 157.4 / 30.0 56intergenic (‑77/+623)iraM/ycgXRpoS stabilzer during Mg starvation, anti‑RssB factor/DUF1398 family protein
Reads supporting (aligned to +/- strand):  ref base . (6/4);  new base C (21/25);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TCTTTATCCTTAAGGATAATACCAGACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATACATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATA  >  NC_000913/1211837‑1212322
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
tCTTTATCCTTAAGGATAATACCAGACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTa                                                                                                                                                                                                                                                          <  2:488167/244‑1 (MQ=255)
        ataagggacagtaCCAGACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTTAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCaaaa                                                                                                                                                                                                                                      <  2:898987/240‑1 (MQ=255)
          tAAGGATAATACCAGACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTTTTGATTTTGCATATGAGTATATTA‑‑ccc                                                                                                                                                                                                                                                   >  2:1172044/1‑237 (MQ=255)
                        gACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAaatt                                                                                                                                                                                                                             >  1:488207/1‑244 (MQ=255)
                          cTTAAACGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTTCCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGGATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAat                                                                                                                                                                                                                    <  1:747141/251‑1 (MQ=255)
                              aCCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑ccc                                                                                                                                                                                                                                                   >  2:1195497/1‑217 (MQ=255)
                                     tAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTtata                                                                                                                                                                                                         >  1:1927343/1‑251 (MQ=255)
                                          ataCTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑c                                                                                                                                                                                                                                                       >  2:302755/1‑203 (MQ=255)
                                           taCTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATTATGGc                                                                                                                                                                                                              >  1:899055/1‑240 (MQ=255)
                                           taCTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CTCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGATTATataaaat                                                                                                                                                                                                   >  1:1888541/1‑251 (MQ=255)
                                               gCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATAtaaaataaaa                                                                                                                                                                                               <  1:2082002/251‑1 (MQ=255)
                                                ccAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTTTATaaaataaaat                                                                                                                                                                                              >  2:747082/1‑251 (MQ=255)
                                                   ggagggagAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATataaaat                                                                                                                                                                                                   <  1:1557280/243‑1 (MQ=255)
                                                               ataGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTAGTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAAtc                                                                                                                                                                               <  2:1888345/251‑1 (MQ=255)
                                                                 agtagATTGATACCAGATTATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAAtctc                                                                                                                                                                             <  2:1488908/251‑1 (MQ=255)
                                                                     gATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACTGTTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGTGTATATTT‑‑CCCCCCCCttcaaa                                                                                                                                                                                                                                      <  1:834582/189‑6 (MQ=255)
                                                                           tACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAATCAAGGAATCTaaa                                                                                                                                                                           >  2:1557135/1‑243 (MQ=255)
                                                                                       aTTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATgag                                                                                                                                                       >  1:1489042/1‑251 (MQ=255)
                                                                                             atatTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCACCCCTTCAAAAAAATAACTTCATTAAAATGATgg                                                                                                                                                                                                               >  2:834515/1‑189 (MQ=255)
                                                                                             atatTACCCAATATGGCTGAAAAAGATATACAACATGTAGGTTGAAGTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTAACCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGTTTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATga                                                                                                                                                        <  2:1933695/244‑1 (MQ=255)
                                                                                                                     gATATACCACATGTAGGTTGAATTACCGTGTCCATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTCTTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTCATTAAAATTATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAtttat                                                                                                                                >  1:1933894/1‑244 (MQ=255)
                                                                                                                           ccACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGAAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTa                                                                                                                   >  2:2189616/1‑251 (MQ=255)
                                                                                                                             acaTGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑cccccccc                                                                                                                                                                                                                                            >  1:55567/1‑127 (MQ=255)
                                                                                                                                tttAGGTTGAATTACCGTGTCAATTACTATCCCATTCATTTGTTATTGCTTTTCCCACGGTTTTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATTGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGAttt                                                                                                              <  2:812302/249‑1 (MQ=255)
                                                                                                                                     gTTGAATTACCGTGTCAATTACTATCCCCTTAATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTa                                                                                                         <  2:1927144/251‑1 (MQ=255)
                                                                                                                                      ttGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAAtt                                                                                                                                                                                             <  1:923975/166‑1 (MQ=255)
                                                                                                                                           ttACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCGCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAACTATTGACTTCGTTTTAGATTTGTTTAGCtata                                                                                                   >  2:1671767/1‑251 (MQ=255)
                                                                                                                                             aCCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAAt                                                                                                 >  1:812365/1‑251 (MQ=255)
                                                                                                                                                  gtCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAg                                                                                                                  <  1:1647344/229‑1 (MQ=255)
                                                                                                                                                  gtCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCGCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCg                                                                                                                        <  1:1929638/223‑1 (MQ=255)
                                                                                                                                                   tCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTtata                                                                                           >  1:878330/1‑251 (MQ=255)
                                                                                                                                                   tCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTtata                                                                                           <  1:1761369/251‑1 (MQ=255)
                                                                                                                                                      aTTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAAATAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGAAGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGAATTCGTTTTAGATTTGTTTAGCTATAATGTTATACAtt                                                                                       >  1:1658274/1‑251 (MQ=255)
                                                                                                                                                          cTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATACATTCaa                                                                                    >  1:49283/1‑251 (MQ=255)
                                                                                                                                                            aTCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATACATTCAAAt                                                                                  <  1:2189836/251‑1 (MQ=255)
                                                                                                                                                              ccACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATg                                                                                                                                                                     >  2:923901/1‑166 (MQ=255)
                                                                                                                                                                          ttATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATAc                                                                                          >  2:1647180/1‑229 (MQ=255)
                                                                                                                                                                          ttATGTCTTATCCCACGGTATTTAATATGGTTCAGTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTCTTGACTTCGTTTTAGATTTGTTTAGCTATAATg                                                                                                >  2:1929439/1‑223 (MQ=255)
                                                                                                                                                                              gTCTTATCCCACGGTATTTAATATGGTTCATTAGGATTTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGAtt                                                                                                                                                    <  2:1053840/167‑1 (MQ=255)
                                                                                                                                                                              gTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATAtaaaataaaa                                                                                                                                                                                               <  1:1997633/124‑1 (MQ=255)
                                                                                                                                                                                  tATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATACATTCAAATGACTGAACATCCTGTAATTaaa                                                            >  2:1276814/1‑251 (MQ=255)
                                                                                                                                                                                         cGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑ccccccccc                                                                                                                                                                                                                                           >  1:2160845/1‑67 (MQ=255)
                                                                                                                                                                                                 aaTATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAATGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATACATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGc                                             <  1:1671933/251‑1 (MQ=255)
                                                                                                                                                                                                      ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATgg                                                                                                                                                                                                               <  2:739830/84‑1 (MQ=255)
                                                                                                                                                                                                      ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGa                                                                                                                                                                       >  2:1997430/1‑124 (MQ=255)
                                                                                                                                                                                                      ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGAc                                                                                                                            >  1:1053922/1‑167 (MQ=255)
                                                                                                                                                                                                      ggTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCACCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATACATTCAAATGACTGAACATCCTGTAATTAAAACATAGCATTTATGCTACtt                                        >  2:1544254/1‑251 (MQ=255)
                                                                                                                                                                                                       gTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAaattaa                                                                                                                                                                                                                           <  2:1961734/72‑1 (MQ=255)
                                                                                                                                                                                                       gTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATACATTCAAATGACTGa                                                                            <  2:213380/214‑1 (MQ=255)
                                                                                                                                                                                                       gTTCATTAGAATGTTTATTTCTTGATTTTGCATATGAGTATCTTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAATAAAATGAGATTTTGTGAAAGCAATAAATGATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATACATTCAAATGACTGAACATCCTGTAATTAAACCATAGCCTTTATGCTACttt                                       <  2:1464819/251‑1 (MQ=255)
                                                                                                                                                                                                          cATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCg                                                                                                                        <  2:1042548/166‑1 (MQ=255)
                                                                                                                                                                                                                   gTTTATCACTTGATTTTGCATATGAGTATATTA‑‑CACCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGTTTATAGAAAATAAAATTTAAAACAAGGAATCCCAATGGATGTTAAACAAAATGAGAGTTTGTGGAAGCAATAAATTATTGACTTCGGTTTAGATTGGTTTAGATCTAATGTTATACATTCAAATGAATGGACATTCTGGAATTAAAACATAGCTTTTATGCTACTTTGTGCCAATTTGc                           >  1:1464953/1‑251 (MQ=255)
                                                                                                                                                                                                                              gATTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGc                                                                                                                                                                                       >  1:739889/1‑84 (MQ=255)
                                                                                                                                                                                                                               aTTTTGCATATGAGTATATTACCCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATatcaaat                                                                                                                                                                                                   >  1:1961933/1‑72 (MQ=255)
                                                                                                                                                                                                                               aTTTTGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATACATTCAAATGACTGAACATCCTGTAATTAAAACATAGcc                                                    >  1:213400/1‑214 (MQ=255)
                                                                                                                                                                                                                                ttttGCATATGAGTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATACATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGt              <  1:1276905/251‑1 (MQ=255)
                                                                                                                                                                                                                                  ttGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATg                                                                                                >  1:1042631/1‑166 (MQ=255)
                                                                                                                                                                                                                                            gTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTa                                                                                                                                                                  <  1:125153/91‑1 (MQ=255)
                                                                                                                                                                                                                                            gTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTg                                                                                                                              <  1:1171717/127‑1 (MQ=255)
                                                                                                                                                                                                                                            gTATATTACCCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTtat                                                                                            <  1:353733/161‑1 (MQ=255)
                                                                                                                                                                                                                                            gTATATTA‑‑CCACCCCCCTCAAAAAAATAAATTAATTAAAATGATGACACATATAAAATAAAATTTAAAGCAAAGAATCATAATGGATTTTAAACAAAATGAGATTTTGTTAAAGCAATAAATTATTGACTACGTTTTAGATTTGTTTAGCTATAATGTTATACATACAAATGACTGACCATCCTGTAAATAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATGATGGTTGCCTTTTtata  <  2:1575554/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gTATATTA‑‑CCCCCCACCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATACATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTtata  <  2:684518/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                       
TCTTTATCCTTAAGGATAATACCAGACTTAACCGGTGTAAATATACTGCCAGGAGGGAGAAATATAGTAGATTGATACCAGATGATCATTTTCATATTACCCCATATGGCTGAAAAAGATATACCACATGTAGGTTGAATTACCGTGTCAATTACTATCCACTTCATTTGTTATGTCTTATCCCACGGTATTTAATATGGTTCATTAGGATGTTTATTTCTTGATTTTGCATATGAGTATATTA‑‑CCCCCCCCTCAAAAAAATAAATTAATTAAAATGATGGCTTATATAAAATAAAATTTAAAGCAAGGAATCTCAATGGATGTTAAACAAAATGAGATTTTGTGAAAGCAATAAATTATTGACTTCGTTTTAGATTTGTTTAGCTATAATGTTATACATTCAAATGACTGAACATCCTGTAATTAAAACATAGCCTTTATGCTACTTTGTGCCAATTTGCTAAACATTATGGTTGCCTTTTTATA  >  NC_000913/1211837‑1212322

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