Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,525,557 T→A 10.6% 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,5570TA10.6% 166.3 / 11.0 66V163D (GTC→GAC) mnmAtRNA 2‑thiouridine(34) synthase MnmA
Reads supporting (aligned to +/- strand):  ref base T (35/24);  new base A (4/3);  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 = 4.66e-01

TTCAAGGCGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGCGAGCGCCGCTTCAGC  >  NC_002947/4525409‑4525705
                                                                                                                                                    |                                                                                                                                                    
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                                                                                                                                      <  5:186478/149‑4 (MQ=255)
                aCTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGGCCGCCGACGGCGTGCAGGaagta                                                                                                                                      <  2:165680/149‑4 (MQ=255)
                aCTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATc                                                                                                                                      <  1:327962/149‑1 (MQ=255)
                     gCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCg                                                                                                                                 <  4:105012/149‑1 (MQ=255)
                     gCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCg                                                                                                                                 >  3:210356/1‑149 (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                         >  5:189355/1‑149 (MQ=255)
                                                                                                                             ccAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCt                         >  2:366416/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)
                                                                                                                                                    |                                                                                                                                                    
TTCAAGGCGTTCCTCGACTACGCCCTGTCGCTGGGTGCCGACCTGATCGCCACCGGCCACTATGTGCGCCGCCGCGATACTGGCGCGCTCACCGAACTGCTCAAGGGCCTGGACCCGAACAAGGACCAGAGCTACTTCCTGCACGCCGTCGGCGGCCAGGAAATCGCCCGCACACTGTTCCCGGTCGGCGAACTGGAAAAGCCCGAAGTACGCGCCATCGCCGAAAAACACGGCCTGGCCACCGCCAAGAAGAAAGATTCCACCGGCATCTGCTTCATTGGCGAGCGCCGCTTCAGC  >  NC_002947/4525409‑4525705

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: