Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 854,655 T→C 7.3% V85A (GTC→GCC)  moeB → molybdopterin synthase adenylyltransferase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947854,6550TC 7.3% 193.8 / 4.9 68V85A (GTC→GCC) moeBmolybdopterin synthase adenylyltransferase subunit
Reads supporting (aligned to +/- strand):  ref base T (32/31);  new base C (2/3);  total (34/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 = 3.23e-01

GGGCTGGGTGGCCTGGGCTCGCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAATTTCTCTACCCGC  >  NC_002947/854507‑854800
                                                                                                                                                    |                                                                                                                                                   
gggCTGGGTGGCCTGGGCTCGCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGt                                                                                                                                                     <  1:331463/149‑1 (MQ=255)
gggCTGGGTGGCCTGGGCTCGCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGt                                                                                                                                                     <  5:246244/149‑1 (MQ=255)
 ggCTGGGTGGCGTGGGCTCGCCAGGCGGGCTGTACCTGGCCGCCGCTGGGGTGGGTGAGCTGCACCTGGCGGACTTCGACACCGTCGACCTGACCACCGTGCAACGCCAGGTCATCCACGTCAGTGCCAGCGTTGGCATGAGCAAGGTc                                                                                                                                                    <  6:226080/149‑1 (MQ=255)
       gTGGCCTGGGCTCGCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCg                                                                                                                                             >  2:349628/1‑149 (MQ=255)
            cTGGTCTGGCCAGTCGCCCTTTACCTGGCCGCCGCTGTCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACTCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGTTCGACTCGGCCtt                                                                                                                                        <  6:259620/149‑1 (MQ=255)
            cTGGGCTCGCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCtt                                                                                                                                        <  4:60295/149‑1 (MQ=255)
               ggCTCGCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCa                                                                                                                                     <  2:221915/149‑1 (MQ=255)
                    gCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCCCGACCGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCt                                                                                                                                >  4:206610/1‑149 (MQ=255)
                    gCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAgcg                                                                                                                                  <  8:31382/147‑1 (MQ=255)
                    gCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCt                                                                                                                                >  6:158693/1‑149 (MQ=255)
                    gCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGCCGAGTCGACCTTGctcat                                                                                                                                  <  1:98340/147‑5 (MQ=255)
                    gCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGCCGAGTCGACCTTGctcat                                                                                                                                  <  2:124303/147‑5 (MQ=255)
                     ccAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCtt                                                                                                                               >  3:64311/1‑149 (MQ=255)
                        gTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCTAGGTCATCCTCGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCa                                                                                                                            <  8:265308/149‑1 (MQ=255)
                           gCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAgg                                                                                                                          <  5:158693/148‑1 (MQ=255)
                            cccTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGc                                                                                                                         >  3:188562/1‑148 (MQ=255)
                             ccTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCaa                                                                                                                       >  2:67965/1‑149 (MQ=255)
                                 tACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGCCGAGTCGACCTTGctc                                                                                                                                    >  1:283864/1‑130 (MQ=255)
                                 tACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGCCGAGTCGACCTTGctc                                                                                                                                    <  2:283864/132‑3 (MQ=255)
                                    cTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAAccc                                                                                                               >  1:241130/1‑149 (MQ=255)
                                          gccgcTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCACCCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCACCCCGGAAGt                                                                                                         >  4:243577/1‑149 (MQ=255)
                                          gccgcTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACTACAGTGCCAGCGTTGGCATGAGCAAGTTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGTAAGt                                                                                                         <  6:354663/149‑1 (MQ=255)
                                                 gCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGcc                                                                                                    <  1:67965/147‑1 (MQ=255)
                                                   gTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGt                                                                                                >  5:175131/1‑149 (MQ=255)
                                                   gTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGt                                                                                                <  7:258838/149‑1 (MQ=255)
                                                         gAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCa                                                                                          >  1:343443/1‑149 (MQ=255)
                                                           gCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCTGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCAcc                                                                                        >  3:50171/1‑149 (MQ=255)
                                                           gCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGcc                                                                                                    <  1:227936/137‑1 (MQ=255)
                                                           gCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGcc                                                                                                    >  2:227936/1‑137 (MQ=255)
                                                            cTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCg                                                                                       >  8:185843/1‑149 (MQ=255)
                                                               cACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGcc                                                                                     <  2:241130/148‑1 (MQ=255)
                                                                       cGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTg                                                                             >  3:352112/1‑148 (MQ=255)
                                                                        gACTTCGACACCTTCGACCTGACAAACCTGCAACGCCAGTTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCTGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGa                                                                           <  7:214461/149‑1 (MQ=255)
                                                                         aCTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGAc                                                                          <  1:349628/149‑1 (MQ=255)
                                                                         aCTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGAc                                                                          <  3:243577/149‑1 (MQ=255)
                                                                                   cGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGCACTCGc                                                                >  3:295390/1‑149 (MQ=255)
                                                                                    gTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCt                                                               >  6:319641/1‑149 (MQ=255)
                                                                                       gACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCTAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGa                                                                        <  1:76685/138‑1 (MQ=255)
                                                                                       gACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGa                                                                        >  2:76685/1‑137 (MQ=255)
                                                                                       gACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGATGACTCGCtggt                                                            <  4:352112/149‑1 (MQ=255)
                                                                                        aCCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCtggtg                                                           <  3:268650/149‑1 (MQ=255)
                                                                                             aCCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGc                                                      >  7:292416/1‑149 (MQ=255)
                                                                                               cAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCtg                                                              >  1:267926/1‑139 (MQ=255)
                                                                                               cAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCtg                                                              <  2:267926/139‑1 (MQ=255)
                                                                                               cAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCgg                                                    <  3:339436/149‑1 (MQ=255)
                                                                                                  ccTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGGTGCGGTgg                                                 >  4:303104/1‑149 (MQ=255)
                                                                                                      cAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTggcggc                                             <  6:175131/149‑1 (MQ=255)
                                                                                                      cAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTggcgg                                              <  3:206610/148‑1 (MQ=255)
                                                                                                       aaCGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAg                                                                       <  2:256802/122‑1 (MQ=255)
                                                                                                       aaCGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAg                                                                       >  1:256802/1‑122 (MQ=255)
                                                                                                       aaCGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGgcggcg                                            <  5:319641/149‑1 (MQ=255)
                                                                                                            cAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGa                                       >  3:20332/1‑149 (MQ=255)
                                                                                                               gTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCt                                    >  7:63938/1‑149 (MQ=255)
                                                                                                                tCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTg                                   >  5:220478/1‑149 (MQ=255)
                                                                                                                  aTCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGt                                 >  2:194319/1‑149 (MQ=255)
                                                                                                                     cACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCt                              >  6:165072/1‑149 (MQ=255)
                                                                                                                      aCGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTg                             <  4:965/149‑1 (MQ=255)
                                                                                                                        gACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGa                           <  2:343443/149‑1 (MQ=255)
                                                                                                                         aCAGTGCCAGCGTGGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGAc                          <  4:64311/149‑1 (MQ=255)
                                                                                                                             tGCCAGCGTTGGCATGAGCAAGGTCGACTCGGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGc                       >  1:34028/1‑149 (MQ=255)
                                                                                                                                 aGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGa                  <  1:342027/149‑1 (MQ=255)
                                                                                                                                 aGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGa                  >  2:300744/1‑149 (MQ=255)
                                                                                                                                 aGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGa                  >  2:342027/1‑149 (MQ=255)
                                                                                                                                     ttGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAAt              <  8:63938/149‑1 (MQ=255)
                                                                                                                                        gCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAATTtc           >  1:117334/1‑149 (MQ=255)
                                                                                                                                          aTGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTgg                                  >  4:350918/1‑124 (MQ=255)
                                                                                                                                          aTGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTgg                                  <  3:350918/124‑1 (MQ=255)
                                                                                                                                             tgCAAGGCCGAGTCGACCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGTCGGTCGACCTGGTGCTGGACTGCTCCGACAATTTCTCTAc      >  6:349624/2‑149 (MQ=255)
                                                                                                                                              gCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAATTTCTCTAcc     <  5:165072/149‑1 (MQ=255)
                                                                                                                                                  ggTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGCCCTGGTGCTGGACTGCTCCGACAATTTCTCTACCcgc  >  6:337674/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGGCTGGGTGGCCTGGGCTCGCCAGTCGCCCTGTACCTGGCCGCCGCTGGCGTGGGTGAGCTGCACCTGGCCGACTTCGACACCGTCGACCTGACCAACCTGCAACGCCAGGTCATCCACGACAGTGCCAGCGTTGGCATGAGCAAGGTCGACTCGGCCTTGCAGCGCTTGCAGGCAATCAACCCGGAAGTCAGCCTGGTTGCTCACCGCCAGGCACTGGACGAGGACTCGCTGGTGGCTGCGGTGGCGGCGGTCGACCTGGTGCTGGACTGCTCCGACAATTTCTCTACCCGC  >  NC_002947/854507‑854800

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: