Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,489,582 G→T 6.0% T193T (ACG→ACT PP_2184 → formate dehydrogenase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,489,5820GT 6.0% 215.8 / 3.8 67T193T (ACG→ACTPP_2184formate dehydrogenase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (26/37);  new base T (2/2);  total (28/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01

CGAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTCG  >  NC_002947/2489435‑2489730
                                                                                                                                                    |                                                                                                                                                    
cGAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTACCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGt                                                                                                                                                     >  1:82132/1‑149 (MQ=255)
 gAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGtt                                                                                                                                                    <  5:319130/149‑1 (MQ=255)
 gAAGTGGTTGCCGCAGTACTCGACTCAGGCCTCCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCGGGGCAGAAGTACGTGGTGTGAAACGCCGACGAGGGTGACTCGGGCACGtt                                                                                                                                                    >  2:67735/1‑149 (MQ=255)
  aaGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTc                                                                                                                                                   <  6:66254/149‑1 (MQ=255)
    gTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGc                                                                                                                                                 <  2:82132/149‑1 (MQ=255)
           ccGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCg                                                                                                                                           >  6:44550/1‑148 (MQ=255)
           ccGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGc                                                                                                                                          <  5:44550/149‑1 (MQ=255)
                gTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCt                                                                                                                                     <  3:257610/149‑1 (MQ=255)
                    tCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGc                                                                                                                                      >  5:236665/1‑144 (MQ=255)
                    tCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGc                                                                                                                                      <  6:236665/144‑1 (MQ=255)
                    tCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATg                                                                                                                                 >  4:362196/1‑149 (MQ=255)
                      gACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGTGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATgg                                                                                                                                >  8:304363/1‑148 (MQ=255)
                      gACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGc                                                                                                                                          <  8:158565/138‑1 (MQ=255)
                      gACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGc                                                                                                                                          >  7:158565/1‑138 (MQ=255)
                         tCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAgg                                                                                                                            >  8:162079/1‑149 (MQ=255)
                           cGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGc                                                                                                                           >  1:326336/1‑148 (MQ=255)
                             gCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGTGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGAc                                                                                                                        <  7:304363/149‑1 (MQ=255)
                             gCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGAc                                                                                                                        >  2:387230/1‑149 (MQ=255)
                             gCCTGCGTGGCCGTGGTGGGGCGGCGTTCCAGGCCGGCATCAAGTGGCGTACCGTGCGGGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGGGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGAc                                                                                                                        <  4:5658/149‑1 (MQ=255)
                              ccTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGAcc                                                                                                                       <  3:105247/149‑1 (MQ=255)
                               cTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGAccc                                                                                                                      <  3:328564/149‑1 (MQ=255)
                                   gTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTc                                                                                                                  <  5:254963/149‑1 (MQ=255)
                                      gCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCctg                                                                                                               <  1:387230/149‑1 (MQ=255)
                                       ccGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGCCTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCctg                                                                                                               <  1:299692/149‑1 (MQ=255)
                                         gtggtggtgCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCctgct                                                                                                             >  7:166195/1‑148 (MQ=255)
                                               gggcggcgTTCCCGGCAGGCATCAAGTGGCGTACCGGGCGTGCTGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGa                                                                                                       <  6:33956/146‑1 (MQ=255)
                                                          ccggccggCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGAGGTGTGCAACGCCGCCGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCAt                                                                                           >  5:277309/1‑149 (MQ=255)
                                                           cggccggcATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCAtt                                                                                          >  6:1279/1‑149 (MQ=255)
                                                           cggccggcATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCAtt                                                                                          >  7:249945/1‑149 (MQ=255)
                                                           cggccggcATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCAtt                                                                                          >  2:32693/1‑149 (MQ=255)
                                                                  cATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGcc                                                                                   <  8:249945/149‑1 (MQ=255)
                                                                   ataaaGTGGCGTACCGTGCGTGATGCTGGCGCCCGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTCACTCCGACACGTTCGCTGATCGCATCCTGATGGAAGGCGACCCCTTCCTGCTGATACAGGACATGATCATTGCAGTCCt                                                                                  >  8:100301/4‑149 (MQ=255)
                                                                   aTCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCt                                                                                  >  1:218079/1‑149 (MQ=255)
                                                                     cAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCg                                                                                >  6:289470/1‑149 (MQ=255)
                                                                         tGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGCGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGt                                                                            <  2:9287/149‑1 (MQ=255)
                                                                         tGGCGGACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGt                                                                            <  2:15473/149‑1 (MQ=255)
                                                                             gTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGt                                                                        <  7:162079/149‑1 (MQ=255)
                                                                                ccGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGcc                                                                     >  6:135230/1‑149 (MQ=255)
                                                                                 cgtgcgtgATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGTGCATGATCATTGCCGGCCTCGCCGTGGGTGCCg                                                                    <  2:326336/149‑1 (MQ=255)
                                                                                  gtgcgtgATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCTGGCCTCGCCGTGGGTGCCGa                                                                   <  8:166195/149‑1 (MQ=255)
                                                                                   tgcgtgATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGCCCCCTTCCTGATGATCGAGGGCATGATCATTTCCGGCCTCGCCGTGGGTGCCGa                                                                   >  6:369161/1‑148 (MQ=255)
                                                                                     cgtgAGGCTGGCGCCGGGCAGACGTACGTGGTGTGCATCGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACaa                                                                <  6:204573/149‑1 (MQ=255)
                                                                                      gtgATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAg                                                               <  3:362196/149‑1 (MQ=255)
                                                                                        gATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACTTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGACCATTGCCGGCCTCGCCGTGGGTGCCGACAAggg                                                             >  4:82996/1‑149 (MQ=255)
                                                                                        gATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACTTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGACCATTGCCGGCCTCGCCGTGGGTGCCGACAAgg                                                              >  4:88320/1‑148 (MQ=255)
                                                                                         aTGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGc                                                            <  5:289470/149‑1 (MQ=255)
                                                                                            cTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCta                                                          >  2:313301/1‑148 (MQ=255)
                                                                                               gcgcCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTatat                                                       <  5:135230/148‑1 (MQ=255)
                                                                                               gcgcCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTatat                                                       <  2:218079/148‑1 (MQ=255)
                                                                                               gcgcCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATc                                                      >  4:182395/1‑149 (MQ=255)
                                                                                                cgcCGGGCAGAAGTACGTGGTGTGCATCGCCGACGAGGGTGACTCGGGCAAGTTCGCTGATCTCATGCTGTTGGAAGGCGACGCCTTCCTGCTGTTGGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCt                                                     <  7:100301/149‑1 (MQ=255)
                                                                                                  ccGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTAt                                                   <  1:32693/149‑1 (MQ=255)
                                                                                                          aaGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGTGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTc                                           >  4:281262/1‑149 (MQ=255)
                                                                                                             tACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACTTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGACCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGa                                        <  3:82996/149‑1 (MQ=255)
                                                                                                             tACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACTTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGACCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGa                                        <  3:88320/149‑1 (MQ=255)
                                                                                                                 tgttgagCAACGCCGACGAGGGTGACTAGGGCACGATCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGGGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTAt                                    <  8:305142/143‑1 (MQ=255)
                                                                                                                     gtgCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATccc                                 >  3:313174/1‑148 (MQ=255)
                                                                                                                      tgCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGa                               <  5:1279/149‑1 (MQ=255)
                                                                                                                        cAACGGCGACGAGGGTGACTCGGGCTCGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTTATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGTTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGAt                              <  1:67735/148‑1 (MQ=255)
                                                                                                                        cAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGAt                              <  1:165553/148‑1 (MQ=255)
                                                                                                                              cgacgaGGGGGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATcc                       <  8:255254/149‑1 (MQ=255)
                                                                                                                                       gTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTc               <  4:313174/149‑1 (MQ=255)
                                                                                                                                       gTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTc               <  1:313301/149‑1 (MQ=255)
                                                                                                                                            tCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGAt                              <  2:393607/129‑1 (MQ=255)
                                                                                                                                            tCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGAt                              >  1:393607/1‑129 (MQ=255)
                                                                                                                                              gggCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGTCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAg         <  6:277309/148‑1 (MQ=255)
                                                                                                                                                gCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGcc      >  8:253630/1‑149 (MQ=255)
                                                                                                                                                    gTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTCg  >  5:177141/1‑149 (MQ=255)
                                                                                                                                                    gTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTCg  >  8:33671/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGAAGTGGTTGCCGCAGTACTCGACTCCGGCCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTCG  >  NC_002947/2489435‑2489730

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: