Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,314,233 A→T 5.3% intergenic (+36/+111) yjdP → / ← phnP putative periplasmic protein/5‑phospho‑alpha‑D‑ribosyl 1,2‑cyclic phosphate phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,314,2330AT 5.3% 996.2 / 32.6 284intergenic (+36/+111)yjdP/phnPputative periplasmic protein/5‑phospho‑alpha‑D‑ribosyl 1,2‑cyclic phosphate phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (133/135);  new base T (8/7);  total (142/142)
Fisher's exact test for biased strand distribution p-value = 7.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.71e-02

Only 100 of 285 total aligned reads displayed.

GCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATC  >  NC_000913/4313966‑4314499
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                          
                   tataACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAgcg                                                                                                                                                                                                                                                     <  2:633685/272‑1 (MQ=255)
                     taACCCTGACCGTGGGCGAGAGTATGAAGATTCGCGCGACGACCGCCGGCAGCGTGAAGTGAGCGACGATTGGCGCCGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATAGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAAc                                                                                                                                                                                                                                                          <  2:215596/265‑1 (MQ=255)
                            gACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGACGAGCGGAGAAGGGAGGATGAGTATGGNGGAGGAGGTTGGGTATGGAGCGGGGATGGCGCCTACACCCAAAGGGGGAAGGGGGCAAACGAGAGCCAGATCCaaaaaca                                                                                                                                                                                                                                      >  3:205007/1‑273 (MQ=255)
                             aCCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGCGCGTCGGCAACGGGAGGATGAAGAGCGGACAATGGAAGATGAGTATGGGCGATGAGGTTGGGTAGTGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCaaaa                                                                                                                                                                                                                                         >  4:373822/1‑272 (MQ=255)
                             aCCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGGATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAg                                                                                                                                                                                                                                                       >  2:538250/1‑260 (MQ=255)
                                       agaCTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACaa                                                                                                                                                                                                                               >  2:691619/1‑274 (MQ=255)
                                           tttGAAGAGTGCCGCGAAGACGGCGGGCAGCGGGGAGTGGGGGACGATCGGCGCAGACAAGAGGACGACCGCCGCCGCCAGTTTGAAGGTCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCgg                                                                                                                                                                                                                           <  4:420677/272‑1 (MQ=255)
                                                       cgcgACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                                                                          >  1:168985/1‑279 (MQ=255)
                                                              acgGCCGGCAGCTAGAAGTGAACGACTGTCTGCGCCAAAAATATGCCGGCCGCCGCCGCCTGTTTGAAGATCAACGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                                                                       <  1:262387/275‑1 (MQ=255)
                                                                ggCTGGAAGGGGGAAGTGGGCGACGATCGGCGCAGAAAAAAGGACGACCGCCGCCGACAGTTTGAAGATCGCGGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                                                                       <  4:691614/273‑1 (MQ=255)
                                                                    ggCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                                                                                  >  4:249521/1‑274 (MQ=255)
                                                                                            ggCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAg                                                                                                                                                                                                                                                       >  2:659783/1‑197 (MQ=255)
                                                                                              cgcAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGcaca                                                                                                                                                                        <  2:274837/274‑1 (MQ=255)
                                                                                                          gacgacCGCCGCGGCCAGTGTGAAGATCGCCGAAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCtcatca                                                                                                                                                           <  2:592654/275‑1 (MQ=255)
                                                                                                             cgccgccgccgCCAGTTTGAAGAGCGCGGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGc                                                                                                                                                        <  4:132401/273‑1 (MQ=255)
                                                                                                               ccgccgccgCCAGTTTGAAGATCGCCTCAGGCAACTGTACTATCTCAAACGCCAGCTTGATAAGGAGCGTCGGCAACTGGAGGCTGACGAGCGGCACATGGAATGTTAGTGAGGCAGACACGTTTTGGTATGGACCAGGGATCACCCCTCACCCAACTCGGAAGTGCGAGAAAACGAGCGCAATCTTCAATAACAAGCACAAGCGGTCCCCCTGCCACACTGGCGACAGGTCTAGGGTGGCGGGGCCAGCGTGCAACTTCCCACtcaacaaacc                                                                                                                                                       >  4:592650/1‑270 (MQ=255)
                                                                                                               ccgccgccgCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGa                                                                                                                                                                                                                                                        <  1:93740/177‑1 (MQ=255)
                                                                                                                        ccAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAAt                                                                                                                                                                                                                                                                           <  7:122453/149‑1 (MQ=255)
                                                                                                                         cAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATg                                                                                                                                                                                                                                                                          <  6:734327/149‑1 (MQ=255)
                                                                                                                          aGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATgg                                                                                                                                                                                                                                                                         >  8:1337992/1‑149 (MQ=255)
                                                                                                                           gTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATggg                                                                                                                                                                                                                                                                        >  6:582123/1‑149 (MQ=255)
                                                                                                                           gTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATggg                                                                                                                                                                                                                                                                        <  5:1598011/149‑1 (MQ=255)
                                                                                                                            tttGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATgggg                                                                                                                                                                                                                                                                       <  5:593207/149‑1 (MQ=255)
                                                                                                                               gAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGgggaga                                                                                                                                                                                                                                                                    <  6:1175228/149‑1 (MQ=255)
                                                                                                                                aaGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCTCCCCATTGgggtgag                                                                                                                                                                                                                                                                   <  6:1642122/149‑1 (MQ=255)
                                                                                                                                aaGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGgggagag                                                                                                                                                                                                                                                                   <  8:618194/149‑1 (MQ=255)
                                                                                                                                aaGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGgggagag                                                                                                                                                                                                                                                                   <  7:723076/149‑1 (MQ=255)
                                                                                                                                 aGATAGCCGAAGGCAGCGGGACGATCGCCAGCGCCAGCTTGAGCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccatcccat                                                                                                                                    <  2:522514/275‑1 (MQ=255)
                                                                                                                                  gATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGgggagagg                                                                                                                                                                                                                                                                  <  8:1200951/148‑1 (MQ=255)
                                                                                                                                  gATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGgggagagg                                                                                                                                                                                                                                                                  >  7:1200951/1‑148 (MQ=255)
                                                                                                                                   aTCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGgggagaggga                                                                                                                                                                                                                                                                >  7:453773/1‑149 (MQ=255)
                                                                                                                                   aTCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGgggagaggga                                                                                                                                                                                                                                                                >  6:829373/1‑149 (MQ=255)
                                                                                                                                    tCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCTCCCCATTGgggtgagg                                                                                                                                                                                                                                                                  <  8:472166/146‑1 (MQ=255)
                                                                                                                                    tCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCTCCCCATTGgggtgagg                                                                                                                                                                                                                                                                  >  7:472166/1‑146 (MQ=255)
                                                                                                                                           aggcagCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGa                                                                                                                                                                                                                                                        <  5:1447573/149‑1 (MQ=255)
                                                                                                                                             gcagCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAgc                                                                                                                                                                                                                                                      <  8:898208/149‑1 (MQ=255)
                                                                                                                                                  tGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAat                                                                                                                                                                                                                                                 <  8:547040/149‑1 (MQ=255)
                                                                                                                                                    gACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAAt                                                                                                                                                                                                                                                                           >  8:1304952/1‑121 (MQ=255)
                                                                                                                                                    gACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGgggaga                                                                                                                                                                                                                                                                    <  8:1429623/128‑1 (MQ=255)
                                                                                                                                                       gATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCTCCCCATTGgggtgaggggc                                                                                                                                                                                                                                                               >  5:1016002/1‑128 (MQ=255)
                                                                                                                                                       gATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCa                                                                                                                                                                                                                                            >  7:278229/1‑149 (MQ=255)
                                                                                                                                                        aTCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAgc                                                                                                                                                                                                                                                      <  5:1305194/138‑1 (MQ=255)
                                                                                                                                                        aTCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAgc                                                                                                                                                                                                                                                      >  6:1305194/1‑138 (MQ=255)
                                                                                                                                                        aTCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCaa                                                                                                                                                                                                                                           <  7:133735/149‑1 (MQ=255)
                                                                                                                                                           gccagcgccagcTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCANTATTCAATATCGAGCACAACCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCCCGCACACTCCCCNCCATCCCGCCACCCCAATCCCCAACNCCCCAAACCCCACCCCacacacccaccccacc                                                                                                       >  3:177921/1‑266 (MQ=255)
                                                                                                                                                            ccagcgccagcTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAGTGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATc                                                                                                                                                                                                                                       >  7:143809/1‑149 (MQ=255)
                                                                                                                                                            ccagcgccagcTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCg                                                                                                                                                                                                                                      >  8:1449205/1‑149 (MQ=255)
                                                                                                                                                              agcgccagcTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGAc                                                                                                             >  2:352466/1‑269 (MQ=255)
                                                                                                                                                                    agcTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACaa                                                                                                                                                                                                                               >  6:12671/1‑149 (MQ=255)
                                                                                                                                                                    agcTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACaa                                                                                                                                                                                                                               >  7:308251/1‑149 (MQ=255)
                                                                                                                                                                          aTCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTc                                                                                                                                                                                                                         >  7:760431/1‑149 (MQ=255)
                                                                                                                                                                                    tCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATggg                                                                                                                                                                                                                                                                        >  4:720256/1‑92 (MQ=255)
                                                                                                                                                                                     cGGCAACTGGAGGATGAAGAGCGGTGAATGGAAGATGAGTTTGTTCGATGAGGATGGGTATGGAGCAGGCATCGAACCTCACACCCATAGAGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACATTGGGTCCTTTCGCCCct                                                                                                                                                                                                              >  6:657268/1‑149 (MQ=255)
                                                                                                                                                                                         aaCTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                                                                          <  6:972337/149‑1 (MQ=255)
                                                                                                                                                                                          aCTGGAGGATGAAGAGCGGAGACTGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                                                                         <  8:1023647/149‑1 (MQ=255)
                                                                                                                                                                                            tGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                                                                       >  8:1058573/1‑149 (MQ=255)
                                                                                                                                                                                            tGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGAGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                                                                       >  8:1423184/1‑149 (MQ=255)
                                                                                                                                                                                               aggaTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                                                                                    <  5:829373/149‑1 (MQ=255)
                                                                                                                                                                                                 gaTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCTCCCCAGAGGGGCGAggg                                                                                                                                                                                                  <  8:760009/149‑1 (MQ=255)
                                                                                                                                                                                                  aTGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGggga                                                                                                                                                                                                                                                                      <  2:717094/80‑1 (MQ=255)
                                                                                                                                                                                                   tGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGtt                                                                                                                                                                                                <  6:488221/149‑1 (MQ=255)
                                                                                                                                                                                                   tGAAGAGCGGAGAATGGAAGATGAGAATGGGCGATGAGGTTGGGTATGAAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGTGTGGGtt                                                                                                                                                                                                >  6:175514/1‑149 (MQ=255)
                                                                                                                                                                                                    gAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGa                                                                                                                                                                                                                                                        >  5:95325/1‑92 (MQ=255)
                                                                                                                                                                                                    gAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAg                                                                                                                                                                                                                                                       >  7:1294043/1‑93 (MQ=255)
                                                                                                                                                                                                    gAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTg                                                                                                                                                                                                           >  8:1151166/1‑137 (MQ=255)
                                                                                                                                                                                                    gAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                                                                               <  5:885441/149‑1 (MQ=255)
                                                                                                                                                                                                      agagCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                                                                             >  6:1619947/1‑149 (MQ=255)
                                                                                                                                                                                                         gCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAgc                                                                                                                                                                                                                                                      >  8:1106775/1‑89 (MQ=255)
                                                                                                                                                                                                         gCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAgc                                                                                                                                                                                                                                                      <  7:1106775/89‑1 (MQ=255)
                                                                                                                                                                                                          cGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCa                                                                                                                                                                                                                                                   <  4:93258/91‑1 (MQ=255)
                                                                                                                                                                                                            gagaATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGt                                                                                                                                                                                                                          >  8:249617/1‑114 (MQ=255)
                                                                                                                                                                                                            gagaATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                                                                       <  8:278229/149‑1 (MQ=255)
                                                                                                                                                                                                                      gATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAg                                                                                                                                                                                              <  6:1381864/132‑1 (MQ=255)
                                                                                                                                                                                                                      gATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAg                                                                                                                                                                                              >  5:1381864/1‑132 (MQ=255)
                                                                                                                                                                                                                               gggCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTcc                                                                                                                                                                    >  5:1583925/1‑149 (MQ=255)
                                                                                                                                                                                                                               gggCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTcc                                                                                                                                                                    <  8:853001/149‑1 (MQ=255)
                                                                                                                                                                                                                                ggCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTccc                                                                                                                                                                   >  5:567945/1‑149 (MQ=255)
                                                                                                                                                                                                                                 gCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCct                                                                                                                                                                  <  7:19662/149‑1 (MQ=255)
                                                                                                                                                                                                                                 gCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAAc                                                                            >  4:653734/1‑235 (MQ=255)
                                                                                                                                                                                                                                  cGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCctc                                                                                                                                                                 <  6:1583925/149‑1 (MQ=255)
                                                                                                                                                                                                                                         gTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCCCAATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                                                                          <  7:1425024/101‑1 (MQ=255)
                                                                                                                                                                                                                                           tGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTg                                                                                                                                                                                                           >  2:76160/1‑98 (MQ=255)
                                                                                                                                                                                                                                                 tGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCct                                                                                                                                                                  >  7:1483495/1‑133 (MQ=255)
                                                                                                                                                                                                                                                   gAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                                                                                   <  7:440672/98‑1 (MQ=255)
                                                                                                                                                                                                                                                   gAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                                                                                   >  8:440672/1‑98 (MQ=255)
                                                                                                                                                                                                                                                    aGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAAt                                                                                                                                               <  5:1670196/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    aGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAAt                                                                                                                                               >  8:1554896/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           aTCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                                                                                    <  5:302798/89‑1 (MQ=255)
                                                                                                                                                                                                                                                           aTCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                                                                                    >  6:302798/1‑89 (MQ=255)
                                                                                                                                                                                                                                                              gCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccatccca                                                                                                                                     >  5:1105448/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               tcccTCTCCCCATTGGGGTGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTccatcccat                                                                                                                                    >  7:896704/2‑149 (MQ=255)
                                                                                                                                                                                                                                                                cccTCTCCCCATTGGGGTGAGGGGGAAAACGAGCGCAATATTCAATATCGCGCACAATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                                                                                     >  5:179697/1‑83 (MQ=37)
                                                                                                                                                                                                                                                                 ccTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                                                                          <  7:1131639/77‑1 (MQ=255)
                                                                                                                                                                                                                                                                  ctctCCCCATTGGGGTGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                                                                          >  5:1140190/1‑76 (MQ=38)
                                                                                                                                                                                                                                                                  cTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGGCGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCa               <  2:150865/263‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCATCCACGCATCAAAc                                                                            >  1:131348/1‑201 (MQ=255)
                                                                                                                                                                                                                                                                        cccAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAAcccc                                                                                                                            >  6:1011520/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         ccAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCa                                                                                                                           >  8:313675/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                          cAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCAc                                                                                                                          >  7:562467/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                          cAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCAc                                                                                                                          >  6:416799/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                          
GCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCGAGACTATGAAGATTCGCGCGACGACGGCTGGCAGCGTGAAGTGAGCGACGATCGGCGCAGACAATATGACGACCGCCGCCGCCAGTTTGAAGATCGCCGCAGGCAGCTGGACGATCGCCAGCGCCAGCTTGATCAGGAGCGTCGGCAACTGGAGGATGAAGAGCGGAGAATGGAAGATGAGTATGGGCGATGAGGTTGGGTATGGAGCAGGCATCGCCCCTCACCCCAATGGGGAGAGGGAGAAAACGAGCGCAATATTCAATATCGAGCACAATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGCCAGCGCGCACATTCCCTCTCATCACGCCACCCCAATCTCCATCCCATCAAACCCCACCTCAAACCCTGACGGTAGTGCGTTTTCCATCAGCCACGCATCAAACTGGTGGCTGATATGGGTCAGAATCACCCGTGGCGAGCGGATAACCTGATTCAGCGCAAGCACGGTATTTAAATC  >  NC_000913/4313966‑4314499

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