Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,849,811 C→T 6.4% R118C (CGC→TGC)  rlmD → 23S rRNA (uracil(1939)‑C(5))‑methyltransferase RlmD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,849,8110CT 6.4% 183.2 / 2.1 62R118C (CGC→TGC) rlmD23S rRNA (uracil(1939)‑C(5))‑methyltransferase RlmD
Reads supporting (aligned to +/- strand):  ref base C (34/24);  new base T (2/2);  total (36/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

GCCCGTGGCAAAGTGGTCGAGGCCCGCCTGGAGCGGTTGCTGCAGGCCGCGC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGTGCGCTGGGACGTCAAGGCGCGTCAGCTGGAGGTGGGCTT  >  NC_002947/1849664‑1849959
                                                                                                                                                    |                                                                                                                                                      
gCCCGTGGCAAAGTGGTCGAGGCCCGCCTGGAGCGGTTGCTGCAGGCCGCGC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAgcg                                                                                                                                                       <  3:40987/149‑1 (MQ=255)
 cccGTGGCAAAGTGGTCGAGGCCCGCCTGGAGCGGTTGCTGCAGGCCGCGC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAgcgc                                                                                                                                                      <  7:277645/149‑1 (MQ=255)
   cGTGGCAAAGTGGTCGAGGCCCGCCTGGAGCGGTTGCTGCAGGCCGCGC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGca                                                                                                                                                     >  3:318884/1‑148 (MQ=255)
    gTGGCAAAGTGGTCGAGGCCCGCCTGGAGCGGTTGCTGCAGGCCGCGC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGcaca                                                                                                                                                   >  4:282596/1‑149 (MQ=255)
        cAAAGTGGTCGAGGCCCGCCTGGAGCGGTTGCTGCAGGCCGCGC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCg                                                                                                                                               <  4:318884/149‑1 (MQ=255)
         aaaGTGGTCGAGGCCCGCCTGGAGCGGTTGCTGCAGGCCGCGC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGc                                                                                                                                              >  3:4588/1‑149 (MQ=255)
         aaaGTGGTCGAGGCCCGCCTGGAGCGGTTGCTGCAGGCCGCGC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGc                                                                                                                                              >  5:278854/1‑149 (MQ=255)
                  gAGGCCCGCCTGGAGCGGTTGCTGCAGGCCGCGCACCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGt                                                                                                                                     >  1:242005/1‑149 (MQ=255)
                  gAGGCCCGCCTGGAGCGGTTGCTGCAGGCCGCGC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGtt                                                                                                                                    >  7:164828/1‑149 (MQ=255)
                       ccGCCTGGAGCGGTTGCTGCAGGCCGCGC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGTGCACACTCGCCGA‑GCAGTTGCAGc                                                                                                                               >  8:42089/1‑149 (MQ=255)
                       ccGCCTGGAGCGGTTGCTGCAGGCCGCGC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGTGCACACTCGCCGA‑GCAGTTGCAGc                                                                                                                               >  8:36288/1‑149 (MQ=255)
                       ccGCCTGGAGCGGTTGCTGCAGGCCGCGC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGc                                                                                                                               >  6:106031/1‑149 (MQ=255)
                                  ggTTGCTGCAGGCCGCGC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGc                                                                                                                    <  2:242005/149‑1 (MQ=255)
                                   gTTGCTGCAGGGCGCGC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTTGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCg                                                                                                                   <  3:282596/149‑1 (MQ=255)
                                   gTTGCTGCAGGCCGCGC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGtt                                                                                                                                    >  5:85598/1‑132 (MQ=255)
                                   gTTGCTGCAGGCCGCGC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGtt                                                                                                                                    <  6:85598/132‑1 (MQ=255)
                                    ttGCTGCAGGCCGCGC‑CCGAGCGCCGTGAGGCGCCATGGCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGt                                                                                                                  <  5:197489/149‑1 (MQ=255)
                                       ctgcAGGCCGCGC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGc                                                                                                                >  6:22800/1‑148 (MQ=255)
                                        tgcAGGCCGCGC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTg                                                                                                                          >  3:353002/1‑137 (MQ=255)
                                        tgcAGGCCGCGC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTg                                                                                                                          <  4:353002/137‑1 (MQ=255)
                                           aGGCCGCGC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGc                                                                                                                        <  8:118558/136‑1 (MQ=255)
                                           aGGCCGCGC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGc                                                                                                                        >  7:118558/1‑136 (MQ=255)
                                               cgcgC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAgg                                                                                                       >  8:139002/1‑149 (MQ=255)
                                                  gc‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGt                                                                                                    <  8:164828/149‑1 (MQ=255)
                                                         gcgcCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTg                                                                                              <  7:298035/149‑1 (MQ=255)
                                                                  ggCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGTCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCt                                                                                       >  7:188554/1‑147 (MQ=255)
                                                                    cgcCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGc                                                                                   <  6:278854/149‑1 (MQ=255)
                                                                     gcCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGTCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCg                                                                                  <  8:188554/149‑1 (MQ=255)
                                                                      ccATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCgg                                                                                 >  6:260831/1‑149 (MQ=255)
                                                                        aTGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGc                                                                               >  8:3053/1‑149 (MQ=255)
                                                                        aTGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGc                                                                               <  7:95715/149‑1 (MQ=255)
                                                                           tCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCg                                                                             <  2:347502/148‑1 (MQ=255)
                                                                                 ctacGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTc                                                                       <  5:260831/148‑1 (MQ=255)
                                                                                 ctacGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCg                                                                      >  2:347632/1‑149 (MQ=255)
                                                                                  tacGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCgg                                                                     >  7:86243/1‑149 (MQ=255)
                                                                                       ccGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTAcc                                                                >  2:103635/1‑149 (MQ=255)
                                                                                         gCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGTGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGTAATTCGGCTACCgg                                                              <  7:42089/149‑1 (MQ=255)
                                                                                         gCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGTGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCgg                                                              <  7:36288/149‑1 (MQ=255)
                                                                                                   cTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGccc                                                   >  3:44519/1‑149 (MQ=255)
                                                                                                    tGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCg                                                  >  2:41666/1‑149 (MQ=255)
                                                                                                    tGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCg                                                  >  5:87337/1‑149 (MQ=255)
                                                                                                     gCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGAGGCAGTTTAACCGCGTGTCCCCCGTGAAACCCCAGCCGAGGGCTCACCCCCAGACCGGGCCCGAAATAGGCTAACCCCGC‑CGGGcccc                                                  >  1:108994/1‑148 (MQ=255)
                                                                                                       aaCCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCgtgt                                               >  4:325172/1‑149 (MQ=255)
                                                                                                          cTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGc                                            >  3:154168/1‑149 (MQ=255)
                                                                                                           tGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGcc                                           >  5:167228/1‑149 (MQ=255)
                                                                                                             cagcaCCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGt                                         <  5:22800/149‑1 (MQ=255)
                                                                                                             cagcaCCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGt                                         <  4:4588/149‑1 (MQ=255)
                                                                                                                 aCCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGTgcgc                                     >  6:274336/1‑149 (MQ=255)
                                                                                                                 aCCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGTgcgc                                     >  7:132238/1‑149 (MQ=255)
                                                                                                                      ccGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGTGCGCTggg                                 <  4:245441/148‑1 (MQ=255)
                                                                                                                      ccGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGTGCGCTggg                                 >  3:245441/1‑148 (MQ=255)
                                                                                                                      ccGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGTGCGCTGGGa                                <  7:139002/149‑1 (MQ=255)
                                                                                                                        gCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGTGCGCTGGGACg                              >  4:310677/1‑149 (MQ=255)
                                                                                                                         cACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCGGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGCACGGGCCCGTGTCGCCGTGCGCTGGGACg                              <  1:214329/149‑1 (MQ=255)
                                                                                                                               gcgcAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGTGCGCTGGGACGTCAAGgc                       <  7:3053/149‑1 (MQ=255)
                                                                                                                                 gcAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGTGCGCTGGGACGTCAAGgcgc                     <  2:372939/149‑1 (MQ=255)
                                                                                                                                       ggCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGTGCGCTGGGACGTCAAGGCGCGTCAGc               >  4:190969/1‑149 (MQ=255)
                                                                                                                                       ggCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGTGCGCTGGGACGTCAACGCGCGTCAGc               >  7:158949/1‑149 (MQ=255)
                                                                                                                                            tGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGTGCGCTGGGACGTCAAGGCGCGTCAGCTGGAg          >  4:283103/1‑149 (MQ=255)
                                                                                                                                            tGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGTGCGCTGGGACGTCAAGGCGCGTCAGCTGGAg          >  4:126478/1‑149 (MQ=255)
                                                                                                                                              aaGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGTGCGCTGGGACGTCAAGGCGCGTCAGCTGGAg          <  4:44519/147‑1 (MQ=255)
                                                                                                                                              aaGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGTGCGCTGGGACGTCAAGGCGCGTCAGCTGGAGGt        <  6:167228/149‑1 (MQ=255)
                                                                                                                                              aaGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGTGCGCTGGGACGTCAAGGCGCGTCAGCTGGAGGt        <  5:106031/149‑1 (MQ=255)
                                                                                                                                              aaGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGTGCGCTGGGACGTCAAGGCGCGTCAGCTGGAGGt        >  4:341633/1‑149 (MQ=255)
                                                                                                                                               agcagcGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGTGCGCTGGGACGTCAAGGCGCGTCAGCTGGAGGTg       <  5:274336/149‑1 (MQ=255)
                                                                                                                                                gcagcGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGTGCGCTGGGACGTCAAGGCGCGTCAGCTGGAGGTgg      <  3:190969/149‑1 (MQ=255)
                                                                                                                                                gcagcGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGTGCGCTGGGACGTCAAGGCGCGTCAGCTGGAGGTgg      <  1:347632/149‑1 (MQ=255)
                                                                                                                                                    cgcACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGAGCCCGTGTCGCCGTGCGCTGGGACGTCAAGGCGCGTCAGCTGGAGGTGGGCtt  >  3:148266/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GCCCGTGGCAAAGTGGTCGAGGCCCGCCTGGAGCGGTTGCTGCAGGCCGCGC‑CCGAGCGCCGTGAGGCGCCATGTCGCTACTACGACCGCTGCGGTGGCTGCAACCTGCAGCACCTGCCGCACGAAGCGCAACTGGCGCTGAAGCAGCGCACACTCGCCGA‑GCAGTTGCAGCGGGTGGCCGGCGTGCAACCCGAGGCGTGGGCTGCCCCCCTGAGCGGGCCGGAATTCGGCTACCGGCGC‑CGGGCCCGTGTCGCCGTGCGCTGGGACGTCAAGGCGCGTCAGCTGGAGGTGGGCTT  >  NC_002947/1849664‑1849959

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