Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,641,625 G→A 8.8% N22N (AAC→AAT) ‡ cmoA ← tRNA 5‑methoxyuridine(34)/uridine‑5‑oxyacetic acid(34) metyltransferase CmoA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,641,6250GA 8.8% 225.1 / 8.2 79N22N (AAC→AAT) ‡cmoAtRNA 5‑methoxyuridine(34)/uridine‑5‑oxyacetic acid(34) metyltransferase CmoA
Reads supporting (aligned to +/- strand):  ref base G (30/42);  new base A (3/4);  total (33/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-01

ACGGCACCCAGCGATGCGCCCAGATCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGCCAAACCGCAGCCGGGATGCATGGTTGGCTGACCAGCGCT  >  NC_002947/1641479‑1641773
                                                                                                                                                  |                                                                                                                                                    
aCGGCACCCAGCGATGCGCCCAGATCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGtt                                                                                                                                                    <  3:194708/149‑1 (MQ=255)
 cGGCACCCAGCGATGCGCCCAGATCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGtt                                                                                                                                                    >  5:131788/1‑148 (MQ=255)
 cGGCACCCAGCGATGCGCCCAGATCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTg                                                                                                                                                   <  6:131788/149‑1 (MQ=255)
  ggCACCCAGCGATGCGCCCAGATCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGa                                                                                                                                                  <  5:62815/149‑1 (MQ=255)
   gCACCCAGCGATGCGCCCAGATCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGaa                                                                                                                                                 <  3:4020/149‑1 (MQ=255)
   gCACCCAGCGATGCGCCCAGATCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGaa                                                                                                                                                 <  2:260217/149‑1 (MQ=255)
    cACCCAGCGATGCGCCCAGATCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTg                                                                                                                                                   >  2:123320/1‑146 (MQ=255)
    cACCCAGCGATGCGCCCAGATCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAg                                                                                                                                                <  5:32930/149‑1 (MQ=255)
    cACCCAGCGATGCGCCCAGATCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAg                                                                                                                                                <  8:240571/149‑1 (MQ=255)
       ccAGCGATGCGCCCAGATCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGAc                                                                                                                                              >  2:170661/1‑148 (MQ=255)
       ccAGCGATGCGCCCAGATCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGAc                                                                                                                                              >  8:303484/1‑148 (MQ=255)
         aGCGATGCGCCCAGATCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACg                                                                                                                                             >  6:101782/1‑147 (MQ=255)
         aGCGATGCGCCCAGATCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCAACGAAGACgtc                                                                                                                                           <  2:39144/149‑3 (MQ=255)
           cGATGCGCCCAGATCGGACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCAACGAAGACgtcgt                                                                                                                                         <  4:315041/149‑5 (MQ=255)
             aTGCGCCCAGATCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCAACGAAGACgtcgtgc                                                                                                                                       <  3:132190/149‑7 (MQ=255)
              gtcgcCCAGATCTTACAGTGCGTTGCCTTGATGGGCGAACCGCGCGTCCAGTACCCCGAGGTTCTCGACAATGTTTTGATAACCGTGCACCTATAGATTGATCATGTCCGGGAACACAAGCACTACGTCTTCGTTTAAGACGAAGTCg                                                                                                                                       <  4:416088/146‑1 (MQ=255)
                 gcCCAGATCGTACAGAGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACAGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGc                                                                                                                                    >  8:189086/1‑148 (MQ=255)
                    cAGCTCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCt                                                                                                                                <  6:55611/149‑1 (MQ=255)
                        tCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTc                                                                                                                            <  8:312790/149‑1 (MQ=255)
                             cAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTg                                                                                                                       <  4:404435/149‑1 (MQ=255)
                                  cGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTgg                                                                                                                  <  1:170661/149‑1 (MQ=255)
                                  cGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTgg                                                                                                                  <  7:34316/149‑1 (MQ=255)
                                  cGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTgg                                                                                                                  <  8:28468/149‑1 (MQ=255)
                                       ccTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGgaga                                                                                                              >  4:378396/1‑148 (MQ=255)
                                        ctggctggGCGAACCGAGCGGCCTGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTGGGGCACCTGCTCCAGTGGTTGGGAGaaa                                                                                                            <  5:193718/149‑1 (MQ=255)
                                           gctggGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAgg                                                                                                         <  2:364131/149‑1 (MQ=255)
                                           gctggGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCATCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAgg                                                                                                         <  3:69124/149‑1 (MQ=255)
                                            ctggGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGc                                                                                                        <  2:373676/149‑1 (MQ=255)
                                              gggCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCgg                                                                                                      <  8:139234/149‑1 (MQ=255)
                                               ggCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGt                                                                                                     >  3:388855/1‑149 (MQ=255)
                                                gCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCAACGAAGACgt                                                                                                                                            >  6:261224/1‑108 (MQ=255)
                                                gCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCAACGAAGACgt                                                                                                                                            <  5:261224/109‑2 (MQ=255)
                                                 cGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCg                                                                                                   <  1:123320/149‑1 (MQ=255)
                                                     ccGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGgtt                                                                                               >  4:352181/1‑149 (MQ=255)
                                                       gAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTGCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGgttgt                                                                                             <  6:313748/149‑1 (MQ=255)
                                                        aGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGgttgtt                                                                                            >  3:186118/1‑149 (MQ=255)
                                                          cGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTg                                                                                          <  1:229637/149‑1 (MQ=255)
                                                             ccAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCt                                                                                        >  2:413006/1‑148 (MQ=255)
                                                                 tACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGGGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACgg                                                                                   >  7:70649/1‑149 (MQ=255)
                                                                  aCCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGc                                                                                  >  5:307656/1‑149 (MQ=255)
                                                                    cccGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCg                                                                                 >  6:432334/1‑148 (MQ=255)
                                                                    cccGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCg                                                                                 >  4:428350/1‑148 (MQ=255)
                                                                       gAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATc                                                                             >  1:400891/1‑149 (MQ=255)
                                                                         ggTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCa                                                                           <  2:147698/149‑1 (MQ=255)
                                                                         ggTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCa                                                                           <  7:30470/149‑1 (MQ=255)
                                                                           ttCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAgg                                                                         <  7:303484/149‑1 (MQ=255)
                                                                                gACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCt                                                                    >  3:87713/1‑149 (MQ=255)
                                                                                gACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCt                                                                    <  3:378396/149‑1 (MQ=255)
                                                                                gACAATGGTTGGATAACCGGGCACCGAGCGCTTGATAATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTAGGGCACCTCCTCCAGTGGTTGGGAGAAAAGGCCGTCGGGTAGTTTGCTAACGGCGAATCCACGCGGCt                                                                    >  7:196863/1‑149 (MQ=255)
                                                                                  cAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTgg                                                                  <  7:189086/149‑1 (MQ=255)
                                                                                            aTAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAg                                                        <  6:307656/149‑1 (MQ=255)
                                                                                               aCCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAgggg                                                     <  4:87713/149‑1 (MQ=255)
                                                                                               aCCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAgggg                                                     <  5:101782/149‑1 (MQ=255)
                                                                                               aCCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAgggg                                                     <  3:428350/149‑1 (MQ=255)
                                                                                               aCCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAggg                                                      <  6:216830/148‑1 (MQ=255)
                                                                                                ccGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATa                                                             <  1:257313/140‑1 (MQ=255)
                                                                                                ccGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATa                                                             >  2:257313/1‑140 (MQ=255)
                                                                                                ccGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAggggg                                                    >  2:131625/1‑149 (MQ=255)
                                                                                                 cGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGgc                                                   <  1:413006/149‑1 (MQ=255)
                                                                                                  gggCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAgg                                                                         <  6:183696/126‑1 (MQ=255)
                                                                                                  gggCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAgg                                                                         >  5:183696/1‑126 (MQ=255)
                                                                                                           gcgcTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGc                                         <  5:432334/149‑1 (MQ=255)
                                                                                                            cgcTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCAtt                                              >  8:99946/1‑143 (MQ=255)
                                                                                                            cgcTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCAtt                                              <  7:99946/143‑1 (MQ=255)
                                                                                                               ttGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGCCaaa                                     >  3:169934/1‑149 (MQ=255)
                                                                                                                 gATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGCCaacca                                   >  1:74865/1‑146 (MQ=255)
                                                                                                                   tCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAG‑AAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCAAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGCCAAACCGCa                                >  1:261975/1‑149 (MQ=255)
                                                                                                                     aTGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGCCAAACCGCa                                <  4:388855/148‑1 (MQ=255)
                                                                                                                         ccGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCGCCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGCCAAACCGCAGCCgg                           >  5:210299/1‑149 (MQ=255)
                                                                                                                           gggCACACACGCACGACGTCTTCGATGAAGTCGAAGTCGGGCACCTGCGCCAGTGGTGGGGAGAAAAGGCCGGCGGGGTGTTTGCTCCCGGCGAATCCAGGCGGCGTGGGATAGGCCGGGGGCGCATTCTAGCCAAACCGCAGCCGGGa                         <  8:196863/149‑1 (MQ=255)
                                                                                                                           gggAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGc                                         >  5:12941/1‑133 (MQ=255)
                                                                                                                           gggAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGc                                         <  6:12941/133‑1 (MQ=255)
                                                                                                                            ggAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGCCAAACCGCAGCCGGGAt                        <  2:400891/149‑1 (MQ=255)
                                                                                                                                 acacGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTg                                                                   >  3:216896/1‑101 (MQ=255)
                                                                                                                                 acacGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTg                                                                   <  4:216896/101‑1 (MQ=255)
                                                                                                                                   ccgcacgACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGGGATAGGCAGGGGGAGCATTCTAGCCAAACCGCAGCCGGGATGCATGGt                 >  6:93593/2‑149 (MQ=255)
                                                                                                                                       acttcgTCTTCAACGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGCCAAACCGCAGCCGGGATGCATGGTTGGc             >  8:227736/5‑149 (MQ=255)
                                                                                                                                        cttcgTCTTCAACGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGCCAAACCGCAGCCGGGATGCATGGTTGGCt            >  2:121470/4‑149 (MQ=255)
                                                                                                                                            gTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGCCAAACCGCAGCCGGGATGCATGGTTGGCTGAc         >  8:125246/1‑148 (MQ=255)
                                                                                                                                                  gTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGCCAAACCGCAGCCGGGATGCATGGTTGGCTGACCAGCGCt  >  4:384484/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ACGGCACCCAGCGATGCGCCCAGATCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGCCAAACCGCAGCCGGGATGCATGGTTGGCTGACCAGCGCT  >  NC_002947/1641479‑1641773

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