Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,525,564 C→G 9.1% G165G (GGC→GGG) ‡ 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,5640CG 9.1% 197.4 / 8.8 66G165G (GGC→GGG) ‡mnmAtRNA 2‑thiouridine(34) synthase MnmA
Reads supporting (aligned to +/- strand):  ref base C (34/25);  new base G (4/2);  total (39/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

CGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGCGAGCGCCGCTTCAGCGACTTCC  >  NC_002947/4525416‑4525712
                                                                                                                                                    |                                                                                                                                                    
cGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCggcggc                                                                                                                                                      <  3:177626/149‑1 (MQ=255)
     cTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAgg                                                                                                                                                  >  1:8711/1‑148 (MQ=255)
     cTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGa                                                                                                                                                 <  5:203912/149‑1 (MQ=255)
         aCTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGGCCGCCGACGGCGTGCAGGaagta                                                                                                                                             <  5:186478/149‑4 (MQ=255)
         aCTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGGCCGCCGACGGCGTGCAGGaagta                                                                                                                                             <  2:165680/149‑4 (MQ=255)
         aCTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATc                                                                                                                                             <  1:327962/149‑1 (MQ=255)
              gCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCg                                                                                                                                        >  3:210356/1‑149 (MQ=255)
              gCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCg                                                                                                                                        <  4:105012/149‑1 (MQ=255)
                 cTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGca                                                                                                                                      >  5:242491/1‑148 (MQ=255)
                      gCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGt                                                                                                                                <  6:56925/149‑1 (MQ=255)
                            tGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCgg                                                                                                                          <  2:167758/149‑1 (MQ=255)
                                gACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCgg                                                                                                                      <  3:53631/149‑1 (MQ=255)
                                gACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCg                                                                                                                       <  7:356437/148‑1 (MQ=255)
                                gACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCg                                                                                                                       >  8:356437/1‑148 (MQ=255)
                                       tCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTg                                                                                                               >  1:375738/1‑149 (MQ=255)
                                         gCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGGTCACCGAACTGCTCAAGGGCCTGGACCCGCACAAGGACCAGAGCTACTTCCTGCACGCCGTCTGCGGACAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTgg                                                                                                              >  8:251077/1‑148 (MQ=255)
                                          ccACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGa                                                                                                             >  2:164411/1‑148 (MQ=255)
                                           cACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGaa                                                                                                            >  7:133922/1‑148 (MQ=255)
                                            aCCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGaaaa                                                                                                          <  2:8711/149‑1 (MQ=255)
                                             ccGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAg                                                                                                         >  5:189929/1‑149 (MQ=255)
                                               ggCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGcc                                                                                                       <  1:184566/149‑1 (MQ=255)
                                                gCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGccc                                                                                                      >  3:1668/1‑149 (MQ=255)
                                                       tgtgCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCTCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTAc                                                                                               >  5:257956/1‑149 (MQ=255)
                                                          gcgcCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGCACTGGAAAAGCCCGAAGTCcgcg                                                                                            >  2:153203/1‑149 (MQ=255)
                                                          gcgcCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTAcgcg                                                                                            >  3:157513/1‑149 (MQ=255)
                                                               gccgcGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCCTc                                                                                       >  4:203307/1‑149 (MQ=255)
                                                                        cTGGCGCGCTCACCGAACTGCTCAAGGGCCTTGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGaaaaa                                                                              <  3:67093/149‑1 (MQ=255)
                                                                          ggCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAaca                                                                            >  3:340560/1‑149 (MQ=255)
                                                                           gcgcgcTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAacac                                                                           <  4:112219/149‑1 (MQ=255)
                                                                            cgcgcTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACg                                                                          <  6:189929/149‑1 (MQ=255)
                                                                             gcgcTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAacacag                                                                         <  5:325207/149‑3 (MQ=255)
                                                                                   aCCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGc                                                                   >  1:286290/1‑149 (MQ=255)
                                                                                     cGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCa                                                                 >  5:139859/1‑149 (MQ=255)
                                                                                            cTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAGGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGcc                                                            >  5:324726/1‑147 (MQ=255)
                                                                                            cTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGcc                                                            >  4:302623/1‑147 (MQ=255)
                                                                                             tCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAACCACGGCATGGCCCCCGCCat                                                          >  8:88025/1‑147 (MQ=255)
                                                                                             tCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACg                                                                          >  8:305663/1‑132 (MQ=255)
                                                                                             tCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACg                                                                          <  7:305663/132‑1 (MQ=255)
                                                                                             tCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGcc                                                            >  1:152743/1‑146 (MQ=255)
                                                                                             tCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGcc                                                            >  8:73081/1‑146 (MQ=255)
                                                                                               aaGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCaag                                                         <  3:289231/147‑1 (MQ=255)
                                                                                                aGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCaag                                                         >  2:318223/1‑146 (MQ=255)
                                                                                                  ggccgcgACCCGAACAACGACCAGAGCTACTTCCAGCACGCCGTCCGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAgaagaa                                                    <  4:213232/143‑1 (MQ=255)
                                                                                                    ccTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAagaaga                                                     >  1:234050/1‑146 (MQ=255)
                                                                                                     cTGGACCCGAACAAGGACCAGAGCCACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCCGCGAACTGCAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAagaaga                                                     >  7:80415/1‑145 (MQ=255)
                                                                                                         aCCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTc                                              >  1:79731/1‑148 (MQ=255)
                                                                                                          cccGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCa                                            <  4:1668/149‑1 (MQ=255)
                                                                                                               aCAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGc                                       <  1:164411/149‑1 (MQ=255)
                                                                                                                 aaGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCAt                                     >  2:340323/1‑149 (MQ=255)
                                                                                                                   ggACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCt                                   >  2:207902/1‑149 (MQ=255)
                                                                                                                   ggACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCt                                   >  8:239274/1‑149 (MQ=255)
                                                                                                                     aCCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGc                                 <  6:242491/149‑1 (MQ=255)
                                                                                                                      ccAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCTCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCt                                <  6:257956/149‑1 (MQ=255)
                                                                                                                      ccAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCt                                >  2:366416/1‑149 (MQ=255)
                                                                                                                      ccAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCt                                >  5:189355/1‑149 (MQ=255)
                                                                                                                       cAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACCCTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCtt                               >  4:330714/1‑149 (MQ=255)
                                                                                                                         gagCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCa                             >  3:149695/1‑149 (MQ=255)
                                                                                                                           gCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCAtt                           <  7:73081/149‑1 (MQ=255)
                                                                                                                            cTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTg                          <  7:239274/149‑1 (MQ=255)
                                                                                                                               tttCCTGGCCGCCGACGGCGTGCAGGAAATTGCCAGCACAATGTTCCCGGTCGGCGCACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGCCCTGTCCACCGCCAAgaagaac                                                   >  2:17770/2‑120 (MQ=255)
                                                                                                                               tttCCTGGCCGCCGACGGCGTGCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGc                        >  3:365916/2‑148 (MQ=255)
                                                                                                                                   ctggccGCCGACGGCGTGCAGGAAATCGCCCGCACATTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATATGCTTCATTGGCGAgcg                   >  3:204983/6‑149 (MQ=255)
                                                                                                                                   ctggccGCCGACGGCGTGCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGCGAgcg                   >  4:76841/6‑149 (MQ=255)
                                                                                                                                           cGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGCGAGCGCCGCtt             >  3:45685/1‑147 (MQ=255)
                                                                                                                                           cGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGCGAGCGCCGCTTCa           <  1:369551/149‑1 (MQ=255)
                                                                                                                                             tCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGCGAGCGCCGCTTCAGc         <  3:302623/149‑1 (MQ=255)
                                                                                                                                              cggcggCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGCGAGCGCCGCTTCAGCg        <  2:375738/149‑1 (MQ=255)
                                                                                                                                                gcggcCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGCGAGCGCCGCTTCAGCGAc      <  6:139859/149‑1 (MQ=255)
                                                                                                                                                    ccAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGCGAGCGCCGCTTCAGCGACTTcc  >  8:340690/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGCGAGCGCCGCTTCAGCGACTTCC  >  NC_002947/4525416‑4525712

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: