Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,525,557 T→A 7.5% V163D (GTC→GAC)  mnmA → tRNA 2‑thiouridine(34) synthase MnmA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,525,5570TA 7.5% 188.5 / 5.1 67V163D (GTC→GAC) mnmAtRNA 2‑thiouridine(34) synthase MnmA
Reads supporting (aligned to +/- strand):  ref base T (31/31);  new base A (2/3);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

GCGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGCGAGCGCC  >  NC_002947/4525415‑4525697
                                                                                                                                               |                                                                                                                                               
gCGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTcggcgg                                                                                                                                         <  7:148177/149‑1 (MQ=255)
 cGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCggcggc                                                                                                                                        <  2:133179/149‑1 (MQ=255)
     ccTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAg                                                                                                                                     >  2:297831/1‑148 (MQ=255)
      cTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAgg                                                                                                                                   <  5:328387/148‑1 (MQ=255)
      cTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAgg                                                                                                                                   >  2:229775/1‑148 (MQ=255)
      cTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAgg                                                                                                                                   >  5:138851/1‑148 (MQ=255)
        cGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGaaa                                                                                                                                >  2:322920/1‑149 (MQ=255)
         gACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGaaa                                                                                                                                <  6:328960/148‑1 (MQ=255)
         gACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGaa                                                                                                                                 <  6:138851/147‑1 (MQ=255)
         gACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGa                                                                                                                                  >  6:366692/1‑146 (MQ=255)
         gACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGa                                                                                                                                  <  5:366692/146‑1 (MQ=255)
           cTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGGCCGCCGACGGCGTGCAGGaagtag                                                                                                                             <  5:215754/149‑5 (MQ=255)
           cTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCg                                                                                                                             >  1:279785/1‑149 (MQ=255)
             aCGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGGCCGCCGAcggcgtgca                                                                                                                                      >  6:223371/1‑135 (MQ=255)
             aCGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGGCCGCCGAcggcgtgca                                                                                                                                      <  5:223371/139‑5 (MQ=255)
              cGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGccc                                                                                                                          >  2:337519/1‑149 (MQ=255)
              cGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGccc                                                                                                                          <  5:230608/149‑1 (MQ=255)
                cccTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCg                                                                                                                         >  4:136298/1‑148 (MQ=255)
                       gCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCAGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGt                                                                                                                 >  7:193578/1‑149 (MQ=255)
                                  atcaggtCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGc                                                                                                      >  2:87912/7‑149 (MQ=255)
                                  aCCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGc                                                                                                      >  3:44283/1‑149 (MQ=255)
                                    cTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGa                                                                                                    >  1:127573/1‑149 (MQ=255)
                                         cGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTg                                                                                                >  6:286966/1‑148 (MQ=255)
                                           ccACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGaa                                                                                             <  2:127573/149‑1 (MQ=255)
                                            cACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGaa                                                                                             >  6:352718/1‑148 (MQ=255)
                                                 gCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGccc                                                                                       >  2:351481/1‑149 (MQ=255)
                                                 gCCACTATGCGCGCCGGCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTACCCGGTCGGCGAACTGGAAAAGccc                                                                                       <  4:256215/149‑1 (MQ=255)
                                                   cACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCg                                                                                      >  4:31332/1‑148 (MQ=255)
                                                    aCTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGaa                                                                                    <  1:229775/149‑1 (MQ=255)
                                                        tgtgCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGTCCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGCGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGt                                                                                  <  1:271967/148‑1 (MQ=255)
                                                            cgccgccgcGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTAcgcg                                                                             <  2:279785/148‑1 (MQ=255)
                                                                gccgcGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATc                                                                        <  3:136298/149‑1 (MQ=255)
                                                                 ccgcTATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCg                                                                       <  4:44283/149‑1 (MQ=255)
                                                                 ccgcGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCg                                                                       >  7:107211/1‑149 (MQ=255)
                                                                  cgcgATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGc                                                                      <  5:309675/149‑1 (MQ=255)
                                                                   gcgATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGcc                                                                     >  3:245301/1‑149 (MQ=255)
                                                                         cTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGaaaaa                                                               <  1:351481/149‑1 (MQ=255)
                                                                            gcgcgcTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAacac                                                            <  5:352718/149‑1 (MQ=255)
                                                                            gcgcgcTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAacac                                                            >  1:64086/1‑149 (MQ=255)
                                                                                 cTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCt                                                       <  1:337519/149‑1 (MQ=255)
                                                                                    aCCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGc                                                    >  7:274507/1‑149 (MQ=255)
                                                                                     ccGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGaaaa                                                                >  5:235329/1‑136 (MQ=255)
                                                                                     ccGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGaaaa                                                                <  6:235329/136‑1 (MQ=255)
                                                                                              tgaacgtccTGTAGCCGAACAAGTACGACAGCTACTTCTTGCAGGCCGCTACTAGCTGCCAGTAAATCGCCCGCACACTGTTCCCGGTCGTCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGcc                                             <  1:297831/142‑1 (MQ=255)
                                                                                                  gggCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACTGCCTGGCCACCGCCAagaaga                                      <  8:193578/149‑1 (MQ=255)
                                                                                                  gggCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAagaaga                                      <  1:322920/149‑1 (MQ=255)
                                                                                                       tGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGAt                                 >  1:157616/1‑149 (MQ=255)
                                                                                                         gACGCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCTGCGAGGAAATCGCCCGCACACTGTTCCCGGGCGGCGAACTGGAAAAGCCCGAAGTACGCGCAATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGAt                                 <  7:172494/147‑1 (MQ=255)
                                                                                                         gACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTc                               <  7:50854/149‑1 (MQ=255)
                                                                                                         gACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTc                               <  3:31332/149‑1 (MQ=255)
                                                                                                            ccGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCAc                            >  7:70248/1‑149 (MQ=255)
                                                                                                             cGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCTGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTTGCCACCGCCAAGAAGAAAGATTCCAcc                           <  8:274507/149‑1 (MQ=255)
                                                                                                              gAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGaagaaaga                                  <  5:9257/141‑1 (MQ=255)
                                                                                                              gAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGaagaaaga                                  >  6:9257/1‑141 (MQ=255)
                                                                                                               aaCAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCaagaag                                       >  8:191043/1‑135 (MQ=255)
                                                                                                               aaCAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCaagaag                                       <  7:191043/135‑1 (MQ=255)
                                                                                                               aaCAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCgg                         >  7:203130/1‑149 (MQ=255)
                                                                                                                    ggACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCt                    <  8:70248/149‑1 (MQ=255)
                                                                                                                      aCCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGc                  <  4:245301/149‑1 (MQ=255)
                                                                                                                       ccAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCt                 >  5:293776/1‑149 (MQ=255)
                                                                                                                          gagCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCa              >  7:285878/1‑149 (MQ=255)
                                                                                                                                tttCCTGGCCGCCGACGGCGTGCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGc         >  6:15988/2‑148 (MQ=255)
                                                                                                                                tttCCTGGCCGCCGACGGCGTGCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGc         <  5:15988/147‑1 (MQ=255)
                                                                                                                                  tCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGCGAg      >  6:248016/1‑149 (MQ=255)
                                                                                                                                     tGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGaagaaa                                    <  1:105371/116‑1 (MQ=255)
                                                                                                                                     tGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGaagaaa                                    >  2:105371/1‑116 (MQ=255)
                                                                                                                                      gCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGCGAGcgcc  >  2:63082/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                               
GCGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGCGAGCGCC  >  NC_002947/4525415‑4525697

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: