Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,535,981 T→G 6.4% L333R (CTG→CGG) ‡ pncB → nicotinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,535,9810TG 6.4% 185.0 / 3.4 62L333R (CTG→CGG) ‡pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (26/32);  new base G (2/2);  total (28/34)
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.56e-01

CCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAGCCGCCAAGACC  >  NC_002947/5535835‑5536126
                                                                                                                                                   |                                                                                                                                                 
ccTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTgc                                                                                                                                                 <  6:22662/149‑1 (MQ=255)
ccTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTgc                                                                                                                                                 <  4:283968/149‑1 (MQ=255)
 cTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTgcg                                                                                                                                                >  4:253329/1‑149 (MQ=255)
     gcCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGa                                                                                                                                            <  8:248629/149‑1 (MQ=255)
      ccATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGAc                                                                                                                                           >  2:122665/1‑149 (MQ=255)
       cATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCAGGGCGCGGaaa                                                                                                                                          <  4:351198/149‑3 (MQ=255)
       cATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCAGGGCGCGGaa                                                                                                                                           >  3:70646/1‑147 (MQ=255)
           aTTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCAGGGCGCGGaaaatt                                                                                                                                       <  4:70646/148‑6 (MQ=255)
              cGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCGATTGGCGATTACCAGATACTGGGTATCGACCCGATGACCAA‑GCCGCTGGGGTTTTCCGGCGGCGTGAACCTGACCCGCTCGCTGGAGATGTTCCGCGCCCTGCGCGGCCGGATCCAc                                                                                                                                   >  4:58454/1‑149 (MQ=255)
                 gCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTc                                                                                                                                <  4:83012/149‑1 (MQ=255)
                     gCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCt                                                                                                                            <  6:124793/149‑1 (MQ=255)
                       cGGTAGCCTGGGCGGAGAAGGCCAATTCCCATTAGCAGAAACTGGGTATCGGCCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCttt                                                                                                                          <  4:85461/149‑1 (MQ=255)
                        ggTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTg                                                                                                                         >  7:205227/1‑149 (MQ=255)
                        ggTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTg                                                                                                                         <  6:2969/149‑1 (MQ=255)
                              cTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCg                                                                                                                   <  6:243879/149‑1 (MQ=255)
                                 ggCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCgg                                                                                                                  <  1:292086/147‑1 (MQ=255)
                                 ggCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCgg                                                                                                                  <  3:53577/147‑1 (MQ=255)
                                      agaAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGGACCCTGGTGTTTTCCGACGGCCACAACCTGACCCGCACGCAGGAAATTTACCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCaa                                                                                                            >  1:283874/1‑149 (MQ=255)
                                      agaAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAac                                                                                                           >  1:111533/1‑149 (MQ=255)
                                      agaAGGCCATTGCCCATTACCACAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAac                                                                                                           <  7:247422/149‑1 (MQ=255)
                                        aaGGCCATTGCACATTACCAGAAACTGGGTAGCGACCCGATGACCAA‑GACGCTGGGGTGTTCCGACGGCCTCCACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAacct                                                                                                         <  7:333276/149‑1 (MQ=255)
                                          ggCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAacctga                                                                                                       >  6:37542/1‑149 (MQ=255)
                                              aTTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAacctgacctg                                                                                                   <  8:203893/149‑1 (MQ=255)
                                              aTTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAacctgacctg                                                                                                   <  6:87404/149‑1 (MQ=255)
                                                  cccATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGcaa                                                                                                <  2:111533/148‑3 (MQ=255)
                                                          cAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAgg                                                                                       >  5:298078/1‑149 (MQ=255)
                                                          cAGAAACTGGGTATCGAACCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTCAACAGCTCGCTCGAAATTTTCCGCGCCCTGAGCGGACGCATCAACGTCAGATTTGGCATCGGCACCAAACTCACCTGCGACATCCCAgg                                                                                       >  4:331668/1‑149 (MQ=255)
                                                           aGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGgt                                                                                      <  7:124713/149‑1 (MQ=255)
                                                           aGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGgt                                                                                      >  5:89826/1‑149 (MQ=255)
                                                              aaCTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGgtgt                                                                                    <  2:316453/148‑1 (MQ=255)
                                                                   ggTAACGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCCCAGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCATAGGCAACGGCACCAACCAGACCTGCGACATCCCAGGTGTGcccac                                                                              >  5:170433/1‑147 (MQ=255)
                                                                          aCCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAAc                                                                       <  1:346450/149‑1 (MQ=255)
                                                                            ccGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACa                                                                      >  1:125216/1‑148 (MQ=255)
                                                                              gATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGAGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCg                                                                   >  8:22091/1‑149 (MQ=255)
                                                                                   ccAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCt                                                               >  1:203433/1‑148 (MQ=255)
                                                                                      a‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTaaa                                                           <  3:58454/149‑1 (MQ=255)
                                                                                        gACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTTGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTaaaa                                                          <  2:110683/149‑1 (MQ=255)
                                                                                           gCTGGTGTTTTCCGACGGCCTCAACCTGTCCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCTTCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCGCATGAACATCGTGCTTAAAATGa                                                       <  2:283874/149‑1 (MQ=255)
                                                                                             tGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGGGATGTCGCAGgtcagg                                                                                  >  6:68339/1‑116 (MQ=255)
                                                                                                 gtTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCg                                                                                                 >  5:9606/1‑101 (MQ=255)
                                                                                                 gtTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCg                                                                                                 <  6:9606/101‑1 (MQ=255)
                                                                                                  ttttCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTg                                                >  5:308015/1‑149 (MQ=255)
                                                                                                       cGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAAc                                            <  6:298078/148‑1 (MQ=255)
                                                                                                              cTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATccc                                                                                          >  2:123825/1‑95 (MQ=255)
                                                                                                              cTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACg                                           >  3:26268/1‑142 (MQ=255)
                                                                                                              cTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACg                                           <  4:26268/142‑1 (MQ=255)
                                                                                                               tCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCg                                   >  5:168773/1‑149 (MQ=255)
                                                                                                                 aaCCTGACCGGCTCTGTGGTATTTTTTCGCGACCTGGGCTGACGCATTATGTTCAGCTTTGTCATCTGTACTTAACTGACCTGCTAAATCCATTGTGTGTGCCCGATGAACATCGTGTTTTACCTGACCGACTGCAACGGCGCGCCGGt                                 <  6:380137/149‑1 (MQ=255)
                                                                                                                  aCCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTg                                <  7:22091/149‑1 (MQ=255)
                                                                                                                    cTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGc                              <  8:205227/149‑1 (MQ=255)
                                                                                                                       aCCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGcc                             >  7:146280/1‑147 (MQ=255)
                                                                                                                         ccGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGCCGCATCACCGTCAGCTTTGGCATCGGCACCCACCTGACCTGCGCCATCCCAGGTGTGGCCCCGATGACCTGCGTGCGGAAAGTGCCCGCCTGCAACGGCGCGCCGGGGGcc                             >  7:317578/1‑145 (MQ=255)
                                                                                                                         ccGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCAAGAAAATCGTGCTTAAAATGACCGACTCCAACGGCGCGCCGGTGGCCaac                          >  8:296890/1‑147 (MQ=255)
                                                                                                                               gctGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATctct                    <  2:125216/148‑1 (MQ=255)
                                                                                                                               gctGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATctct                    <  5:37542/148‑1 (MQ=255)
                                                                                                                               gctGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTg                   <  6:89826/149‑1 (MQ=255)
                                                                                                                                 tGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGAc                 >  4:89293/1‑149 (MQ=255)
                                                                                                                                         gtCCGCGCAGGGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCCTGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAgccgc         >  5:329071/2‑149 (MQ=255)
                                                                                                                                               gcCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGTCACCAACCTGACCTGCGACATCCCAGGTGTGGACCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAGCCGCCAAGAc   <  3:331668/149‑1 (MQ=255)
                                                                                                                                                 ccTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGgcg                                        <  7:4783/110‑1 (MQ=255)
                                                                                                                                                 ccTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGgcg                                        >  8:4783/1‑110 (MQ=255)
                                                                                                                                                  cTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGCTCTCTGACGAAGCCGCCAAGAcc  >  6:164441/1‑147 (MQ=255)
                                                                                                                                                   tGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAgccgc         <  7:340836/139‑1 (MQ=255)
                                                                                                                                                   tGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAgccgc         >  8:340836/1‑139 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCTGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAA‑GACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAGCCGCCAAGACC  >  NC_002947/5535835‑5536126

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