Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,525,550 C→G 7.7% H161D (CAC→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,5500CG 7.7% 204.4 / 6.9 65H161D (CAC→GAC) ‡mnmAtRNA 2‑thiouridine(34) synthase MnmA
Reads supporting (aligned to +/- strand):  ref base C (37/23);  new base G (2/3);  total (39/26)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

TCAAGTTCAAGGCGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGC  >  NC_002947/4525404‑4525690
                                                                                                                                                  |                                                                                                                                            
tCAAGTTCAAGGCGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGAACCGAACAAGGACCAGAGCTACTTCCTGCAc                                                                                                                                            >  2:327962/1‑149 (MQ=255)
  aaGTTCAAGGCGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGc                                                                                                                                          >  5:56925/1‑149 (MQ=255)
     ttCAAGGCGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGt                                                                                                                                       >  4:53631/1‑149 (MQ=255)
            cGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCggcggc                                                                                                                                <  3:177626/149‑1 (MQ=255)
             gTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGGCCGCCGAcggcgtgc                                                                                                                               <  3:28625/149‑4 (MQ=255)
                 cTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAgg                                                                                                                            >  1:8711/1‑148 (MQ=255)
                 cTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGa                                                                                                                           <  5:203912/149‑1 (MQ=255)
                     aCTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGGCCGCCGACGGCGTGCAGGaagta                                                                                                                       <  2:165680/149‑4 (MQ=255)
                     aCTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGGCCGCCGACGGCGTGCAGGaagta                                                                                                                       <  5:186478/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                                                                                                 >  8:356437/1‑148 (MQ=255)
                                            gACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCg                                                                                                 <  7:356437/148‑1 (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)
                                                                                                                                                  |                                                                                                                                            
TCAAGTTCAAGGCGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGC  >  NC_002947/4525404‑4525690

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