Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 307,866 +C intergenic (+881/‑988) 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_002947307,8651.C96.6% 192.4 / 0.0 59intergenic (+880/‑989)hslO/PP_0254molecular chaperone HslO/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (27/30);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

ATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAGTA  >  NC_002947/307723‑308006
                                                                                                                                                |                                                                                                                                                
aTGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCTCTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGccgtta                                                                                                                                            >  2:51883/1‑144 (MQ=255)
 tGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTAc                                                                                                                                           <  7:98525/149‑1 (MQ=255)
      tcGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGAt                                                                                                                                      >  2:257396/1‑149 (MQ=255)
        ggTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGAt                                                                                                                                      <  2:313217/147‑1 (MQ=255)
            cAGAACTTCCTGCTGCCGGCTGGCAATGAAATGACTAAGCACGGCGCGCCCAAAAACGGCGAAGAAGGAGATGTGACATTGTTGTTCGGCCTGACCGGCACTGGAAAGAAAACCCTTACGGGTGACGAAAGCCGTTACCGGATCGGTga                                                                                                                                >  7:55085/1‑149 (MQ=255)
            cAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTga                                                                                                                                <  6:286534/149‑1 (MQ=255)
             aGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTgac                                                                                                                               <  7:160371/149‑1 (MQ=255)
                 cTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTgacgac                                                                                                                            >  7:297858/1‑148 (MQ=255)
                  ttCCTGCTTCCGGCTGCCGATGTACTTCCTATGCACTTCGCGTCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGAGTACCCTTTCGGCTGACGAAAGccgtt                                                                                                                                             <  8:219254/130‑5 (MQ=255)
                  ttCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGccgtt                                                                                                                                             >  7:219254/1‑126 (MQ=255)
                      tgctgcCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCAc                                                                                                                      >  1:197204/1‑149 (MQ=255)
                        ctgccggctgcctATGTACTGGCTATGCACTGCGCTGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACg                                                                                                                     <  3:44569/148‑1 (MQ=255)
                        ctgccggctgccgATGTACTGCCTATGAACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACg                                                                                                                     <  3:200307/148‑1 (MQ=255)
                         tgccggctgccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑Gt                                                                                                                 >  2:176081/1‑149 (MQ=255)
                         tgccggctgccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑Gt                                                                                                                 >  5:230999/1‑149 (MQ=255)
                             ggctgccgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTggg                                                                                                             <  6:39477/149‑1 (MQ=255)
                                   cgATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAg                                                                                                       <  1:137386/149‑1 (MQ=255)
                                        tACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGtt                                                                                                  <  4:256610/149‑1 (MQ=255)
                                        tACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGtt                                                                                                  <  1:257396/149‑1 (MQ=255)
                                         aCTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTg                                                                                                 >  5:256417/1‑149 (MQ=255)
                                          cTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCCAGACTACCCTTTCGGCTGACGAAAGCAGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGt                                                                                                >  3:232093/1‑149 (MQ=255)
                                          cTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGt                                                                                                >  5:141949/1‑149 (MQ=255)
                                          cTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGt                                                                                                >  8:329939/1‑149 (MQ=255)
                                             ccTATGCTCTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCtt                                                                                             <  1:51883/149‑1 (MQ=255)
                                             ccTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCt                                                                                              <  2:339131/148‑1 (MQ=255)
                                             ccTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCt                                                                                              >  1:339131/1‑148 (MQ=255)
                                              cTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTc                                                                                            <  8:277046/149‑1 (MQ=255)
                                              cTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTc                                                                                            <  7:286884/149‑1 (MQ=255)
                                                 tGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCaa                                                                                          <  8:297858/148‑1 (MQ=255)
                                                    aCTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGCCGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATg                                                                                      >  5:306838/1‑149 (MQ=255)
                                                           ggCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAggcg                                                                               >  3:355306/1‑149 (MQ=255)
                                                           ggCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAggcg                                                                               >  6:270174/1‑149 (MQ=255)
                                                              cAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGgct                                                                            <  5:219147/149‑1 (MQ=255)
                                                                 cATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCAGGGAGTTGGGGCGAAGGCGTTGTATTCAAAATGCAAGGCGgct                                                                            >  1:120442/1‑148 (MQ=255)
                                                                     ggCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACg                                                                      <  3:38192/148‑1 (MQ=255)
                                                                      gCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGcc                                                                    <  5:270174/149‑1 (MQ=255)
                                                                            aagGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTg                                                               >  5:28126/1‑148 (MQ=255)
                                                                            aagGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTg                                                               >  8:250675/1‑148 (MQ=255)
                                                                             agGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGTTCGCTGACGACGTGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGTCGGCTGCTACGCCAAGTGCa                                                             <  8:55085/149‑1 (MQ=255)
                                                                             agGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCa                                                             <  2:197204/149‑1 (MQ=255)
                                                                              ggCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACACCAAGTGCa                                                             >  3:286681/1‑148 (MQ=255)
                                                                                  aTGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGAc                                                        <  4:355306/149‑1 (MQ=255)
                                                                                       aCCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGt                                                    <  1:176081/148‑1 (MQ=255)
                                                                                       aCCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGt                                                    <  6:230999/148‑1 (MQ=255)
                                                                                         cTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTcc                                                  >  5:28077/1‑148 (MQ=255)
                                                                                          ttGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTcc                                                  >  1:316779/1‑147 (MQ=255)
                                                                                            gTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCgag                                               <  6:141949/148‑1 (MQ=255)
                                                                                                ttcGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGtt                                                                                                  >  7:22166/1‑93 (MQ=255)
                                                                                                ttcGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGtt                                                                                                  <  8:22166/93‑1 (MQ=255)
                                                                                                 tcGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAACCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAAc                                         >  7:20946/1‑149 (MQ=255)
                                                                                                    gCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCAGTCCGAgaagaa                                          >  8:322774/1‑145 (MQ=255)
                                                                                                    gCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCAGTCCGAgaagaa                                          <  7:322774/145‑1 (MQ=255)
                                                                                                     ccTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGt                                                                                                   <  8:180032/87‑1 (MQ=255)
                                                                                                     ccTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGt                                                                                                   >  7:180032/1‑87 (MQ=255)
                                                                                                         tCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCgg                                  <  6:28077/148‑1 (MQ=255)
                                                                                                         tCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCg                                   <  2:316779/147‑1 (MQ=255)
                                                                                                                        aCTACCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGCAGACCGAGCCGGTCCTCTCGAAAGCCAt                  >  2:107198/1‑149 (MQ=255)
                                                                                                                           aCCCTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATc                 >  6:142414/1‑147 (MQ=255)
                                                                                                                           aCCCTTTCGGCTGACGAAAGCCGATACCTGATCGGTGACGACGAGCACG‑‑GTAGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATc                 >  6:136877/1‑147 (MQ=255)
                                                                                                                              cTTTCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAAAGAGCAGGTCATCTGGAAAGCCATCAAGCa            >  2:236328/1‑149 (MQ=255)
                                                                                                                               tttCGGCTGACGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCAc           >  1:228367/1‑149 (MQ=255)
                                                                                                                                        aCGAAAGCCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGGATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAGTa  <  8:20946/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
ATGTTCTCGGTGCAGAACTTCCTGCTGCCGGCTGCCGATGTACTGCCTATGCACTGCGCGGCCAACATCGGCGAAGAAGGCGATGTGACCTTGTTCTTCGGCCTGTCCGGCACTGGCAAGACTACCCTTTCGGCTGACGAAAGCGTTACCTGATCGGTGACGACGAGCACG‑‑GTTGGGGCGAAGGCGTTGTCTTCAACATGGAAGGCGGCTGCTACGCCAAGTGCATCGACCTGTCCGAGAAGAACGAGCCGGTCATCTGGAAAGCCATCAAGCACGGCGCAGTA  >  NC_002947/307723‑308006

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: