Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,346,776 G→T 8.6% intergenic (‑34/‑20) 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,7760GT 8.6% 277.8 / 14.6 81intergenic (‑34/‑20)gmr/ymiBcyclic‑di‑GMP phosphodiesterase; csgD regulator; modulator of RNase II stability/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (33/41);  new base T (3/4);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

CTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGG  >  NC_000913/1346628‑1346913
                                                                                                                                                    |                                                                                                                                         
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                                                                                                                                        <  3:1254526/141‑1 (MQ=255)
           gtggtAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACCGCGGACTCCGCTGtt                                                                                                                                        >  4:1254526/1‑141 (MQ=255)
           gtggtAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATaaaa                                                                                                                                <  1:1137052/149‑1 (MQ=255)
           gtggtAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATaaaa                                                                                                                                <  2:449096/149‑1 (MQ=255)
           gtggtAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATaaaa                                                                                                                                <  2:832774/149‑1 (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                                                                                                                                >  4:632011/1‑138 (MQ=255)
                      gcgcAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATaaaa                                                                                                                                <  3:632011/138‑1 (MQ=255)
                      gcgcAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATg                                                                                                                     >  3:1416459/1‑149 (MQ=255)
                       cgcAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGa                                                                                                                    <  1:1112395/149‑1 (MQ=255)
                         cAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACt                                                                                                                  <  4:1539023/149‑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                                                                                                                                     <  2:1297441/124‑1 (MQ=255)
                               aTCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaa                                                                                                                                     >  1:1297441/1‑124 (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)
                                                   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                                                                                                                        <  4:366070/111‑1 (MQ=255)
                                                         aTGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaacaac                                                                                                                        >  3:366070/1‑111 (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                                                                           <  2:1742987/149‑1 (MQ=255)
                                                                gtgCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGtttttttt                                                                           <  4:1751101/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                                                           >  1:1759810/1‑133 (MQ=255)
                                                                                                aCTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTaa                                                           <  2:1759810/133‑1 (MQ=255)
                                                                                                     cTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaacaac                                                                                                                        >  3:323370/1‑67 (MQ=255)
                                                                                                     cTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaacaac                                                                                                                        <  4:323370/67‑1 (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)
                                                                                                                                          cGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTgg  <  4:1440835/148‑1 (MQ=255)
                                                                                                                                          cGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTgg  >  3:1440835/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                         
CTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGG  >  NC_000913/1346628‑1346913

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