Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,657,944 C→A 6.1% intergenic (‑199/‑53) tktA ← / → sahR transketolase A/methionine metabolism transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,657,9440CA 6.1% 422.8 / 8.7 130intergenic (‑199/‑53)tktA/sahRtransketolase A/methionine metabolism transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (63/59);  new base A (5/3);  total (68/62)
Fisher's exact test for biased strand distribution p-value = 7.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 131 total aligned reads displayed.

CATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTG  >  NC_002947/5657795‑5658089
                                                                                                                                                     |                                                                                                                                                 
cATTTTCACTCACTTCGGGCGCTTGTAGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAACGTGTAACTTGAGGCTTTCGGCTTGCAGCCCATCGCCAATACCAAAAc                                                                                                                                                   >  2:228126/1‑150 (MQ=255)
     tCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACttttt                                                                                                                                              <  2:312215/150‑1 (MQ=255)
     tCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACttttt                                                                                                                                              <  2:524050/150‑1 (MQ=255)
           aCTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATAtt                                                                                                                                        >  2:1923170/1‑150 (MQ=255)
           aCTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATAtt                                                                                                                                        <  1:1919402/150‑1 (MQ=255)
            cTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTg                                                                                                                                       >  1:2039270/1‑150 (MQ=255)
             tGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCGGCTCATCGCCAATATCAAAAAGTTTTGATATTgg                                                                                                                                      <  1:1001460/150‑1 (MQ=255)
               ggggCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGAGACTCGGGGCTCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGcc                                                                                                                                    >  2:1288801/1‑150 (MQ=255)
                 ggCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCtt                                                                                                                                  >  2:504050/1‑150 (MQ=255)
                  gcgcTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGTCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTg                                                                                                                                 >  1:1741631/1‑150 (MQ=255)
                  gcgcTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTg                                                                                                                                 >  2:303132/1‑150 (MQ=255)
                      ttGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAg                                                                                                                             >  1:1767678/1‑150 (MQ=255)
                       tGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGTCTCGGGGCGCCACAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTatg                                                                                                                            <  2:926190/150‑3 (MQ=255)
                       tGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAgg                                                                                                                            <  1:526284/150‑1 (MQ=255)
                         gggCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGc                                                                                                                          <  2:244614/150‑1 (MQ=255)
                         gggCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGc                                                                                                                          >  1:352692/1‑150 (MQ=255)
                          ggCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCa                                                                                                                         >  2:2333461/1‑150 (MQ=255)
                           gCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAgg                                                                                                                            >  2:1196328/1‑146 (MQ=255)
                           gCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCAt                                                                                                                        >  1:1945027/1‑150 (MQ=255)
                           gCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCAt                                                                                                                        >  1:277017/1‑150 (MQ=255)
                            cAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTTGCCTTGCTAGAGCATc                                                                                                                       >  2:1177082/1‑150 (MQ=255)
                             aaaCATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCCAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTATGGCATCg                                                                                                                      >  2:539236/1‑150 (MQ=255)
                             aaaCATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCg                                                                                                                      >  1:1300873/1‑150 (MQ=255)
                              aaCATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGa                                                                                                                     >  1:1726890/1‑150 (MQ=255)
                              aaCATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGa                                                                                                                     >  2:419300/1‑150 (MQ=255)
                               aCATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGAt                                                                                                                    <  2:1726890/150‑1 (MQ=255)
                               aCATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGAt                                                                                                                    <  1:1139283/150‑1 (MQ=255)
                                cATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATg                                                                                                                   <  2:1866875/150‑1 (MQ=255)
                                   ttAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCt                                                                                                                <  1:1204122/150‑1 (MQ=255)
                                     aagcgccaagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGt                                                                                                              <  1:1698977/150‑1 (MQ=255)
                                      agcgccaagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTc                                                                                                             <  2:352692/150‑1 (MQ=255)
                                        cgccaagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTa                                                                                                           <  1:1345994/150‑1 (MQ=255)
                                         gccaagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGa                                                                                                                     >  1:153333/1‑139 (MQ=255)
                                         gccaagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGa                                                                                                                     <  2:153333/139‑1 (MQ=255)
                                          ccaagcgccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGa                                                                                                         <  2:2223983/150‑1 (MQ=255)
                                                 ccGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCt                                                                                                                >  2:506279/1‑136 (MQ=255)
                                                  cGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGAGATGCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTgccgcc                                                                                                 >  1:1171716/1‑150 (MQ=255)
                                                     cTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCa                                                                                              >  1:1146165/1‑150 (MQ=255)
                                                          aGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGaag                                                                                         <  1:977011/150‑2 (MQ=255)
                                                           gCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGAc                                                                                                        >  1:330684/1‑134 (MQ=255)
                                                           gCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAAAGTTTTGATATTgg                                                                                                                                      <  2:1402322/104‑1 (MQ=255)
                                                           gCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAAAGTTTTGATATTgg                                                                                                                                      >  1:1402322/1‑104 (MQ=255)
                                                            ccGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCt                                                                                       >  1:1422728/1‑150 (MQ=255)
                                                              gCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCtgc                                                                                     <  1:1583721/150‑1 (MQ=255)
                                                                aaCTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAgg                                                                                                                            <  2:1178898/109‑1 (MQ=255)
                                                                aaCTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAgg                                                                                                                            >  1:1178898/1‑109 (MQ=255)
                                                                  cTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCACCATGAATCtgcgtgc                                                                                 <  1:1288801/150‑1 (MQ=255)
                                                                   tGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGaa                                                                                          >  2:744059/1‑140 (MQ=255)
                                                                   tGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCtgcgtgcg                                                                                <  2:1146165/150‑1 (MQ=255)
                                                                   tGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCtgcgtgcg                                                                                >  2:2066323/1‑150 (MQ=255)
                                                                      agaagCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCaa                                                                             <  2:1741631/150‑1 (MQ=255)
                                                                       gaagCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAATCTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCt                                                                                                                <  1:1068258/114‑1 (MQ=255)
                                                                       gaagCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCt                                                                                                                >  2:1068258/1‑114 (MQ=255)
                                                                          gCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGa                                                                         >  2:940301/1‑150 (MQ=255)
                                                                             gagCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATcc                                                                      <  1:504050/150‑1 (MQ=255)
                                                                               gcgcGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATcccc                                                                    <  2:2042048/150‑1 (MQ=255)
                                                                                   gACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGc                                                                >  1:931971/1‑150 (MQ=255)
                                                                                   gACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGc                                                                >  1:3992/1‑150 (MQ=255)
                                                                                      tCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAAc                                                             <  1:1831694/150‑1 (MQ=255)
                                                                                        ggggCGCCCCAACGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAAGATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAAcgc                                                           <  1:2299593/150‑1 (MQ=255)
                                                                                          ggCGCCCCAAAGCGTGTAACTTGCGGCTTTGGGCTTGCAGCTCATAGCCAATAGCAAAACTTTTTGATATTGGCCTTACTAGGACATCGATTCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGAAACGCGa                                                         >  1:1156525/1‑150 (MQ=255)
                                                                                               cccAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACAt                                                    >  2:1515018/1‑150 (MQ=255)
                                                                                                ccAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACAtt                                                   <  1:228126/150‑1 (MQ=255)
                                                                                                ccAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACAtt                                                   <  2:931971/150‑1 (MQ=255)
                                                                                                 cAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAg                                                                                                                             <  1:478487/75‑1 (MQ=255)
                                                                                                  aaaGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCtgcgt                                                                                   >  1:828858/1‑116 (MQ=255)
                                                                                                  aaaGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCtgcgt                                                                                   <  2:828858/116‑1 (MQ=255)
                                                                                                  aaaGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTgg                                                 <  2:2039270/150‑1 (MQ=255)
                                                                                                    aGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCa                                               <  1:539236/150‑1 (MQ=255)
                                                                                                    aGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCACCAATATCAAAACTTTTTGATATTGGCCTGGGCGCTGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCa                                               <  1:419300/150‑1 (MQ=255)
                                                                                                         gtAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCt                                          >  1:253770/1‑150 (MQ=255)
                                                                                                           aaCTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTTGCAGCCCtgt                                        >  2:145661/1‑150 (MQ=255)
                                                                                                           aaCTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCtgt                                        <  1:2157725/150‑1 (MQ=255)
                                                                                                           aaCTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCtgt                                        <  1:2158127/150‑1 (MQ=255)
                                                                                                             cTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGc                                      <  2:1314590/150‑1 (MQ=255)
                                                                                                              ttGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCa                                     >  2:830345/1‑150 (MQ=255)
                                                                                                              ttGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCa                                     >  2:1630244/1‑150 (MQ=255)
                                                                                                                gAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCaaa                                   <  1:1515018/150‑1 (MQ=255)
                                                                                                                gAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCaaa                                   <  1:2066323/150‑1 (MQ=255)
                                                                                                                          gCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCg                         >  1:2285001/1‑150 (MQ=255)
                                                                                                                           cTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGa                        <  1:863473/150‑1 (MQ=255)
                                                                                                                            ttGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTgg                           >  2:2051087/1‑146 (MQ=255)
                                                                                                                            ttGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGAt                       >  1:65112/1‑150 (MQ=255)
                                                                                                                                 gCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGAtgcgct                  >  1:7027/1‑150 (MQ=255)
                                                                                                                                  cTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCACTGTGCAAAGCCAGTGGCGAtgcgctg                 >  2:1101581/1‑150 (MQ=255)
                                                                                                                                   tCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTg                                                  >  1:755836/1‑116 (MQ=255)
                                                                                                                                     aTCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGAtgcgctgcgc              <  2:1152545/150‑1 (MQ=255)
                                                                                                                                       cGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCt            >  2:436990/1‑150 (MQ=255)
                                                                                                                                        gCCAATATCAAAAAGTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTGGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTg           >  2:2359481/1‑150 (MQ=255)
                                                                                                                                        gCCAATATCAAAAAGTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTg           >  1:1983864/1‑150 (MQ=255)
                                                                                                                                        gCCAATATCAAAAAGTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTg           >  1:2374069/1‑150 (MQ=255)
                                                                                                                                         ccAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATc                                                                       <  1:1543498/89‑1 (MQ=255)
                                                                                                                                         ccAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATc                                                                       >  2:1543498/1‑89 (MQ=255)
                                                                                                                                          cAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAAcgc                                                           <  2:279427/100‑1 (MQ=255)
                                                                                                                                          cAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCt                                          >  1:2228425/1‑117 (MQ=255)
                                                                                                                                          cAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCt                                          <  2:2228425/117‑1 (MQ=255)
                                                                                                                                          cAATATCAAAAAGTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGaa         <  2:1983864/150‑1 (MQ=255)
                                                                                                                                            atatCAAAAAGTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACg       >  2:443602/1‑150 (MQ=255)
                                                                                                                                                cAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTAtt   <  2:65112/150‑1 (MQ=255)
                                                                                                                                                 aaaaCTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTg  >  2:830564/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CATTTTCGCCCACTGGGGGCGCTTGGGGCAAACATTTAAGCGCCAAGCGCCGGCTTTAAGCCGCAACTGAAGAAGCAGAGCGCGACTCGGGGCGCCCCAAAGCGTGTAACTTGAGGCTTTCGGCTTGCAGCTCATCGCCAATATCAAAACTTTTTGATATTGGCCTTGCTAGGGCATCGATGCCTGTCTAGACTGCCGCCCCATGAATCTGCGTGCGCAATCGATCCCCGAGCAACGCGAAACATTGGCAGCCCTGTGCAAAGCCAGTGGCGATGCGCTGCGCCTGAACGTATTG  >  NC_002947/5657795‑5658089

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