Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,891,506 C→G 6.1% intergenic (+23/‑109) tnaB → / → mdtL tryptophan transporter of low affinity/multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,891,5060CG 6.1% 756.8 / 29.5 218intergenic (+23/‑109)tnaB/mdtLtryptophan transporter of low affinity/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (114/90);  new base G (7/7);  total (121/97)
Fisher's exact test for biased strand distribution p-value = 7.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 221 total aligned reads displayed.

TCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTGCCGGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGGCAGCTGCG  >  NC_000913/3891249‑3891776
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                              
        tCCCGAACGGCTTTATTTACGGGATCGGCGGTGCCGGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCCGTGCTTGCAATATAAGCTCGCAAGAAGTTTCCCACTCAGACGTTCACGGTCTCGGGCGGCAATCTTATTCCGGCTACTGTCATTCTCTTTGGTATAACCGTCATTTTGTGCTGGTTCGGCACCGTCTTTCAAGTGTTACTCACTTTTGTCTTAATCTTTCTCGAAGCCACCCTATCCCTGActtctttaccctcccgtaa                                                                                                                                                                                                                                                        >  4:132922/1‑261 (MQ=255)
                            gggATCGGCGGTGCCGGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTTCACGGACTGGGGCGGCAATCTTATTCCGGCGCTTGTCATTCTCTTTGGTATAACCGTGATTTTTTGCTGGTTCGGCCACGTCTTTAACGTTTTCCCTAAATTTGGCTTAATCCTTCAAGAAGCCACCCTTTCTCTGGCTTCTTGCCTCTCAGGAAATCCCTTATGTCGAAATGGCa                                                                                                                                                                                                                                    >  4:250282/1‑274 (MQ=255)
                                                 tgCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTACATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTTCCTCTCAGGTATTCACTTATGTCCAAATGGCAACTAGCCCGATCCTCCTTcacc                                                                                                                                                                                                              >  2:173132/1‑275 (MQ=255)
                                                                       ttattCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtc                                                                                                                                                                                                                        <  4:93784/243‑1 (MQ=255)
                                                                              cGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTactg                                                                                                                                                                                  >  2:150060/1‑271 (MQ=255)
                                                                                                             tttCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATCGTCATTCTCTTTGGTATAACCGTGACTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTc                                                                                                                                                                                                                                                                                <  5:1187296/149‑1 (MQ=255)
                                                                                                               tCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGAAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGc                                                                                                                                                                                                                                                                              >  8:1253752/1‑149 (MQ=255)
                                                                                                                 ccAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTg                                                                                                                                                                                                                                                                            >  7:853827/1‑149 (MQ=255)
                                                                                                                  cAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTgg                                                                                                                                                                                                                                                                           <  5:862067/149‑1 (MQ=255)
                                                                                                                  cAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTgg                                                                                                                                                                                                                                                                           >  5:358650/1‑149 (MQ=255)
                                                                                                                    aTCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGct                                                                                                                                                                                                                                                                         >  7:1428015/1‑149 (MQ=255)
                                                                                                                     tCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGctt                                                                                                                                                                                                                                                                        <  5:1433792/149‑1 (MQ=255)
                                                                                                                      cAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATGCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttc                                                                                                                                                                                                                                                                       <  5:1124107/149‑1 (MQ=255)
                                                                                                                       aGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGcttct                                                                                                                                                                                                                                                                      >  7:809349/1‑149 (MQ=255)
                                                                                                                       aGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGcttct                                                                                                                                                                                                                                                                      >  8:55598/1‑149 (MQ=255)
                                                                                                                        gATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGcttctt                                                                                                                                                                                                                                                                     >  6:849018/1‑149 (MQ=255)
                                                                                                                        gATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGcttctt                                                                                                                                                                                                                                                                     >  8:621313/1‑149 (MQ=255)
                                                                                                                         aTGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTg                                                                                                                                                                                                                                                                    >  6:564979/1‑149 (MQ=255)
                                                                                                                            ttCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGctt                                                                                                                                                                                                                                                                        <  8:482154/142‑1 (MQ=255)
                                                                                                                            ttCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGc                                                                                                                                                                                                                                                                   >  7:785999/1‑147 (MQ=255)
                                                                                                                            ttCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCct                                                                                                                                                                                                                                                                 >  7:1424247/1‑149 (MQ=255)
                                                                                                                                   tCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTAc                                                                                                                            >  2:108005/1‑275 (MQ=255)
                                                                                                                                       gggCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATc                                                                                                                                                                                                                                                      <  6:493439/149‑1 (MQ=255)
                                                                                                                                       gggCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATc                                                                                                                                                                                                                                                      >  5:1341039/1‑149 (MQ=255)
                                                                                                                                        ggcggcAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCa                                                                                                                                                                                                                                                     <  5:1425559/149‑1 (MQ=255)
                                                                                                                                         gcggcAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCAc                                                                                                                                                                                                                                                    >  8:221283/1‑149 (MQ=255)
                                                                                                                                           ggcAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCACCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACtt                                                                                                                                                                                                                                                  <  8:1379535/149‑1 (MQ=255)
                                                                                                                                             cAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCCACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCCGGACGCTTTAGCCTATCTCCAGCTTTTAGACTTGATGTTTTGCCCTGGTTTTACttttctccgc                                                                                                                  >  4:135885/1‑268 (MQ=255)
                                                                                                                                              aaTCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTa                                                                                                                                                                                                                                                 >  8:75111/1‑147 (MQ=255)
                                                                                                                                                  ttATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCa                                                                                                                                                                                                                                           >  7:250323/1‑149 (MQ=255)
                                                                                                                                                  ttATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGa                                                                                                             <  4:108005/275‑1 (MQ=255)
                                                                                                                                                    aTTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCaaa                                                                                                                                                                                                                                         <  6:200425/149‑1 (MQ=255)
                                                                                                                                                        cGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACt                                                                                                                                                                                                                                                   >  5:1235101/1‑135 (MQ=255)
                                                                                                                                                        cGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACt                                                                                                                                                                                                                                                   <  6:1235101/135‑1 (MQ=255)
                                                                                                                                                        cGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGc                                                                                                                                                                                                                                     <  7:566658/149‑1 (MQ=255)
                                                                                                                                                        cGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGc                                                                                                                                                                                                                                     >  5:1377451/1‑149 (MQ=255)
                                                                                                                                                         ggCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCa                                                                                                                                                                                                                                    >  5:435332/1‑149 (MQ=255)
                                                                                                                                                                gTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGcc                                                                                                                                                                                                                             <  8:685007/149‑1 (MQ=255)
                                                                                                                                                                 tCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCt                                                                                                                                                                                                                            <  6:977726/149‑1 (MQ=255)
                                                                                                                                                                   aTTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGa                                                                                                                                                                                                                          >  5:237246/1‑149 (MQ=255)
                                                                                                                                                                   aTTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGa                                                                                                                                                                                                                          >  5:747727/1‑149 (MQ=255)
                                                                                                                                                                    ttCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAt                                                                                                                                                                                                                         <  8:1425652/149‑1 (MQ=255)
                                                                                                                                                                       tctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACGTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcct                                                                                                                                                                                                                      >  6:30744/1‑149 (MQ=255)
                                                                                                                                                                       tctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcct                                                                                                                                                                                                                      >  7:204462/1‑149 (MQ=255)
                                                                                                                                                                         tttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATgg                                                                                                                                                                                                                                      <  8:409906/131‑1 (MQ=255)
                                                                                                                                                                         tttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATgg                                                                                                                                                                                                                                      >  7:409906/1‑131 (MQ=255)
                                                                                                                                                                         tttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcctcc                                                                                                                                                                                                                    <  8:1137746/149‑1 (MQ=255)
                                                                                                                                                                            ggTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGa                                                                                                             <  2:317112/249‑1 (MQ=255)
                                                                                                                                                                               attaCCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCg                                                                                              <  4:3346/258‑1 (MQ=255)
                                                                                                                                                                               ataACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTcacc                                                                                                                                                                                                              >  7:1034901/1‑149 (MQ=255)
                                                                                                                                                                                   ccGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCCAATGGCAACTCGCCTGATCCTCTTTCACCACGTATGCTTTGCGTCCCCTTACTATCAGGACGCTTTAGCCCTTGTCCCGCTTTTTGATTTGTAGTTTTTCCCTGGGTTTTCTTTATc                                                                                                                      >  2:290261/1‑233 (MQ=255)
                                                                                                                                                                                            tgtgCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGc                                                                                                                                                                                                 <  5:1004551/149‑1 (MQ=255)
                                                                                                                                                                                             gtgCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTgaa                                                                                                                                                                                                                                                                  >  7:1160130/1‑81 (MQ=37)
                                                                                                                                                                                                  ggTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAGGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCAcc                                                                                                                                                                                           <  6:17195/149‑1 (MQ=255)
                                                                                                                                                                                                      cGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGc                                                                                                                                                                                                                              <  7:770686/110‑1 (MQ=255)
                                                                                                                                                                                                           aCGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCa                                                                                                                                                                                  <  5:251237/149‑1 (MQ=255)
                                                                                                                                                                                                           aCGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCa                                                                                                                                                                                  <  5:155768/149‑1 (MQ=255)
                                                                                                                                                                                                            cGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAg                                                                                                                                                                                 <  8:1162567/149‑1 (MQ=255)
                                                                                                                                                                                                              tCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGa                                                                                                                                                                               <  6:747727/149‑1 (MQ=255)
                                                                                                                                                                                                               cTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGAc                                                                                                                                                                              >  8:1341964/1‑149 (MQ=255)
                                                                                                                                                                                                               cTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGAc                                                                                                                                                                              <  8:672358/149‑1 (MQ=255)
                                                                                                                                                                                                                   aaCGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTAACTCGTTGGTTTACAGCGCATCGCCGCCGATCTCAATGCCATCGCAGCGCAGtt                                             >  1:1926/1‑274 (MQ=255)
                                                                                                                                                                                                                      gtgtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAg                                                                                                                                                                       <  6:142634/149‑1 (MQ=255)
                                                                                                                                                                                                                        gtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTGTAGcc                                                                                                                                                                     >  6:1146687/1‑149 (MQ=255)
                                                                                                                                                                                                                         ttACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCttttt                                                                                                                                                      >  4:198976/1‑163 (MQ=255)
                                                                                                                                                                                                                         ttACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGa                                                                                                         >  4:68107/1‑208 (MQ=255)
                                                                                                                                                                                                                            ccTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATg                                                                                                                                                                 >  7:8812/1‑149 (MQ=255)
                                                                                                                                                                                                                            ccTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATg                                                                                                                                                                 >  8:581420/1‑149 (MQ=255)
                                                                                                                                                                                                                             cTAAATTTGGCTAAATCCTTCAAGAAGCCGGCCATTCGCGGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGt                                  <  2:132922/275‑1 (MQ=255)
                                                                                                                                                                                                                             cTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATg                                                                                                     <  2:68107/208‑1 (MQ=255)
                                                                                                                                                                                                                              tAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTc                                                                                                                                                               >  5:669303/1‑149 (MQ=255)
                                                                                                                                                                                                                                      gCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCtttt                                                                                                                                                       <  5:894100/149‑1 (MQ=255)
                                                                                                                                                                                                                                       cTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCttttt                                                                                                                                                      <  6:602838/149‑1 (MQ=255)
                                                                                                                                                                                                                                        tAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTg                                                                                                                                                     >  6:396917/1‑149 (MQ=255)
                                                                                                                                                                                                                                          aaTCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGAt                                                                                                                                                   <  6:1341039/149‑1 (MQ=255)
                                                                                                                                                                                                                                           aTCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGAtt                                                                                                                                                  <  5:814071/149‑1 (MQ=255)
                                                                                                                                                                                                                                            tCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGAttt                                                                                                                                                 <  7:724992/149‑1 (MQ=255)
                                                                                                                                                                                                                                            tCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGAttt                                                                                                                                                 <  6:832165/149‑1 (MQ=255)
                                                                                                                                                                                                                                              cTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCACTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGt                                                                                                                                               >  6:1041649/1‑149 (MQ=255)
                                                                                                                                                                                                                                               ttCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAGATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTa                                                                                                                                              >  8:790472/1‑149 (MQ=255)
                                                                                                                                                                                                                                               ttCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCtttt                                                                                                                                                       <  7:1184550/140‑1 (MQ=255)
                                                                                                                                                                                                                                               ttCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTa                                                                                                                                              >  7:1274057/1‑149 (MQ=255)
                                                                                                                                                                                                                                               ttCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACGATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTa                                                                                                                                              <  6:1188303/149‑1 (MQ=255)
                                                                                                                                                                                                                                                  aagaagCCAGCGATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGtt                                                                                                                                           <  7:947320/149‑1 (MQ=255)
                                                                                                                                                                                                                                                  aagaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGtt                                                                                                                                           >  8:348876/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    gaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGtttt                                                                                                                                         >  7:322328/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    gaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGtttt                                                                                                                                         >  6:698888/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    gaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGtttt                                                                                                                                         <  8:783611/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     aagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCGTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTg                                                                                                                                        <  6:441288/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      agccagcgAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGGTCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGAAGCTTTAg                                                                                                                                                                       <  6:794208/117‑1 (MQ=255)
                                                                                                                                                                                                                                                        ccagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGccc                                                                                                                                     <  6:343252/149‑1 (MQ=255)
                                                                                                                                                                                                                                                         cagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCt                                                                                                                                    <  7:1375274/149‑1 (MQ=255)
                                                                                                                                                                                                                                                         cagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCt                                                                                                                                    >  5:44492/1‑149 (MQ=255)
                                                                                                                                                                                                                                                           gccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGt                                                                                                                                                                                               <  3:26271/88‑1 (MQ=255)
                                                                                                                                                                                                                                                            ccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGt                                                                                                                                 <  5:396917/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            ccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGtt                                                                                                                                <  8:1013089/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             caTTGGCTGGTTTCTTGCCTCTCTGGAAATCACTTATGCCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGGCAGCTGCg  <  4:301444/275‑1 (MQ=255)
                                                                                                                                                                                                                                                               ttCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGtttt                                                                                                                              >  8:896247/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 cGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTAc                                                                                                                            <  7:489145/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 cGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTAc                                                                                                                            <  8:204462/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                              
TCTGATCTTCCCGAACGGCTTTATTTACGGGATCGGCGGTGCCGGGCTGTGCGCCACCATCTGGGCGGTCATTATTCCCGCAGTGCTTGCAATCAAAGCTCGCAAGAAGTTTCCCAATCAGATGTTCACGGTCTGGGGCGGCAATCTTATTCCGGCGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACGCTTTAGCCCATGTCCCGCTTTTTGATTTGTAGTTTTGCCCTGGTTTTACTTTATCCCGCCGGGATTGATATGTACCTCGTTGGTTTACCGCGCATCGCCGCCGATCTCAATGCCAGCGAAGCGCAGTTGCATATTGCGTTCTCCGTATATCTGGCGGGGATGGCAGCTGCG  >  NC_000913/3891249‑3891776

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