Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,204,963 G→A 6.6% W123* (TGG→TAG)  yehL → putative hexameric AAA+ MoxR family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,204,9630GA 6.6% 193.5 / 2.4 60W123* (TGG→TAG) yehLputative hexameric AAA+ MoxR family ATPase
Reads supporting (aligned to +/- strand):  ref base G (24/32);  new base A (2/2);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGGCCGCGGTTTGATGTTGGTGGGGGAACCGGGTACCGCAAAATCTCTTCTTTCTGAATTACTGGCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAACGCGTACACCGCTGGAAGTGCAGGACTGTCTGTTAG  >  NC_000913/2204818‑2205106
                                                                                                                                                 |                                                                                                                                               
aCCGGCCGCGGTTTGATGTTGGTGGGGGAACCGGGTACCGCAAAATCTCTTCTTTCTGAATTACTGGCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTgg                                                                                                                                                >  7:325290/1‑147 (MQ=255)
           tttGATGTTGGTGGGGGAACCGGGTACCGCAAAATCTCTTCTTTCTGAATTACTGGCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTg                                                                                                                                    <  4:250548/148‑1 (MQ=255)
           tttGATGTTGGTGGGGGAACCGGGTACCGCAAAATCTCTTCTTTCTGAATTACTGGCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTg                                                                                                                                    <  3:99457/148‑1 (MQ=255)
            ttGATGTTGGTGGGGGAACCGGGTACCGCAAAATCTCTTCTTTCTGAATTACTGGCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGc                                                                                                                                   >  2:326171/1‑148 (MQ=255)
                     gtggGGGAACCGGGTACCGCAAAATCTCTTCTTTCTGAATTACTGGCACCCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTatcaatc                                                                                                                          <  1:194758/148‑1 (MQ=255)
                     gtggGGGAACCGGGTACCGCAAAATCTCTTCTTTCTGAATTACTGGCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTatcaat                                                                                                                           <  5:253017/147‑1 (MQ=255)
                       gggggAACCGGGTACCGCAAAATCTCTTCTTTCTGAATTACTGGCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCAt                                                                                                                        <  8:302406/148‑1 (MQ=255)
                           gAACCGGGTACCGCAAAATCTCGTCTTTCTGAATTACTGGCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGCACTACGCCCTGCTTATCAATCATGGTc                                                                                                                    >  2:244069/1‑148 (MQ=255)
                                        cAAAATCTCTTCTTTCTGAATTACTGGCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCCCCGAAGATCAAAACAAGTATGGCTGCAACTACGCCCTGCTTATCAATCATGGTCCGACAACGGAAgc                                                                                                       >  3:286089/1‑148 (MQ=255)
                                           aaTCTCTTCTTTCTGAATTACTGGCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCt                                                                                                    <  7:167847/148‑1 (MQ=255)
                                            aTCTCTTCTTTCTGAATTACTGGCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTg                                                                                                   <  8:325290/148‑1 (MQ=255)
                                                     ttctGAATTACTGGCAACCGCTATTTGGGGGGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGc                                                                                           <  4:286089/147‑1 (MQ=255)
                                                        tGAATTACTGGCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTAGa                                                                                                                                               <  5:254523/92‑1 (MQ=255)
                                                        tGAATTACTGGCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTAGa                                                                                                                                               <  5:250295/92‑1 (MQ=255)
                                                        tGAATTACTGGCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTAGa                                                                                                                                               >  6:250295/1‑92 (MQ=255)
                                                        tGAATTACTGGCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTAGa                                                                                                                                               >  6:254523/1‑92 (MQ=255)
                                                          aaTTACTGGCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCt                                                                                     >  4:67585/1‑148 (MQ=255)
                                                           aTTACTGGCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTAAGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCtt                                                                                    >  4:275696/1‑148 (MQ=255)
                                                                tGGCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCtt                                                                                                                                 <  2:117890/98‑1 (MQ=255)
                                                                tGGCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCtt                                                                                                                                 >  1:117890/1‑98 (MQ=255)
                                                                 ggCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAg                                                                              >  7:66852/1‑148 (MQ=255)
                                                                  gCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAg                                                                              <  8:66852/147‑1 (MQ=255)
                                                                  gCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTAAGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAg                                                                              <  3:275696/147‑1 (MQ=255)
                                                                    aaCCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGa                                                                           <  5:152195/148‑1 (MQ=255)
                                                                        gCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGGCCCTGCCCCGCTTTATCAGGGAATgc                                                                       <  1:274647/148‑1 (MQ=255)
                                                                         cTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCATCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATgcg                                                                      >  3:214882/1‑148 (MQ=255)
                                                                               gcggcgATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCACGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACg                                                                 <  1:244069/147‑1 (MQ=255)
                                                                                cggcgATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACgg                                                                >  7:176152/1‑147 (MQ=255)
                                                                                            gTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGt                                                   >  6:305832/1‑148 (MQ=255)
                                                                                                aCCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTc                                                >  8:296141/1‑147 (MQ=255)
                                                                                                    tCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGaa                                             >  2:305954/1‑146 (MQ=255)
                                                                                                        gggCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATa                                       <  3:41871/148‑1 (MQ=255)
                                                                                                          gCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAAc                                     >  5:18427/1‑148 (MQ=255)
                                                                                                           cGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAAcg                                    <  8:176152/148‑1 (MQ=255)
                                                                                                                   cTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATa                                       >  8:313417/1‑137 (MQ=255)
                                                                                                                   cTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATa                                       <  7:313417/137‑1 (MQ=255)
                                                                                                                     accaccGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATa                                       >  1:60229/1‑135 (MQ=255)
                                                                                                                     accaccGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATa                                       <  2:60229/135‑1 (MQ=255)
                                                                                                                     accaccGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAACGCGTACACCGc                          >  3:168949/1‑148 (MQ=255)
                                                                                                                       caccGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAACGCGTACACCGCTg                        <  1:305954/148‑1 (MQ=255)
                                                                                                                         ccGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAACGCGTACACCGCTgg                       >  3:159745/1‑147 (MQ=255)
                                                                                                                         ccGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAACGCGTACACCGCTGGa                      <  3:100117/148‑1 (MQ=255)
                                                                                                                         ccGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAACGCGTACACCGCTGGa                      <  4:168949/148‑1 (MQ=255)
                                                                                                                         ccGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGGCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAACGCGTACACCGCTGGa                      <  1:326171/148‑1 (MQ=255)
                                                                                                                            aaGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAACGCGTACACCGCTGGAAGt                   >  6:234624/1‑148 (MQ=255)
                                                                                                                              gATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAACGCGTACACCGCTGGAAGTGc                 <  7:296141/148‑1 (MQ=255)
                                                                                                                              gATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAACGCGTACACCGCTGGAAGTGc                 <  7:40121/148‑1 (MQ=255)
                                                                                                                              gATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCATCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAACGCGTACACCGCTGGAAGTGc                 <  4:214882/148‑1 (MQ=255)
                                                                                                                               atcaaatcaaGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAACGCGTACACCGCTGGAAGTg                  >  8:203909/1‑146 (MQ=255)
                                                                                                                               atcaaatcaaGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAACGCGTACACCGCTGGAAGTGCa                <  5:126585/148‑1 (MQ=255)
                                                                                                                                tcaaatcaaGTATGGCTGGAACTACGCCCTGCTTAACAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCATTATCAGGGAATGCGCGACGGCAAAATCGGCCGGTTCGAAGAGATAACGCGTACACCGCAGGAAGTGCAg               >  6:126585/1‑148 (MQ=255)
                                                                                                                                 caaatcaaGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAACGCGTACACCGCTGGAAGTGCAg               >  7:297863/1‑147 (MQ=255)
                                                                                                                                  aaatcaaGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAACGCGTACACCGCTGGAAGTGCAgg              <  7:203909/147‑1 (MQ=255)
                                                                                                                                  aaatcaaGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAACGCGTACACCGCTGGAAGTGCAgg              <  3:222419/147‑1 (MQ=255)
                                                                                                                                    atcaaGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGCCGGCAAAATCGTCCGTTTCGAAGAGATAACGCGTACACCGCTGGAAGTGCTGGAct           >  5:133952/1‑148 (MQ=255)
                                                                                                                                             ggCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAACGCGt                                 >  6:53293/1‑117 (MQ=255)
                                                                                                                                             ggCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAACGCGt                                 <  5:53293/117‑1 (MQ=255)
                                                                                                                                              gCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCACAATCGTCCGTTTCGAAGAGATAACGCGTACACCGCTGGAAGTGCAGGACTGTCTGTTAg  <  4:159745/147‑1 (MQ=255)
                                                                                                                                              gCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAACGCGTACACCGCTGGAAGTGCAGGACTGTCTGTTAg  <  5:305832/147‑1 (MQ=255)
                                                                                                                                               cTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAACGCGTACACCGCTGGAAg                    <  7:4396/128‑1 (MQ=255)
                                                                                                                                               cTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAACGCGTACACCGCTGGAAg                    >  8:4396/1‑128 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
ACCGGCCGCGGTTTGATGTTGGTGGGGGAACCGGGTACCGCAAAATCTCTTCTTTCTGAATTACTGGCAACCGCTATTAGCGGCGATGCCGGGTTAACCATCCAGGGCGGGGCATCTACCACCGAAGATCAAATCAAGTATGGCTGGAACTACGCCCTGCTTATCAATCATGGTCCGTCAACGGAAGCGCTGGTCCCTGCCCCGCTTTATCAGGGAATGCGCGACGGCAAAATCGTCCGTTTCGAAGAGATAACGCGTACACCGCTGGAAGTGCAGGACTGTCTGTTAG  >  NC_000913/2204818‑2205106

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