Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,206,4350CT49.1% 39.9 / 43.9 57intergenic (‑179/‑28)ttdR/ttdAtranscriptional activator of ttdABT/L‑tartrate dehydratase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (2/27);  new base T (28/0);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 3.09e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCT  >  NC_000913/3206286‑3206581
                                                                                                                                                      |                                                                                                                                                  
aGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTcc                                                                                                                                                   <  2:663475/151‑1 (MQ=255)
   gtATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAg                                                                                                                                                 <  2:615261/150‑1 (MQ=255)
   gtATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGa                                                                                                                                                <  2:955769/151‑1 (MQ=255)
     aTTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGa                                                                                                                                                >  2:155218/1‑149 (MQ=255)
     aTTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAAc                                                                                                                                               >  1:165963/1‑150 (MQ=255)
       ttCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGAATTTTTCAGaaa                                                                                                                                              >  1:711089/1‑147 (MQ=255)
       ttCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACaa                                                                                                                                             >  2:462038/1‑150 (MQ=255)
         cGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAAt                                                                                                                                           <  2:1029001/150‑1 (MQ=255)
          gCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAGTTTTTcttgagaaa                                                                                                                                            >  1:159698/1‑140 (MQ=255)
              aaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCCCACAAAAACaaa                                                                                                                                      >  1:892606/1‑150 (MQ=255)
              aaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAAAAAAACaaa                                                                                                                                      >  2:1143129/1‑150 (MQ=255)
                    gCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTgg                                                                                                                                <  2:247814/150‑1 (MQ=255)
                     ccATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCACACAaaaaaaaaaa                                                                                                                                    >  1:1106912/1‑145 (MQ=255)
                       aTTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGtt                                                                                                                            <  1:374861/151‑1 (MQ=255)
                       aTTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGtt                                                                                                                            <  2:976917/151‑1 (MQ=255)
                       aTTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAAGTTTTTTAAGCAAAAAAAAAAAATGGAggt                                                                                                                            >  1:365513/1‑149 (MQ=255)
                          cccTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTTCGACAaaaaaaaaaaatagggtt                                                                                                                            >  1:1088265/1‑140 (MQ=255)
                              gCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCatg                                                                                                                      <  1:307339/150‑1 (MQ=255)
                                 ctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCATACAAAAAAAAAAAATGGAGtttgcttg                                                                                                                      >  1:712180/1‑141 (MQ=255)
                                      gcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGa                                                                                                              <  1:1134865/150‑1 (MQ=255)
                                        aaTATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGAAAAATAAAAAAACTGGTGATGGCATGGTGAGCGaa                                                                                                             >  1:517878/1‑149 (MQ=255)
                                          tatTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAAAAAAAAAAAAGGAGTTGCCATGATGAGCGgaa                                                                                                            >  1:994294/1‑148 (MQ=255)
                                             tcattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaat                                                                                                       <  1:479231/150‑1 (MQ=255)
                                                    tATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAAc                                                                                                >  1:350765/1‑150 (MQ=255)
                                                    tATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCa                                                                                                  <  1:914317/148‑1 (MQ=255)
                                                     aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAg                                                                                              <  2:787486/151‑1 (MQ=255)
                                                        ttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGc                                                                                            <  2:32644/150‑1 (MQ=255)
                                                         tCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCa                                                                                           <  2:705642/150‑1 (MQ=255)
                                                          cgcgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTTAAAAAAAAAAAAAAAATTGAGGTTGCCTGATGGGGGAAAGGAAAAAACaaaa                                                                                               >  1:518148/1‑143 (MQ=255)
                                                            cgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGa                                                                                       <  2:1061253/151‑1 (MQ=255)
                                                             gTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGaa                                                                                      <  2:90803/151‑1 (MQ=255)
                                                                  ttCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAATGGAGTTGTCATGATGAGCGAAAGAAATAAACAAAAAGCAGGGaaaa                                                                                    >  1:63178/1‑146 (MQ=255)
                                                                   tCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAACGGAGTTGGCCTGATGAGGGGAAGTAAAAAAAAACAAGCAATGAATaa                                                                                   >  2:1110795/1‑148 (MQ=255)
                                                                     aCGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGa                                                                              <  1:369717/151‑1 (MQ=255)
                                                                       gAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAATATAAAAAATGGAGTTTGCACGATGAGCGAAAGGAATAAAAAACAAGAAGTGAAAAAGt                                                                                 >  1:1061524/1‑146 (MQ=255)
                                                                           tGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGAGGTGTCATGATGAGCGGAAAGAATAAACAAAAGGCAGGGAATaa                                                                                   >  2:649532/1‑140 (MQ=255)
                                                                           tGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAAAAAAAAATGGAGTTTGCATGATTAGGGGAAAGAATAAACAAAAAGCaaggaaaaaatt                                                                                >  1:239370/1‑132 (MQ=255)
                                                                           tGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAAAAAAAAAAAAAAAAGGGGTTGGCCTGATGAGGGAAAGAAATAAAAAACAAGCAGTGAATAAATTGaaa                                                                            >  1:1111945/1‑145 (MQ=255)
                                                                               aaTTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGt                                                                    <  1:274385/151‑1 (MQ=255)
                                                                                  tACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCg                                                                  <  2:1035058/150‑1 (MQ=255)
                                                                                          cACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATTGAGTTGCCCTGATGAGCGGAAAGAATAAGCAAAAAGCAGTGAATAAGGTGGAAGAGATTTTCGCTAACttt                                                          >  2:422868/1‑150 (MQ=255)
                                                                                          cACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAAAACAAAAAATGGAGTTTCCCCGATGAGGGGAAGGAATAAAAAACAAGCCAAGGAAAAAATGAAAGAGATTTTCGCTAAAttt                                                          >  1:771676/1‑150 (MQ=255)
                                                                                                       aaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGAAAAAAAAAAAAAACTGAGGTGTCCTGGTGGGGGGAAGGAATAAGCAACAAGCAGGGGATAAGGTGGCAGAGAATGTCGCTAACttttaccccccgattt                                             >  2:922843/1‑137 (MQ=255)
                                                                                                             aaaCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAATATAAAAAATTGAGGTGGCCTGATGAGGGGAAGGAAAAAAAAAAAAGCAGGGAAAAAAATTAAAGAGGATGTCGCCAAATTTTTCGCCAAGAtttttttcc                                       >  1:99868/1‑144 (MQ=255)
                                                                                                                  acaTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAAAACAAAAAATGGAGGTTGCATGATGAGCGGAAGTAATAAACAACAGGCAGGGAATAAATTGACAGAGATTGTCGCTAAATTTTCCGCCATGATTTCTACCCCaaag                                  >  2:63299/1‑148 (MQ=255)
                                                                                                                    aTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTTTGAAAAATTAAAAAACAGGAGGTGTCATGATGAGGGGAAAAAAAAAAAAAAAAGGAGGGAATAAAATTACAAAGAATGTTGCTAAATTTTCCGCCATGATTTCTAcccccaagc                                 >  1:687433/1‑143 (MQ=255)
                                                                                                                    aTTTTTAGTGGTAGTTTGTGGCGGCGAATTTTTTCAGACAAAAAAAAAAAACGGGGGTGCCATGATGGGCGAAAAAAAAAAGCAACACGCCGTGAATAAGATGAAAGTGACTATCCGTTACTTTAACCCCAGGAATTTTAcccacacgcc                                >  2:177926/1‑143 (MQ=255)
                                                                                                                        ttAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCAGGAGGAGCGAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtgat                            <  2:234453/150‑1 (MQ=255)
                                                                                                                           gTGGTATTTTGTGGCGGTGATTTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAATTTTACCGCCATGATTTCTACCCGAATGCCtgatg                           <  1:451864/148‑1 (MQ=255)
                                                                                                                            tGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACg                        <  1:580649/150‑1 (MQ=255)
                                                                                                                             ggTAGTTTGTGGCGGTGAATTTTTCCAGAAAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgt                       <  1:780434/150‑1 (MQ=255)
                                                                                                                                aGTTTGTGGCGGTGAATTTTTTCCAACAAAAACAAAAAATGGGGGTTCCCTGATGAGGGGAAGGAATAAAAAACAAGCAGGGAATAAATTGACAGAGATTGTCGCTAACTTTACCCCCATGATTTTTTCCCGAAAGCCTGATGACgggggg                   >  1:861181/1‑149 (MQ=255)
                                                                                                                                 gTTTGGGGGGGTGAATTTTTCAAGAAAAATAAAAAAACTGGAGTTGCCATGATGAGCGAAAGTAAAAAACAACAGGTTGTAAATAGGTTGACAGAGATTGTCGCTATTTTGACCGCCATGATTTCTACCCGAATGCCTGATGACgtggt                    <  1:377943/149‑1 (MQ=255)
                                                                                                                                  tttGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgtggt                    <  2:251670/148‑1 (MQ=255)
                                                                                                                                     gtgGCGGTGAATTTTTCCAGACAAATAAAAAAACTGGTGTTGCAAAGATGAGAGAAAGTAAAAAGCAACAGGCAGTAAATAAGTTGAAAGAGTTTGCCGCTATCTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATaa              <  1:552581/151‑1 (MQ=255)
                                                                                                                                       ggCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATaaa             <  2:763797/150‑1 (MQ=255)
                                                                                                                                       ggCGGTGAATTTTTCCAGACAAATAAAAAAACTGGAGTTGCCTTGATGAGCAAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGTTTGTCGCTAATTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAAc            <  1:462020/151‑1 (MQ=255)
                                                                                                                                                aTTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAAc      >  1:384924/1‑149 (MQ=255)
                                                                                                                                                  ttttCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCt  >  2:997037/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
AGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCT  >  NC_000913/3206286‑3206581

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