Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,762,926 T→C 8.2% F57L (TTC→CTC)  PP_3326 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,762,9260TC 8.2% 164.9 / 4.7 61F57L (TTC→CTC) PP_3326hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/32);  new base C (2/3);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

CTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAA  >  NC_002947/3762780‑3763063
                                                                                                                                                   |                                                                                                                                          
cTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTc                                                                                                                                          <  1:11054/149‑1 (MQ=255)
 ttCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTc                                                                                                                                          >  1:243540/1‑148 (MQ=255)
   cccTGTTGGCCCTGATGGTCAGTGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑c                                                                                                                                      >  4:163378/1‑149 (MQ=255)
   cccTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACGGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGa                                                                                                                                        <  1:282051/148‑1 (MQ=255)
   cccTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCg                                                                                                                                         <  2:239928/147‑1 (MQ=255)
    ccTGTTGGCCCTGATGGTCAGCGTCAACGCGGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑Cg                                                                                                                                     >  5:48784/1‑149 (MQ=255)
    ccTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑c                                                                                                                                      <  2:120121/148‑1 (MQ=255)
     cTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGt                                                                                                                                    >  5:126100/1‑149 (MQ=255)
        ttGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGa                                                                                                                                  <  6:342236/148‑1 (MQ=255)
               tGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGaa                                                                                                                           >  3:126767/1‑148 (MQ=255)
                   ggTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATg                                                                                                                       >  1:251244/1‑149 (MQ=255)
                     tCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATc                                                                                                                    <  8:166479/149‑1 (MQ=255)
                     tCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATc                                                                                                                    <  5:199009/149‑1 (MQ=255)
                     tCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATc                                                                                                                    >  6:332851/1‑149 (MQ=255)
                     tCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATc                                                                                                                    <  6:126100/149‑1 (MQ=255)
                         cGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCgg                                                                                                                <  7:3884/149‑1 (MQ=255)
                             aaCGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCC‑AGGCGCTCGACCGTGACCGGGATCATGTTCCCGGTCa                                                                                                             <  1:241160/148‑1 (MQ=38)
                                   ggCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGa                                                                                                      <  2:114507/149‑1 (MQ=255)
                                      gcCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATGTTCCCGGTCACgtcgaacg                                                                                                    >  4:75965/1‑142 (MQ=38)
                                       cccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATGTTCCCGGTCAC‑GTCGAacgc                                                                                                  <  2:251244/149‑1 (MQ=38)
                                           cGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATGTTCCCGGTCAC‑GTCGAACGcccag                                                                                              <  3:75965/149‑4 (MQ=37)
                                            ggCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGc                                                                                             <  2:203370/149‑1 (MQ=255)
                                               cTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGa                                                                                          <  4:126767/149‑1 (MQ=255)
                                                  tgGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGAt                                                                                         <  5:280056/147‑1 (MQ=255)
                                                    gACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATg                                                                                     <  5:332851/149‑1 (MQ=255)
                                                    gACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATg                                                                                     >  7:34612/1‑149 (MQ=255)
                                                         aggagCGACGTGGCAATATCGAGGCCGTGGACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGtg                                                                                <  3:280017/149‑1 (MQ=255)
                                                          ggagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGtg                                                                                <  1:295843/148‑1 (MQ=255)
                                                           gatcgtcGTTGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGTGAACATGATCCCGGTCACGGTCGAGCGCGTGGCGGATCATGCCGtgct                                                                              <  5:192089/143‑1 (MQ=255)
                                                           gagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGtgct                                                                              >  2:253640/1‑149 (MQ=255)
                                                           gagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGtgct                                                                              <  6:48784/149‑1 (MQ=255)
                                                            agCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGtgctg                                                                             <  1:253640/149‑1 (MQ=255)
                                                            agCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGtg                                                                                >  3:81548/1‑146 (MQ=255)
                                                            agCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGtg                                                                                <  4:81548/146‑1 (MQ=255)
                                                                 cGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAAcc                                                                        <  1:132275/149‑1 (MQ=255)
                                                                      cAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAGCATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCACCCACTTa                                                                   >  1:153195/1‑149 (MQ=255)
                                                                           tCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCg                                                              >  6:2324/1‑149 (MQ=255)
                                                                            cgtggCCGTACACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGc                                                             >  8:136470/4‑149 (MQ=255)
                                                                                   gtgtACGGCCACGGCGCCGAGGATAACGCCTACAGAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGATCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCgg                                                       >  1:47095/1‑149 (MQ=255)
                                                                                       aCGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACtt                                                  >  4:236357/1‑149 (MQ=255)
                                                                                        cGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTg                                                 >  2:48356/1‑149 (MQ=255)
                                                                                                gcgcCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtgctg                                         <  7:136470/149‑1 (MQ=255)
                                                                                                   ccGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGctgctg                                      >  8:103797/1‑149 (MQ=255)
                                                                                                   ccGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGctgctg                                      <  1:229105/149‑1 (MQ=255)
                                                                                                    cGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGctgctg                                      >  2:199421/1‑148 (MQ=255)
                                                                                                    cGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTgg                                     <  5:2324/149‑1 (MQ=255)
                                                                                                     gAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtgctgc                                        <  8:157070/145‑1 (MQ=255)
                                                                                                     gAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtgctgc                                        >  7:157070/1‑145 (MQ=255)
                                                                                                     gAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGa                                    <  7:100931/149‑1 (MQ=255)
                                                                                                      aGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGAc                                   <  1:48356/149‑1 (MQ=255)
                                                                                                        gATAACGCCTACA‑AGGCGGAAAAGATCAAGCTCGCCTGGGCGTTCGA‑CGTGACCGGGAACCTGATACCGGTCACGGTCGAGCGCCTTGCGGATCATGCCGTGCTGGAACCACTTACCCCGACTGCGGCACTTGCCGTGCTGCTGGACaa                                 >  7:231958/1‑149 (MQ=255)
                                                                                                          tAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACg                               <  3:236357/149‑1 (MQ=255)
                                                                                                                   cA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGt                      <  2:243540/149‑1 (MQ=255)
                                                                                                                   cA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGATCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTg                        <  2:47095/147‑1 (MQ=255)
                                                                                                                          ggAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGCGAACAAGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACACCGCCTGCGTCACTTGCCCTGATGCTGGACAACGGGCCTTGGTCACAgc                >  3:85172/1‑149 (MQ=255)
                                                                                                                           gAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATc                                                                                        >  7:175137/1‑76 (MQ=255)
                                                                                                                           gAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATc                                                                                        <  8:175137/76‑1 (MQ=255)
                                                                                                                            aaaaGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAgcgc              >  6:137111/1‑149 (MQ=255)
                                                                                                                            aaaaGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAgcgc              <  8:34612/149‑1 (MQ=255)
                                                                                                                             aaaGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGGACATGATACCGGTCACGGTGGAGCGCCTGGCGGGTCATGGCGTGCTGGAACCACTTACCCCGCCTGCGGCCCTTGCCGGGCGGCTGGCCAACGGGGCGTGGTCACAgcggc             >  8:308908/1‑147 (MQ=255)
                                                                                                                                        gcgcCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCaa  >  2:200796/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
CTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAA  >  NC_002947/3762780‑3763063

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: ‑