Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,346,770 G→C 8.2% intergenic (‑28/‑26) 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,7700GC 8.2% 293.5 / 16.0 85intergenic (‑28/‑26)gmr/ymiBcyclic‑di‑GMP phosphodiesterase; csgD regulator; modulator of RNase II stability/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (34/44);  new base C (3/4);  total (37/48)
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.94e-01

GGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGG  >  NC_000913/1346624‑1346913
                                                                                                                                                  |                                                                                                                                               
ggTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTc                                                                                                                                               <  4:1000909/149‑1 (MQ=255)
    cTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCcgcg                                                                                                                                           <  3:824812/149‑1 (MQ=255)
     tgtgGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCgg                                                                                                                                          <  1:168505/149‑1 (MQ=255)
     tgtgGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCgg                                                                                                                                          >  2:1137052/1‑149 (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                                                                                                                                <  2:449096/149‑1 (MQ=255)
               gtggtAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATaaaa                                                                                                                                <  1:1137052/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                                                                                                                                     >  1:1297441/1‑124 (MQ=255)
                                   aTCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaa                                                                                                                                     <  2:1297441/124‑1 (MQ=255)
                                      gctgctTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGaaa                                                                                                         <  3:950404/149‑1 (MQ=255)
                                        tgctTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGaaa                                                                                                         <  4:425872/147‑1 (MQ=255)
                                        tgctTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGaaa                                                                                                         >  3:425872/1‑147 (MQ=255)
                                           tttCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAAc                                                                                                    <  4:1416459/149‑1 (MQ=255)
                                            ttCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTAGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGTTTAAATAAAAGGAACAACATGACTTACATGAAATTAAcg                                                                                                   >  2:1606574/1‑149 (MQ=255)
                                                       gCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATg                                                                                                                     <  2:1543089/120‑1 (MQ=255)
                                                       gCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATg                                                                                                                     >  1:1543089/1‑120 (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:657745/1‑149 (MQ=255)
                                                                 aTTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGttttt                                                                              >  3:71198/1‑149 (MQ=255)
                                                                  ttGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGtttttt                                                                             >  1:768266/1‑149 (MQ=255)
                                                                    gtgCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGtttttttt                                                                           <  2:1742987/149‑1 (MQ=255)
                                                                    gtgCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGtttttttt                                                                           <  4:71198/149‑1 (MQ=255)
                                                                    gtgCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGtttttttt                                                                           <  4:1751101/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                                                                       <  2:1770585/108‑1 (MQ=255)
                                                                                                                 tttCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTtatat                                                                       >  1:1770585/1‑108 (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                                                                        <  2:1659345/80‑1 (MQ=255)
                                                                                                                                            aGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTtata                                                                        >  1:1659345/1‑80 (MQ=255)
                                                                                                                                              cGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTgg  <  4:1440835/148‑1 (MQ=255)
                                                                                                                                              cGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTgg  >  3:1440835/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGG  >  NC_000913/1346624‑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: