Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 305,155 C→A 8.8% intergenic (+13/+16) hslR → / ← PP_0251 heat shock protein Hsp15/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947305,1550CA 8.8% 179.5 / 6.0 57intergenic (+13/+16)hslR/PP_0251heat shock protein Hsp15/membrane protein
Reads supporting (aligned to +/- strand):  ref base C (23/29);  new base A (2/3);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CGAAGAAAGTGTGCGTCGCCGTGAACAGGCGGCTGAAATGCGTAAGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAGCGCATCAACGCCTCGCTGCCTTTGGCCGCCAGCGGCGTGTAGTAGCTCCAGCCCAGCGCCA  >  NC_002947/305010‑305301
                                                                                                                                                 |                                                                                                                                                  
cGAAGAAAGTGTGCGTCGCCGTGAACAGGCGGCTGAAATGCGTAAGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTgcg                                                                                                                                                   >  6:294214/1‑147 (MQ=255)
cGAAGAAAGTGTGCGTCGCCGTGAACAGGCGGCTGAAATGCGTAAGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTgcg                                                                                                                                                   >  7:148897/1‑147 (MQ=255)
 gaagaaAGTGTGCGTCGCCGTGAACAGGCGGCTGAAATGCGTAAGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAg                                                                                                                                                >  1:138847/1‑149 (MQ=255)
 gaagaaAGTGTGCGTCGCCGTGAACAGGCGGCTGAAATGCGTAAGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAg                                                                                                                                                >  7:38432/1‑149 (MQ=255)
    gaaAGTGTGCGTCGCCGTGAACAGGCGGCTGAAATGCGTAAGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGccc                                                                                                                                             <  6:240887/149‑1 (MQ=255)
        gtgtGCGTCGCCGTGAACAGGCGGCTGAAATGCGTAAGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCGCAGCAGCCCCatt                                                                                                                                         <  3:164033/149‑3 (MQ=255)
           tgCGTCGCCGTGAACAGGCGGCTGAAATGCGTAAGGCCGGGGCCATGGGCGTGTCCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCaaaaaa                                                                                                                                      <  6:207548/149‑1 (MQ=255)
            gCGTCGCCGTGAACAGGCGGCTGAAATGCGTAAGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAAt                                                                                                                                     <  7:291941/149‑1 (MQ=255)
                  ccGTGAACAGGCGGCTGAAATGCGTAAGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCGCAGCAGCCCCa                                                                                                                                           >  3:289555/1‑137 (MQ=255)
                  ccGTGAACAGGCGGCTGAAATGCGTAAGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCGCAGCAGCCCCa                                                                                                                                           <  4:289555/137‑1 (MQ=255)
                   cGTGAGCAGGCGGCTGAAATGCGTAAGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAgcg                                                                                                                               >  6:279948/1‑148 (MQ=255)
                   cGTGAACAGGCGGCTGAAATGCGTAAGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAgcgc                                                                                                                              <  2:291392/149‑1 (MQ=255)
                   cGTGAACAGGCGGCTGAAATGCGTAAGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAg                                                                                                                                 <  2:161934/146‑1 (MQ=255)
                   cGTGAACAGGCGGCTGAAATGCGTAAGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAg                                                                                                                                 >  1:161934/1‑146 (MQ=255)
                                cTGAAATGCGTAAGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCAAGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTc                                                                                                                 >  5:146870/1‑149 (MQ=255)
                                  gAAATGCGTAAGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCGCAGCAGCCCCatt                                                                                                                                         <  5:185886/123‑3 (MQ=255)
                                  gAAATGCGTAAGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCGCAGCAGCCCCatt                                                                                                                                         >  6:185886/1‑121 (MQ=255)
                                   aaaTGCGTAAGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTc                                                                                                              >  3:147530/1‑149 (MQ=255)
                                      tGCGTAAGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGa                                                                                                           <  4:147530/149‑1 (MQ=255)
                                          tAAGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCaa                                                                                                       <  5:279948/149‑1 (MQ=255)
                                           aaGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAAc                                                                                                      >  5:169118/1‑149 (MQ=255)
                                            aGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACGAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAAc                                                                                                      <  5:294214/148‑1 (MQ=255)
                                              gCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACgg                                                                                                    >  5:196087/1‑148 (MQ=255)
                                                 ggggCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTaa                                                                                                 >  8:248126/1‑148 (MQ=255)
                                                 ggggCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTaa                                                                                                 <  7:248126/148‑1 (MQ=255)
                                                      cATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGa                                                                                           <  1:37157/149‑1 (MQ=255)
                                                        tGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGa                                                                                           >  7:241994/1‑147 (MQ=255)
                                                          ggCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGAcag                                                                                        <  6:196087/148‑1 (MQ=255)
                                                            cGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTcc                                                                                                             >  5:139172/1‑125 (MQ=255)
                                                            cGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTcc                                                                                                             <  6:139172/125‑1 (MQ=255)
                                                                accaccGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTg                                                                                 <  8:38432/149‑1 (MQ=255)
                                                                  caccGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGaa                                                                               >  1:321777/1‑149 (MQ=255)
                                                                   accGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGaac                                                                              <  2:252550/149‑1 (MQ=255)
                                                                            cggccgACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCt                                                                     >  8:151408/1‑149 (MQ=255)
                                                                               ccgACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTc                                                                  >  1:159279/1‑149 (MQ=255)
                                                                                 gACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCa                                                                <  2:138847/149‑1 (MQ=255)
                                                                                   ccAAGAAGCAGCCGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAgc                                                              <  8:148897/149‑1 (MQ=255)
                                                                                    cAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGaa                                                                               >  4:67312/1‑131 (MQ=255)
                                                                                    cAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGaa                                                                               <  3:67312/131‑1 (MQ=255)
                                                                                    cAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAgcg                                                             >  4:199427/1‑149 (MQ=255)
                                                                                    cAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAgcg                                                             >  3:53507/1‑149 (MQ=255)
                                                                                       gaagCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCGACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAGCGCAt                                                          <  5:140661/149‑1 (MQ=255)
                                                                                           cagcGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAGCGCATCAAc                                                      <  3:199427/149‑1 (MQ=255)
                                                                                                     aGATTCACCAGCTGCATGGAAGTTACGAATGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCCACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAGCGCATCAACGCCTCGCTGc                                            <  7:84609/149‑1 (MQ=255)
                                                                                                       aTTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAGCGCATCAACGCCTCGCTGCCt                                          >  1:159417/1‑149 (MQ=255)
                                                                                                           aCCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAGCGCATCAACGCCTCGCTGCCTTTgg                                      <  4:53507/149‑1 (MQ=255)
                                                                                                             cAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAGCGCATCAACGCCTCGCTGCCTTTGgcc                                    <  7:151408/149‑1 (MQ=255)
                                                                                                                   cATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAGCGCATCAACGCCTCGCTGCCTTTGGCCGCCAgc                              <  2:321777/149‑1 (MQ=255)
                                                                                                                    aTGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAGCGCATCAACGCCTCGCTGCCTTTGGCCGCCAgcg                             >  6:130968/1‑149 (MQ=255)
                                                                                                                     tGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCt                                                                     >  7:202104/1‑108 (MQ=255)
                                                                                                                     tGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCt                                                                     <  8:202104/108‑1 (MQ=255)
                                                                                                                          gTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAGCGCATCAACGCCTCGCTGCCTTTGGCCGCCAGCGGCgtg                        >  3:59234/1‑148 (MQ=255)
                                                                                                                              cGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAGCGCATCAACGCCTCGCTGCCTTTGGCCGCCAGCGGCGTgtagta                   <  5:167076/149‑1 (MQ=255)
                                                                                                                              cGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAGCGCATCAACGCCTCGCTGCCTTTGGCCGCCAGCGGCGTgtagt                    <  6:169118/148‑1 (MQ=255)
                                                                                                                                aTTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAGCGCATCAACGCCTCGCTGCCTTTGGCCGCCAGCGGCgt                         >  7:55965/1‑141 (MQ=255)
                                                                                                                                aTTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAGCGCATCAACGCCTCGCTGCCTTTGGCCGCCAGCGGCgt                         <  8:55965/141‑1 (MQ=255)
                                                                                                                                   gAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAGCGCATCAACGCCTCGCTGCCTTTGGCCGCCAGCGGCGTGTAGTAGCTc               <  6:146870/148‑1 (MQ=255)
                                                                                                                                     aTGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAGCGCATCAACGCCTCGCTGCCTTTGGCCGCCAGCGGCGTGTAGTAGCTccag            >  5:298503/1‑149 (MQ=255)
                                                                                                                                          gctgctGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAGCGCATCAACGCCTCGCTGCCTTTGGCCGCCAGCGGCGTGTAGTAGCTccagcccag       <  2:159279/149‑1 (MQ=255)
                                                                                                                                               tgcGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAGCGCATCAACGCCTCGCTGCCTTTGGCCGCCAGCGGCGTGTAGTAGCTCCAGCCCAgcgc    >  1:300328/1‑147 (MQ=255)
                                                                                                                                                 cgcAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAGCGCATCAACGCCTCGCTGCCTTTGGCCGCCAGCGGCGTGTAGTAGCTCCAGCCCAGCGCCa  <  2:159417/147‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CGAAGAAAGTGTGCGTCGCCGTGAACAGGCGGCTGAAATGCGTAAGGCCGGGGCCATGGGCGTGACCACCGATGGCCGGCCGACCAAGAAGCAGCGGCGGCAGATTCACCAGCTGCATGGAAGTTACGATTGAATGGGGCTGCTGCGCAGCCCCAAAAAATCTAGCGCAACACCGCCATTCGTCCGACCAACGGTAACGTGGACAGCAGCTTGAACAACGGCGCTGTCCAGCGCATCAACGCCTCGCTGCCTTTGGCCGCCAGCGGCGTGTAGTAGCTCCAGCCCAGCGCCA  >  NC_002947/305010‑305301

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