Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 308,739 A→G 7.0% intergenic (+1754/‑115) hslO → / → PP_0254 molecular chaperone HslO/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947308,7390AG 7.0% 164.3 / 2.6 58intergenic (+1754/‑115)hslO/PP_0254molecular chaperone HslO/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/28);  new base G (2/2);  total (27/31)
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.84e-01

GGCTGACAAAGCGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGAGCCAGCCAGATGACCGAAGCCCCTCAGCG  >  NC_002947/308591‑308873
                                                                                                                                                     |                                                                                                                                       
ggCTGACAAAGCGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCa                                                                                                                                         <  7:71519/149‑1 (MQ=255)
   tGACAAAGCGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAggg                                                                                                                                     >  3:38267/1‑149 (MQ=255)
     aCAAAGCGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAgggg                                                                                                                                    <  1:162502/148‑1 (MQ=255)
     aCAAAGCGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAggg                                                                                                                                     <  8:110666/147‑1 (MQ=255)
        aaGCGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGggcggc                                                                                                                                <  8:62261/149‑1 (MQ=255)
         aGCGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCTGGGGGCGGCt                                                                                                                               <  4:147058/149‑1 (MQ=255)
           cGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCttt                                                                                                                             <  2:33757/149‑1 (MQ=255)
                      gaGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGt                                                                                                                    >  3:141602/1‑147 (MQ=255)
                             ccAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCtt                                                                                                            <  1:45685/148‑1 (MQ=255)
                               aaaGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGcc                                                                                                         >  3:108873/1‑149 (MQ=255)
                               aaaGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGcc                                                                                                         >  3:125246/1‑149 (MQ=255)
                                    tCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGa                                                                                                    <  8:22642/149‑1 (MQ=255)
                                        gCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAg                                                                                                                                       >  4:238519/1‑110 (MQ=255)
                                         ccAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTgg                                                                                                                      >  6:33456/1‑126 (MQ=255)
                                          cAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCGCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGaa                                                                                                   >  1:56434/1‑145 (MQ=255)
                                          cAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGaa                                                                                                   <  2:56434/144‑1 (MQ=255)
                                          cAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGc                                                                                              >  2:51488/1‑149 (MQ=255)
                                              cTGTTCATCGAAAACTTCAACAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCg                                                                                           <  1:51488/149‑1 (MQ=255)
                                                   cATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAg                                                                                      <  4:141602/148‑1 (MQ=255)
                                                              tCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGa                                                                               <  8:84876/144‑1 (MQ=255)
                                                              tCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCAGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGa                                                                               >  7:84876/1‑144 (MQ=255)
                                                                aagaagTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCTACGAGAAAGCCGCGCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGTGGCCCGCATAGGCAGCGAAAGTCTg                                                                        <  7:31142/149‑1 (MQ=255)
                                                                aagaagTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAg                                                                                      >  2:66428/1‑135 (MQ=255)
                                                                aagaagTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAACCCGCCCCCTGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAg                                                                                      <  1:66428/135‑1 (MQ=255)
                                                                         gAAGTTTTCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCt                                                                <  2:4475/148‑1 (MQ=255)
                                                                                   aCGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTc                                                     <  6:143234/149‑1 (MQ=255)
                                                                                    cGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCa                                                    <  1:124777/149‑1 (MQ=255)
                                                                                    cGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCa                                                    <  5:276504/149‑1 (MQ=255)
                                                                                         aTCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCa                                                 >  6:80788/1‑148 (MQ=255)
                                                                                             aGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCTGCATCCAt                                            <  4:125246/149‑1 (MQ=255)
                                                                                             aGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCAt                                            >  6:246977/1‑149 (MQ=255)
                                                                                                ccGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCACGCATAGGCAGCGAAAGTGTGGTTTTTCTGTTTCACTGTCAGCATcctgtgc                                         >  4:9179/1‑144 (MQ=255)
                                                                                                 cGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGt                                                                      <  8:207276/119‑1 (MQ=255)
                                                                                                 cGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGt                                                                      >  7:207276/1‑119 (MQ=255)
                                                                                                  gCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATcc                                              <  3:178646/142‑1 (MQ=255)
                                                                                                  gCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATcc                                              >  4:178646/1‑142 (MQ=255)
                                                                                                      gCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGAc                                   <  5:80788/149‑1 (MQ=255)
                                                                                                          gCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGaa                                      >  1:25511/1‑142 (MQ=255)
                                                                                                          gCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGaa                                      <  2:25511/142‑1 (MQ=255)
                                                                                                           cAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACaa                                 >  8:78158/1‑146 (MQ=255)
                                                                                                           cAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAaga                               <  7:78158/148‑1 (MQ=255)
                                                                                                            aGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCa                                                    <  5:84706/126‑1 (MQ=255)
                                                                                                            aGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCa                                                    >  6:84706/1‑126 (MQ=255)
                                                                                                              cTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCATGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGagcca                           <  3:173733/149‑1 (MQ=255)
                                                                                                              cTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGTCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGa                                       <  7:224802/137‑1 (MQ=255)
                                                                                                              cTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGa                                       >  8:224802/1‑137 (MQ=255)
                                                                                                                   aGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTGTTTTGTGGGTGAGGGCTTGCCCGGGAAGGGGGCCGCATAGGCAGCGAAGGGCTGGTTTTTGTGTTTCACTGGGCGCCTCGAGTGCGAAGCCGCGagccaggctg                      >  8:123037/1‑145 (MQ=255)
                                                                                                                    gCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGagc                             <  4:155932/141‑1 (MQ=255)
                                                                                                                    gCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGagc                             >  3:155932/1‑141 (MQ=255)
                                                                                                                     ccGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAagag                              <  3:113001/139‑1 (MQ=255)
                                                                                                                     ccGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTATGTTTCACTTTAAGCATCCATTGCGAAGACAagag                              >  4:113001/1‑139 (MQ=255)
                                                                                                                                 gCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGAGCCAGCCAGATGAc                 >  5:89527/1‑140 (MQ=255)
                                                                                                                                 gCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGAGCCAGCCAGATGAc                 <  6:89527/140‑1 (MQ=255)
                                                                                                                                  ccacaaaaaaGCCGCCCCTGGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGAGCCAGCCAGATGACCGAAGCCCCt       >  7:97806/8‑149 (MQ=255)
                                                                                                                                    acaaaaaaGCCGCCCCTGGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGaa                                      >  5:222490/6‑116 (MQ=255)
                                                                                                                                    acaaaaaaGCCGCCCCTGGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGaa                                      <  6:222490/111‑1 (MQ=255)
                                                                                                                                     cGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGAGCCAGCCAGa                     >  5:263991/1‑132 (MQ=255)
                                                                                                                                     cGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGAGCCAGCCAGa                     <  6:263991/132‑1 (MQ=255)
                                                                                                                                        gaAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGCCAAGAGCCAGCCAGATGACCGAAGCCCCTCAgcg  >  3:218739/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                       
GGCTGACAAAGCGGCCTACGACGAGGCTGCCAAAGGTCTGGCCAAGCTGTTCATCGAAAACTTCAAGAAGTTCGAAGTTTCCGACGCCATCAAGGCCGCTGGCCCGCAGCTGTAAGCCGAGCTCAGCCGCAACGAGAAAGCCGCCCCCAGGGGCGGCTTTTTTGTGGGTTAGGGCTTGCCCGCGAAGGGGCCCGCATAGGCAGCGAAAGTCTGGTTTTTCTGTTTCACTGTCAGCATCCATTGCGAAGACAAGAGCCAGCCAGATGACCGAAGCCCCTCAGCG  >  NC_002947/308591‑308873

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: