Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,346,764 A→C 7.0% intergenic (‑22/‑32) gmr ← / → ymiB cyclic‑di‑GMP phosphodiesterase; csgD regulator; modulator of RNase II stability/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,346,7640AC 7.0% 348.8 / 17.3 100intergenic (‑22/‑32)gmr/ymiBcyclic‑di‑GMP phosphodiesterase; csgD regulator; modulator of RNase II stability/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (40/53);  new base C (3/4);  total (43/57)
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.07e-01

Only 100 of 104 total aligned reads displayed.

ACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAAT  >  NC_000913/1346617‑1346905
                                                                                                                                                   |                                                                                                                                             
aCGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAg                                                                                                                                              <  4:1890673/149‑1 (MQ=255)
aCGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAg                                                                                                                                              <  2:99682/149‑1 (MQ=255)
  gATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCg                                                                                                                                            <  1:1126919/149‑1 (MQ=255)
    tCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGa                                                                                                                                          <  1:532260/149‑1 (MQ=255)
    tCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGa                                                                                                                                          <  4:376775/149‑1 (MQ=255)
    tCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGa                                                                                                                                          <  2:1418619/149‑1 (MQ=255)
    tCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGa                                                                                                                                          <  3:68265/149‑1 (MQ=255)
    tCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGa                                                                                                                                          >  1:1742987/1‑149 (MQ=255)
           cTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCcgcg                                                                                                                                   <  3:824812/149‑1 (MQ=255)
            tgtgGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCgg                                                                                                                                  >  2:1137052/1‑149 (MQ=255)
            tgtgGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCgg                                                                                                                                  <  1:168505/149‑1 (MQ=255)
             gtgGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGt                                                                                                                                 >  4:1045577/1‑149 (MQ=255)
             gtgGTTACGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGt                                                                                                                                 >  2:1334978/1‑149 (MQ=255)
              tgGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACCGCGGACTCCGCTGtt                                                                                                                                <  1:302839/149‑1 (MQ=255)
              tgGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGtt                                                                                                                                >  1:775540/1‑149 (MQ=255)
               ggTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTta                                                                                                                               >  2:560730/1‑149 (MQ=255)
                 tttCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaa                                                                                                                             <  2:945902/149‑1 (MQ=255)
                     ggtggtAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaataaa                                                                                                                         <  1:560730/149‑1 (MQ=255)
                      gtggtAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACCGCGGACTCCGCTGtt                                                                                                                                >  4:1254526/1‑141 (MQ=255)
                      gtggtAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACCGCGGACTCCGCTGtt                                                                                                                                <  3:1254526/141‑1 (MQ=255)
                      gtggtAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATaaaa                                                                                                                        <  2:832774/149‑1 (MQ=255)
                      gtggtAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATaaaa                                                                                                                        <  2:449096/149‑1 (MQ=255)
                      gtggtAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATaaaa                                                                                                                        <  1:1137052/149‑1 (MQ=255)
                         gtAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGTTTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCg                                                                                                                                            >  4:1545513/1‑126 (MQ=255)
                         gtAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCg                                                                                                                                            <  3:1545513/126‑1 (MQ=255)
                         gtAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCg                                                                                                                                            <  3:1836960/126‑1 (MQ=255)
                         gtAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCg                                                                                                                                            >  4:1836960/1‑126 (MQ=255)
                         gtAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGa                                                                                                                     <  3:711549/149‑1 (MQ=255)
                          tAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaa                                                                                                                    <  2:775540/149‑1 (MQ=255)
                            gAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaaca                                                                                                                  <  3:1045577/149‑1 (MQ=255)
                            gAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaaca                                                                                                                  <  1:1334978/149‑1 (MQ=255)
                              aaaGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaacaac                                                                                                                >  3:1539023/1‑149 (MQ=255)
                               aaGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAacaaca                                                                                                               >  3:284080/1‑149 (MQ=255)
                                 gcgcAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATaaaa                                                                                                                        <  3:632011/138‑1 (MQ=255)
                                 gcgcAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATaaaa                                                                                                                        >  4:632011/1‑138 (MQ=255)
                                 gcgcAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATg                                                                                                             >  3:1416459/1‑149 (MQ=255)
                                  cgcAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGa                                                                                                            <  1:1112395/149‑1 (MQ=255)
                                    cAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACt                                                                                                          <  4:1539023/149‑1 (MQ=255)
                                     aaaaCATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCg                                                                                                                                            >  4:357726/1‑114 (MQ=255)
                                     aaaaCATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCg                                                                                                                                            <  3:357726/114‑1 (MQ=255)
                                      aaaCATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTa                                                                                                        <  4:284080/149‑1 (MQ=255)
                                       aaCATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACCGCGGACTCCGCTGTtaac                                                                                                                             >  2:1417133/1‑126 (MQ=255)
                                       aaCATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACCGCGGACTCCGCTGTtaac                                                                                                                             >  1:1060475/1‑126 (MQ=255)
                                       aaCATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACCGCGGACTCCGCTGTtaac                                                                                                                             <  2:1060475/127‑2 (MQ=255)
                                       aaCATCGCTGCTTTCCGTCAACCGCCAGTATGGACTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACCGCGGACTCCGCTGTtaac                                                                                                                             <  1:1417133/127‑2 (MQ=255)
                                          aTCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaa                                                                                                                             >  1:1297441/1‑124 (MQ=255)
                                          aTCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaa                                                                                                                             <  2:1297441/124‑1 (MQ=255)
                                             gctgctTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGaaa                                                                                                 <  3:950404/149‑1 (MQ=255)
                                               tgctTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGaaa                                                                                                 >  3:425872/1‑147 (MQ=255)
                                               tgctTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGaaa                                                                                                 <  4:425872/147‑1 (MQ=255)
                                                  tttCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAAc                                                                                            <  4:1416459/149‑1 (MQ=255)
                                                   ttCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTAGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGTTTAAATAAAAGGAACAACATGACTTACATGAAATTAAcg                                                                                           >  2:1606574/1‑149 (MQ=255)
                                                       gTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGa                                                                                                                                          >  4:1779924/1‑98 (MQ=255)
                                                       gTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGa                                                                                                                                          <  3:1779924/98‑1 (MQ=255)
                                                        tCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACa                                                                                                                                               >  2:613101/1‑92 (MQ=255)
                                                        tCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACa                                                                                                                                               <  1:613101/92‑1 (MQ=255)
                                                              gCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATg                                                                                                             >  1:1543089/1‑120 (MQ=255)
                                                              gCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATg                                                                                                             <  2:1543089/120‑1 (MQ=255)
                                                                 aGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGTTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCAt                                                                             <  1:1606574/149‑1 (MQ=255)
                                                                   tATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACAt                                                                                                              <  1:688806/114‑1 (MQ=255)
                                                                   tATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTACGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACAt                                                                                                              >  2:688806/1‑114 (MQ=255)
                                                                    aTGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaacaac                                                                                                                >  3:366070/1‑111 (MQ=255)
                                                                    aTGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaacaac                                                                                                                <  4:366070/111‑1 (MQ=255)
                                                                      ggttggtgCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGtt                                                                         >  4:284974/6‑148 (MQ=255)
                                                                      ggATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACTGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGtt                                                                         <  3:284974/148‑1 (MQ=255)
                                                                       gATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATaaaa                                                                                                                        <  1:74925/100‑1 (MQ=255)
                                                                       gATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATaaaa                                                                                                                        >  2:74925/1‑100 (MQ=255)
                                                                        aTTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGttttt                                                                      >  3:71198/1‑149 (MQ=255)
                                                                        aTTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGttttt                                                                      >  3:657745/1‑149 (MQ=255)
                                                                         ttGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGtttttt                                                                     >  1:768266/1‑149 (MQ=255)
                                                                           gtgCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGtttttttt                                                                   <  4:1751101/149‑1 (MQ=255)
                                                                           gtgCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGtttttttt                                                                   <  2:1742987/149‑1 (MQ=255)
                                                                           gtgCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGtttttttt                                                                   <  4:71198/149‑1 (MQ=255)
                                                                            tgCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTta                                                                  >  2:1796822/1‑149 (MQ=255)
                                                                                     agaAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATAtatttat                                                         <  1:1796822/149‑1 (MQ=255)
                                                                                      gaAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATAtatttatt                                                        <  4:657745/149‑1 (MQ=255)
                                                                                             gTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCa                                                                              >  4:1705341/1‑120 (MQ=255)
                                                                                             gTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGATTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCa                                                                              <  3:1705341/120‑1 (MQ=255)
                                                                                                 aaCGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGt                                             >  1:1154809/1‑149 (MQ=255)
                                                                                                           aCTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTaa                                                   <  2:1759810/133‑1 (MQ=255)
                                                                                                           aCTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTaa                                                   >  1:1759810/1‑133 (MQ=255)
                                                                                                                cTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaacaac                                                                                                                <  4:323370/67‑1 (MQ=255)
                                                                                                                cTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaacaac                                                                                                                >  3:323370/1‑67 (MQ=255)
                                                                                                                  aaCGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGt                            >  3:886472/1‑149 (MQ=255)
                                                                                                                        tttCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTtatat                                                               >  1:1770585/1‑108 (MQ=255)
                                                                                                                        tttCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTtatat                                                               <  2:1770585/108‑1 (MQ=255)
                                                                                                                              gatgCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGaaa                <  4:886472/149‑1 (MQ=255)
                                                                                                                              gatgCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGaaa                <  2:1154809/149‑1 (MQ=255)
                                                                                                                                tgCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGaaa                                                                                                 <  1:1197254/66‑1 (MQ=255)
                                                                                                                                tgCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGaaa                                                                                                 >  2:1197254/1‑66 (MQ=255)
                                                                                                                                tgCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGt              <  2:768266/149‑1 (MQ=255)
                                                                                                                                 gCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGaa                                                                                                  <  3:674332/64‑1 (MQ=255)
                                                                                                                                 gCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGaa                                                                                                  >  4:674332/1‑64 (MQ=255)
                                                                                                                                   atatCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAAc                                                                                 >  3:658879/1‑79 (MQ=255)
                                                                                                                                   atatCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAAc                                                                                 <  4:658879/79‑1 (MQ=255)
                                                                                                                                       cctccGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTta                                                                  <  4:107360/90‑1 (MQ=255)
                                                                                                                                       cctccGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTta                                                                  >  3:107360/1‑90 (MQ=255)
                                                                                                                                            ggTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAAt  >  4:21177/1‑149 (MQ=255)
                                                                                                                                                   aGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTtata                                                                >  1:1659345/1‑80 (MQ=255)
                                                                                                                                                   aGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTtata                                                                <  2:1659345/80‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
ACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAAT  >  NC_000913/1346617‑1346905

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