Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,005,680 G→A 7.0% H79Y (CAC→TAC) ‡ amaB ← L‑piperidine‑6‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,005,6800GA 7.0% 1144.6 / 49.8 371H79Y (CAC→TAC) ‡amaBL‑piperidine‑6‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (198/147);  new base A (16/10);  total (214/157)
Fisher's exact test for biased strand distribution p-value = 8.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 379 total aligned reads displayed.

GCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGCCGAGCAGTTTCACCGAGGCGATCTGGCTGC  >  NC_002947/6005531‑6005829
                                                                                                                                                     |                                                                                                                                                     
gcgcTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCATCGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTtgtg                                                                                                                                                       <  2:678201/150‑1 (MQ=255)
gcgcTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTtgtg                                                                                                                                                       <  2:2083714/150‑1 (MQ=255)
gcgcTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTtgtg                                                                                                                                                       <  2:6725279/150‑1 (MQ=255)
   cTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTc                                                                                                                                                    <  2:4983755/150‑1 (MQ=255)
    tCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCa                                                                                                                                                   >  2:7257725/1‑150 (MQ=255)
    tCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCa                                                                                                                                                   >  1:3742500/1‑150 (MQ=255)
    tCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCa                                                                                                                                                   >  2:5850308/1‑150 (MQ=255)
     cGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCAc                                                                                                                                                  >  1:4982955/1‑150 (MQ=255)
         aGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACgcag                                                                                                                                              <  1:4057124/150‑1 (MQ=255)
           cGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACgcagca                                                                                                                                            <  2:112864/150‑1 (MQ=255)
           cGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACgcagca                                                                                                                                            <  2:3456940/150‑1 (MQ=255)
                 tCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGc                                                                                                                                      >  1:8985940/1‑150 (MQ=255)
                    ggCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGa                                                                                                                                   >  2:5401862/1‑150 (MQ=255)
                     gCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCGTTGTGCTCACGCAGCACTTCGCCGaa                                                                                                                                  >  1:2385608/1‑150 (MQ=255)
                     gCCGTAACGCTGACGCGACCGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCGTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGCGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGag                                                                                                                                  <  1:6564972/150‑2 (MQ=255)
                      ccGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAAc                                                                                                                                 <  1:3253230/150‑1 (MQ=255)
                        gTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAg                                                                                                                               <  2:628784/150‑1 (MQ=255)
                        gTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAg                                                                                                                               <  1:9136848/150‑1 (MQ=255)
                        gTAAAGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAg                                                                                                                               >  2:2732252/1‑150 (MQ=255)
                         taaaGCTGACGCGACAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAgg                                                                                                                              <  1:5338015/147‑1 (MQ=255)
                              cTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCAc                                                                                                                         <  2:6332053/150‑1 (MQ=255)
                              cTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCAc                                                                                                                         <  1:4596497/150‑1 (MQ=255)
                              cTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCAc                                                                                                                         <  1:726692/150‑1 (MQ=255)
                                ggCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCa                                                                                                                       <  2:9007433/150‑1 (MQ=255)
                                 gcgcGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAg                                                                                                                      <  1:6878109/150‑1 (MQ=255)
                                     gACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCg                                                                                                                  <  1:8809961/150‑1 (MQ=255)
                                       cAGACCAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAGGCGGATCAGTTCGcc                                                                                                                <  1:486510/150‑1 (MQ=255)
                                         ggCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCAc                                                                                                              >  2:6946360/1‑150 (MQ=255)
                                          gCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACg                                                                                                             >  2:661891/1‑150 (MQ=255)
                                          gCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCGGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCGTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCCGTTCGCCACg                                                                                                             >  2:8899569/1‑150 (MQ=255)
                                          accAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGATTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAGGCGGATCAGTTCGCCACg                                                                                                             >  2:1227258/2‑150 (MQ=255)
                                            cAACTGCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAGGCGGATCAGTTCGCCACgct                                                                                                           >  1:5429020/1‑148 (MQ=255)
                                                tgCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAGGCGGATCAGTTCGCCACgctttgg                                                                                                       >  1:1112323/2‑144 (MQ=255)
                                                tgCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAGGCGGATCAGTTCGCCACgctttgg                                                                                                       >  1:8501728/2‑144 (MQ=255)
                                                tgCGAAATCGCAGATGTCGATCATTTCCTGTACTTCGCCGAGGCCTTCCGAAGGAATCTTGCCAGCTTCGAGCGAAACGAGACGGCCGAGGTCTTCCTTGGAAGCACGCAGTTCTTCACCGAGCAGGCGGATCAGTTCGCCACgctttgg                                                                                                       >  1:8053803/2‑144 (MQ=255)
                                                  cGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAg                                                                                                     >  2:7833141/1‑150 (MQ=255)
                                                    caGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGcc                                                                                                   <  2:6136108/149‑1 (MQ=255)
                                                       tCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGc                                                                                                >  2:8272086/1‑150 (MQ=255)
                                                         gCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGGTCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGcac                                                                                              >  2:7430297/1‑150 (MQ=255)
                                                         gCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGcac                                                                                              >  2:4336601/1‑150 (MQ=255)
                                                           aGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCacac                                                                                            <  2:2341594/150‑1 (MQ=255)
                                                              tGTCGATCATTTCCTGCACTTCACCCAGGCGTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTgc                                                                                         <  2:188966/150‑1 (MQ=255)
                                                                  gATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCa                                                                                     <  2:8492101/150‑1 (MQ=255)
                                                                  gATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCa                                                                                     >  1:4517229/1‑150 (MQ=255)
                                                                   aTCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAg                                                                                    >  1:6693343/1‑150 (MQ=255)
                                                                       tttCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCtt                                                                                >  2:4427807/1‑150 (MQ=255)
                                                                         tCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCg                                                                              >  1:3554512/1‑150 (MQ=255)
                                                                          ccTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGa                                                                             >  1:7455613/1‑150 (MQ=255)
                                                                           cTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGaa                                                                            <  1:836361/150‑1 (MQ=255)
                                                                            tGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAg                                                                           >  1:2221255/1‑150 (MQ=255)
                                                                            tGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAg                                                                           >  1:4192324/1‑150 (MQ=255)
                                                                               aCTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGgcg                                                                        >  1:1864007/1‑150 (MQ=255)
                                                                                gatCACCCAGGAAGGCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACCGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCCGGCTTCGAAGgcgc                                                                       >  1:3466592/3‑150 (MQ=255)
                                                                                 ttCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGct                                                                      >  1:2750159/1‑150 (MQ=255)
                                                                                 ttCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGct                                                                      >  2:798477/1‑150 (MQ=255)
                                                                                 ttCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGct                                                                      >  1:3135282/1‑150 (MQ=255)
                                                                                   cccccAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGctct                                                                    <  2:9066150/148‑1 (MQ=255)
                                                                                      ccAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTggg                                                                 >  2:1664238/1‑150 (MQ=255)
                                                                                      ccAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTggg                                                                 >  1:1404072/1‑150 (MQ=255)
                                                                                      ccAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTggg                                                                 >  2:4135394/1‑150 (MQ=255)
                                                                                      ccAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTggg                                                                 >  2:318610/1‑150 (MQ=255)
                                                                                      ccAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTggg                                                                 >  2:3109379/1‑150 (MQ=255)
                                                                                       cAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGc                                                                <  1:3005162/150‑1 (MQ=255)
                                                                                            cTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGAt                                                           >  1:6296694/1‑150 (MQ=255)
                                                                                            cTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGAt                                                           >  1:6215664/1‑150 (MQ=255)
                                                                                               ccTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGa                                                        >  2:1065951/1‑150 (MQ=255)
                                                                                                 tGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATg                                                      >  2:3898517/1‑150 (MQ=255)
                                                                                                   ggTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCg                                                    <  1:7594715/150‑1 (MQ=255)
                                                                                                          ttGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATg                                             >  2:2935621/1‑150 (MQ=255)
                                                                                                          ttGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATg                                             >  1:5189231/1‑150 (MQ=255)
                                                                                                          ttGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATg                                             >  2:8741146/1‑150 (MQ=255)
                                                                                                            gCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGt                                           >  2:4618974/1‑150 (MQ=255)
                                                                                                                  ttCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGc                                     <  2:7708379/150‑1 (MQ=255)
                                                                                                                  ttCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGc                                     >  1:10147/1‑150 (MQ=255)
                                                                                                                  ttCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGc                                     >  1:5634737/1‑150 (MQ=255)
                                                                                                                     gatcgaGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCtt                                  >  1:4679071/1‑150 (MQ=255)
                                                                                                                     gatcgaGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCtt                                  <  1:7933350/150‑1 (MQ=255)
                                                                                                                     gatcgaGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCtt                                  >  2:8861174/1‑150 (MQ=255)
                                                                                                                       tcgaGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTg                                 <  1:2610947/149‑1 (MQ=255)
                                                                                                                       tcgaGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGc                                >  2:6270283/1‑150 (MQ=255)
                                                                                                                       tcgaGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGc                                <  2:6473716/150‑1 (MQ=255)
                                                                                                                       tcgaGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGc                                >  1:5974413/1‑150 (MQ=255)
                                                                                                                         gagaCCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGCCg                              >  1:7754941/1‑150 (MQ=255)
                                                                                                                          agaCCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGGGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGCCGa                             <  1:2876910/150‑1 (MQ=255)
                                                                                                                          agaCCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGCCGa                             <  1:4430774/150‑1 (MQ=255)
                                                                                                                          agaCCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGCCGa                             <  1:5850308/150‑1 (MQ=255)
                                                                                                                              cAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGCCGAGCAg                         >  2:4872825/1‑150 (MQ=255)
                                                                                                                              cAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGCCGAGCAg                         >  1:3831484/1‑150 (MQ=255)
                                                                                                                                  tCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGCCGAGCAGTTTc                     <  1:1758509/150‑1 (MQ=255)
                                                                                                                                  tCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGCCGAGCAGTTTc                     >  2:3836445/1‑150 (MQ=255)
                                                                                                                                      cGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGCCGAGCAGTTTCACCg                 >  1:6853205/1‑150 (MQ=255)
                                                                                                                                       gAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGCCGAGCAGTTTCACCGa                <  2:5625327/150‑1 (MQ=255)
                                                                                                                                        aGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGCCGAGCAGTTTCACCGAg               <  1:6422639/150‑1 (MQ=255)
                                                                                                                                          gTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGCCGAGCAGTTTCACCGAGGc             <  1:8564838/150‑1 (MQ=255)
                                                                                                                                          gTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGCCGAGCAGTTTCACCGAGGc             >  2:5532714/1‑150 (MQ=255)
                                                                                                                                          gTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGCCGAGCAGTTTCACCGAGGc             >  1:3890105/1‑150 (MQ=255)
                                                                                                                                             ggCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGCCGAGCAGTTTCACCGAGGCGAt          >  1:1512780/1‑150 (MQ=255)
                                                                                                                                                 ttGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGCCGAGCAGTTTCACCGAGGCGATCTgg      <  2:2215712/150‑1 (MQ=255)
                                                                                                                                                 ttGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGCCGAGCAGTTTCACCGAGGCGATCTgg      <  2:2507337/150‑1 (MQ=255)
                                                                                                                                                     gCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGCCGAGCAGTTTCACCGAGGCGATCTGGCTGc  >  2:5204765/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGCTCGGAAGCGATGGTCAGGCCGTACAGCTGGCGCGACAGGCCAACGGCGAAGTCGCAGATGTCGATCATTTCCTGCACTTCACCCAGGCCTTCCTGGGTGATCTTGCCGGCTTCGATCGAGACCAGCTCGCCGAGGTCGGCCTTGTGCTCACGCAGCACTTCGCCGAACAGGCGCACCAGTTCGCCACGGCGCGGAGCCGGCACACTGCGCCAGGCTTCGAAGGCGCTCTGGGCCTGATCGATGCGGGCGATGGTCTCGGCCTTGCCGAGCAGTTTCACCGAGGCGATCTGGCTGC  >  NC_002947/6005531‑6005829

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