Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,346,770 G→C 5.3% 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 5.3% 267.5 / 5.4 76intergenic (‑28/‑26)gmr/ymiBcyclic‑di‑GMP phosphodiesterase; csgD regulator; modulator of RNase II stability/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (32/40);  new base C (2/2);  total (34/42)
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.18e-01

GAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCG  >  NC_000913/1346516‑1346992
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                              
gAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCACTGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGa                                                                                                                                                                                                            >  2:343583/1‑275 (MQ=255)
        gCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaacaac                                                                                                                                                                                                       <  2:40293/272‑1 (MQ=255)
               tCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACa                                                                                                                                                                                             <  2:84251/275‑1 (MQ=255)
                ggACGTTCAGACTCATCATCAGGCTGGAGGTGATGACCGTAATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACAt                                                                                                                                                                                            <  1:55364/274‑1 (MQ=255)
                cGACGGTCAGACTCATCATCAGGCGTGACGTGCTGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACAt                                                                                                                                                                                            <  4:418963/275‑1 (MQ=255)
                                                                                                            ggTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGt                                                                                                                                                                                                                               <  7:948358/148‑1 (MQ=255)
                                                                                                                   tgGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCGGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGATAAACGCATTGTTTGTTTTATATCTTTATTTGTAATCCAGTTTTGGGAAAAAGCCCGTTTTCAAACGAAAGTCAGTTAAAAAAt                                                                                         >  4:191959/1‑275 (MQ=255)
                                                                                                                   tgGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCgg                                                                                                                                                                                                                         <  8:617486/147‑1 (MQ=255)
                                                                                                                   tgGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCgg                                                                                                                                                                                                                         >  7:617486/1‑147 (MQ=255)
                                                                                                                       ttCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaa                                                                                                                                                                                                                    >  8:555437/1‑148 (MQ=255)
                                                                                                                           gtggtAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaataaa                                                                                                                                                                                                                <  5:394989/148‑1 (MQ=255)
                                                                                                                                  aaaaGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaataaa                                                                                                                                                                                                                <  5:299217/141‑1 (MQ=255)
                                                                                                                                  aaaaGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaataaa                                                                                                                                                                                                                >  6:299217/1‑141 (MQ=255)
                                                                                                                                    aaGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaataaa                                                                                                                                                                                                                >  8:293822/1‑139 (MQ=255)
                                                                                                                                    aaGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaataaa                                                                                                                                                                                                                <  7:293822/139‑1 (MQ=255)
                                                                                                                                         cAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaa                                                                                                                                                                                                                    >  6:486943/1‑130 (MQ=255)
                                                                                                                                         cAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaa                                                                                                                                                                                                                    <  5:486943/130‑1 (MQ=255)
                                                                                                                                          aaaaCATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACt                                                                                                                                                                                                 <  8:793104/148‑1 (MQ=255)
                                                                                                                                            aaCATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTa                                                                                                                                                                                               <  5:184495/148‑1 (MQ=255)
                                                                                                                                              cATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACa                                                                                                                                                                                             <  6:770642/148‑1 (MQ=255)
                                                                                                                                              cATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACa                                                                                                                                                                                             <  6:48138/148‑1 (MQ=255)
                                                                                                                                                 cGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACCGCGGACTCcgctg                                                                                                                                                                                                                         <  6:881946/117‑3 (MQ=255)
                                                                                                                                                 cGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACCGCGGACTCcgctg                                                                                                                                                                                                                         >  5:881946/1‑115 (MQ=255)
                                                                                                                                                 cGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCgggg                                                           >  2:406799/1‑275 (MQ=255)
                                                                                                                                                  gctgctTTCCGTCAACCGCCATGATGGATTTTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGa                                                          <  2:100281/275‑1 (MQ=255)
                                                                                                                                                  gctgctTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaaca                                                                                                                                                                                                         >  5:989849/1‑132 (MQ=255)
                                                                                                                                                  gctgctTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaaca                                                                                                                                                                                                         <  6:989849/132‑1 (MQ=255)
                                                                                                                                                   ctgctTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGt                                                                                                                                                                                                                        <  7:1137792/116‑1 (MQ=255)
                                                                                                                                                   ctgctTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGt                                                                                                                                                                                                                        >  8:1137792/1‑116 (MQ=255)
                                                                                                                                                   ctgctTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGa                                                                                                                                                                                                            <  5:548162/128‑1 (MQ=255)
                                                                                                                                                   ctgctTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGa                                                                                                                                                                                                            >  6:548162/1‑128 (MQ=255)
                                                                                                                                                        ttCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTaa                                                                                                                                                                                    <  5:860303/147‑1 (MQ=255)
                                                                                                                                                        ttCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTaa                                                                                                                                                                                    >  6:860303/1‑147 (MQ=255)
                                                                                                                                                                     cAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCcg                                                                                                                                                                                                                            <  5:833861/94‑1 (MQ=255)
                                                                                                                                                                     cAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCcg                                                                                                                                                                                                                            >  6:833861/1‑94 (MQ=255)
                                                                                                                                                                     cAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGtt                                                                                                                                                                                                                       >  7:1113215/1‑99 (MQ=255)
                                                                                                                                                                     cAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGtt                                                                                                                                                                                                                       <  8:1113215/99‑1 (MQ=255)
                                                                                                                                                                               tgtgCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGc                             >  2:25627/1‑275 (MQ=255)
                                                                                                                                                                                gtgCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGttttttt                                                                                                                                                           <  7:555437/148‑1 (MQ=255)
                                                                                                                                                                                     aTCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaataaa                                                                                                                                                                                                                >  7:245689/1‑90 (MQ=255)
                                                                                                                                                                                     aTCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaataaa                                                                                                                                                                                                                <  8:245689/90‑1 (MQ=255)
                                                                                                                                                                                     aTCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCAt                                                                                                                                                                    <  8:764826/134‑1 (MQ=255)
                                                                                                                                                                                     aTCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCAt                                                                                                                                                                    >  7:764826/1‑134 (MQ=255)
                                                                                                                                                                                          agaAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACCGCGGACTCCGCTGTta                                                                                                                                                                                                                      <  8:1122057/79‑1 (MQ=38)
                                                                                                                                                                                          agaAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACCGCGGACTCCGCTGTta                                                                                                                                                                                                                      >  7:1122057/1‑79 (MQ=39)
                                                                                                                                                                                          agaAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTta                                                                                                                                                         <  5:1086797/140‑1 (MQ=255)
                                                                                                                                                                                          agaAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTta                                                                                                                                                         >  6:1086797/1‑140 (MQ=255)
                                                                                                                                                                                            aaaGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGa                                                                                                                                                                                                            <  5:591533/87‑1 (MQ=255)
                                                                                                                                                                                            aaaGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGa                                                                                                                                                                                                            >  6:591533/1‑87 (MQ=255)
                                                                                                                                                                                                ttGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGaa                                                                                                                                                                                         <  8:874535/102‑1 (MQ=255)
                                                                                                                                                                                                ttGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGaa                                                                                                                                                                                         >  7:874535/1‑102 (MQ=255)
                                                                                                                                                                                                     cAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaa                                                                                                                                                                                                                     >  5:395712/1‑69 (MQ=255)
                                                                                                                                                                                                     cAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaa                                                                                                                                                                                                                     <  6:395712/69‑1 (MQ=255)
                                                                                                                                                                                                       aCGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGaaa                                                                                                                                                                                        >  6:532594/1‑96 (MQ=255)
                                                                                                                                                                                                       aCGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGaaa                                                                                                                                                                                        <  7:925222/96‑1 (MQ=255)
                                                                                                                                                                                                       aCGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGaaa                                                                                                                                                                                        <  5:532594/96‑1 (MQ=255)
                                                                                                                                                                                                       aCGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGaaa                                                                                                                                                                                        >  8:925222/1‑96 (MQ=255)
                                                                                                                                                                                                        cGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGaa                                                                                                                                                                                         <  5:493786/94‑1 (MQ=255)
                                                                                                                                                                                                        cGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGaa                                                                                                                                                                                         >  6:493786/1‑94 (MQ=255)
                                                                                                                                                                                                           tgtgGACTCCCTAACGGTTTTCATGATGCAGATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGgcg  <  4:406798/274‑1 (MQ=255)
                                                                                                                                                                                                                      tAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAacaaca                                                                                                                                                                                                      <  5:751747/67‑1 (MQ=255)
                                                                                                                                                                                                                      tAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAacaaca                                                                                                                                                                                                      >  6:751747/1‑67 (MQ=255)
                                                                                                                                                                                                                      tAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAAcgg                                                                                                                                                                                 >  5:551383/1‑88 (MQ=255)
                                                                                                                                                                                                                      tAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAAcgg                                                                                                                                                                                 <  6:551383/88‑1 (MQ=255)
                                                                                                                                                                                                                       aaCGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCa                                                                                                                                                                     >  5:479860/1‑99 (MQ=255)
                                                                                                                                                                                                                       aaCGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCa                                                                                                                                                                     <  6:479860/99‑1 (MQ=255)
                                                                                                                                                                                                                                  tgatgCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCAtt                                                                                                                                                                   <  8:953799/90‑1 (MQ=255)
                                                                                                                                                                                                                                  tgatgCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCAtt                                                                                                                                                                   >  7:953799/1‑90 (MQ=255)
                                                                                                                                                                                                                                        atatCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATAtattta                                                                                                                                                 <  6:544020/102‑1 (MQ=255)
                                                                                                                                                                                                                                        atatCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATAtattta                                                                                                                                                 >  5:544020/1‑102 (MQ=255)
                                                                                                                                                                                                                                        atatCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATAtattt                                                                                                                                                  >  6:943854/1‑101 (MQ=255)
                                                                                                                                                                                                                                        atatCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATAtattt                                                                                                                                                  <  5:943854/101‑1 (MQ=255)
                                                                                                                                                                                                                                             ctccGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCAt                                                                                                                                                                    >  8:426373/1‑78 (MQ=255)
                                                                                                                                                                                                                                             ctccGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCAt                                                                                                                                                                    <  7:426373/78‑1 (MQ=255)
                                                                                                                                                                                                                                              tccGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGt                                                                                                                                                                 >  6:694379/1‑80 (MQ=255)
                                                                                                                                                                                                                                              tccGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGt                                                                                                                                                                 <  5:694379/80‑1 (MQ=255)
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                              
GAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAGCAGCCTGTTCGGCAGATAACTGCAATGTACGATCAGGTTCTGTGGTTTCGGTGGTAGAAAAGCGCAAAACATCGCTGCTTTCCGTCAACCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAAACGAAAGTCAGTTAAAAAATCTGCCTGGATATAACGAAGGTAGAGCGGGGAAATAAACGGCCCATCCATGAGGAATGGGCCGTGAAAGGAGATTACGCGACCGGGCG  >  NC_000913/1346516‑1346992

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