Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,641,627 T→C 9.0% N22D (AAC→GAC) ‡ 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,6270TC 9.0% 209.3 / 9.6 78N22D (AAC→GAC) ‡cmoAtRNA 5‑methoxyuridine(34)/uridine‑5‑oxyacetic acid(34) metyltransferase CmoA
Reads supporting (aligned to +/- strand):  ref base T (30/41);  new base C (3/4);  total (33/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

ACGGCACCCAGCGATGCGCCCAGATCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGCCAAACCGCAGCCGGGATGCATGGTTGGCTGACCAGCGCTTT  >  NC_002947/1641479‑1641775
                                                                                                                                                    |                                                                                                                                                    
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                                                                                                                                                  <  8:240571/149‑1 (MQ=255)
    cACCCAGCGATGCGCCCAGATCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAg                                                                                                                                                  <  5:32930/149‑1 (MQ=255)
       ccAGCGATGCGCCCAGATCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGAc                                                                                                                                                >  8:303484/1‑148 (MQ=255)
       ccAGCGATGCGCCCAGATCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGAc                                                                                                                                                >  2:170661/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                                                                                                                    <  8:28468/149‑1 (MQ=255)
                                  cGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTgg                                                                                                                    <  1:170661/149‑1 (MQ=255)
                                  cGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTgg                                                                                                                    <  7:34316/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                                                                                   >  4:428350/1‑148 (MQ=255)
                                                                    cccGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCg                                                                                   >  6:432334/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:378396/149‑1 (MQ=255)
                                                                                gACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCt                                                                      >  3:87713/1‑149 (MQ=255)
                                                                                gACAATGGTTGGATAACCGGGCACCGAGCGCTTGATAATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTAGGGCACCTCCTCCAGTGGTTGGGAGAAAAGGCCGTCGGGTAGTTTGCTAACGGCGAATCCACGCGGCt                                                                      >  7:196863/1‑149 (MQ=255)
                                                                                  cAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTgg                                                                    <  7:189086/149‑1 (MQ=255)
                                                                                            aTAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAg                                                          <  6:307656/149‑1 (MQ=255)
                                                                                               aCCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAgggg                                                       <  5:101782/149‑1 (MQ=255)
                                                                                               aCCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAgggg                                                       <  4:87713/149‑1 (MQ=255)
                                                                                               aCCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAgggg                                                       <  3:428350/149‑1 (MQ=255)
                                                                                               aCCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAggg                                                        <  6:216830/148‑1 (MQ=255)
                                                                                                ccGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATa                                                               >  2:257313/1‑140 (MQ=255)
                                                                                                ccGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATa                                                               <  1:257313/140‑1 (MQ=255)
                                                                                                ccGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAggggg                                                      >  2:131625/1‑149 (MQ=255)
                                                                                                 cGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGgc                                                     <  1:413006/149‑1 (MQ=255)
                                                                                                  gggCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAgg                                                                           >  5:183696/1‑126 (MQ=255)
                                                                                                  gggCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAgg                                                                           <  6:183696/126‑1 (MQ=255)
                                                                                                           gcgcTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGc                                           <  5:432334/149‑1 (MQ=255)
                                                                                                            cgcTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCAtt                                                <  7:99946/143‑1 (MQ=255)
                                                                                                            cgcTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCAtt                                                >  8:99946/1‑143 (MQ=255)
                                                                                                               ttGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGCCaaa                                       >  3:169934/1‑149 (MQ=255)
                                                                                                                 gATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGCCaacca                                     >  1:74865/1‑146 (MQ=255)
                                                                                                                   tCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCAAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGCCAAACCGCa                                  >  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                                           <  6:12941/133‑1 (MQ=255)
                                                                                                                           gggAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGc                                           >  5:12941/1‑133 (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)
                                                                                                                                                    tGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGCCAAACCGCAGCCGGGATGCATGGTTGGCTGACCAGCGCttt  >  1:149409/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACGGCACCCAGCGATGCGCCCAGATCGTACAGTGCGGTGCCTGGCTGGGCGAACCGAGCGGCCAGTACCCCGAGGTTCTCGACAATGGTTGGATAACCGGGCACCGAGCGCTTGATCATGTCCGGGAACACACGCACGACGTCTTCGTTGAAGACGAAGTCGGGCACCTGCTCCAGTGGTTGGGAGAAAAGGCGGTCGGGTTGTTTGCTCACGGCGAATCCAGGCGGCTGGTGATAGGCAGGGGGCGCATTCTAGCCAAACCGCAGCCGGGATGCATGGTTGGCTGACCAGCGCTTT  >  NC_002947/1641479‑1641775

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: