Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,426,178 A→T 8.5% intergenic (+96/‑276) ydaY → / → tmpR pseudogene, Rac prophage;Phage or Prophage Related/Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,426,1780AT 8.5% 634.1 / 39.9 190intergenic (+96/‑276)ydaY/tmpRpseudogene, Rac prophage;Phage or Prophage Related/Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein
Reads supporting (aligned to +/- strand):  ref base A (82/91);  new base T (7/9);  total (90/100)
Fisher's exact test for biased strand distribution p-value = 8.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-02

Only 100 of 190 total aligned reads displayed.

TTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAAT  >  NC_000913/1425913‑1426451
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                 
                                    agcagAAGCAGGATACGATAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACGCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTCCCTTTGCCTTTATCCTTGACAGAGATAACTGAATATTTAGCCACTAATCCGGTattgatttaac                                                                                                                                                                                                                                      >  2:285193/1‑271 (MQ=255)
                                        gAAGCAGGATACGATAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCAc                                                                                                                                                                                                                                                          <  2:147497/251‑1 (MQ=255)
                                                aTACGATAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATAtt                                                                                                                                                                                                                                                                 >  2:636337/1‑236 (MQ=255)
                                                     aTAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTTGCTCTGACAGGTCTAAATCAATATTTAGCACCTCATCAGGTATTGATTGCGCTCGATGAATTTGAACCa                                                                                                                                                                                                                     >  4:528099/1‑275 (MQ=255)
                                                      tAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAg                                                                                                                                                                                                                    >  2:310931/1‑275 (MQ=255)
                                                                          aggggttgtTATCGGGGTGTGGGGGAGAAAGGGGGGAAAAGGAAAGAAGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGAc                                                                                                                                                                                                <  4:147833/272‑1 (MQ=255)
                                                                                      cGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGGGATCttt                                                                                                                                                                                                            >  2:590297/1‑251 (MQ=255)
                                                                                      cGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGGGATCGTTGCACTGGATGACCAGTATTTTCAg                                                                                                                                                                                    >  2:99495/1‑275 (MQ=255)
                                                                                           ttCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAg                                                                                                                                                                                    >  2:175149/1‑270 (MQ=255)
                                                                                           tgcGGGGGGGAAAGGGGGGAAAAGGCAAGAGGGAAGCGAGAAAAACTCAGGTGGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGc                                                                                                                                                                                <  4:634769/272‑1 (MQ=255)
                                                                                           gtCGGGGGGGAAAGGGGGGAAAAGGGAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGc                                                                                                                                                                                <  4:175148/273‑1 (MQ=255)
                                                                                                             gAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATAtt                                                                                                                                                                                                                                                                 >  4:299236/1‑175 (MQ=255)
                                                                                                                    cAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTa                                                                                                                                                       <  2:672396/274‑1 (MQ=255)
                                                                                                                        acagaggtgcgagaaaaTCCGGGTGCGACCAGTAACGGAAGCAGAAAAAGACCCGGTGGGTAATGGGGGTGTGTAGACGCATTCAGTAATATCCCATGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTGTGCCTTTATCGCTGACAGAGATAAATGAATATGTAGGTACTCATCCGGGATTGATTTAACGCGATGAATTTGAAGCAGTGATCGTTGCACTGGATGACCAGTATTTTCAGGAGCAGGGTGTGTAGTTGTTAATTACGTACAct                                                                                                                                                   <  4:430418/261‑1 (MQ=255)
                                                                                                                             aGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATa                                                                                                                                                                                                                                                                           >  7:715689/1‑149 (MQ=255)
                                                                                                                               cGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaa                                                                                                                                                                                                                                                                         <  8:124799/149‑1 (MQ=255)
                                                                                                                                gagaAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaag                                                                                                                                                                                                                                                                        <  7:354183/149‑2 (MQ=255)
                                                                                                                                gagaAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATaaa                                                                                                                                                                                                                                                                         <  6:3389/148‑1 (MQ=255)
                                                                                                                                   aaaaaCTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGa                                                                                                                                                                                                                                                                      >  6:151429/1‑148 (MQ=255)
                                                                                                                                    aaaaCTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGaa                                                                                                                                                                                                                                                                     >  6:439202/1‑148 (MQ=255)
                                                                                                                                        cTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGCCAGAGATAAATGAATAttt                                                                                                                                                                                                                                                                >  8:616241/1‑149 (MQ=255)
                                                                                                                                             gTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGcc                                                                                                                                                                                                                                                            <  6:85810/148‑1 (MQ=255)
                                                                                                                                                gCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGGGTTGTTGTACGGCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTTGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTc                                                                                                                           <  4:461085/274‑1 (MQ=255)
                                                                                                                                                 ccACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCa                                                                                                                                                                                                                                                       <  7:5557/149‑1 (MQ=255)
                                                                                                                                                  cACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATAttt                                                                                                                                                                                                                                                                <  8:678968/139‑1 (MQ=255)
                                                                                                                                                  cACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATAttt                                                                                                                                                                                                                                                                >  7:678968/1‑139 (MQ=255)
                                                                                                                                                   aCCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCAt                                                                                                                                                                                                                                                      <  5:500270/148‑1 (MQ=255)
                                                                                                                                                    ccGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATGCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATc                                                                                                                                                                                                                                                     >  6:641509/1‑148 (MQ=255)
                                                                                                                                                          ccGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTat                                                                                                                                                                                                                                               <  5:4026/148‑1 (MQ=255)
                                                                                                                                                            ggAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCAc                                                                                                                                                                                                                                                          <  7:240299/135‑1 (MQ=255)
                                                                                                                                                                ccAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgattg                                                                                                                                                                                                                                         <  5:151639/148‑1 (MQ=255)
                                                                                                                                                                  aGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAg                                                                                                                                                                                                                                       <  5:119028/148‑1 (MQ=255)
                                                                                                                                                                   gAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAgc                                                                                                                                                                                                                                      >  5:345481/1‑148 (MQ=255)
                                                                                                                                                                      aTAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGcgcg                                                                                                                                                                                                                                   <  5:641509/148‑1 (MQ=255)
                                                                                                                                                                          aCCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGa                                                                                                                                                                                                                               >  6:49398/1‑148 (MQ=255)
                                                                                                                                                                           cccGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGaa                                                                                                                                                                                                                              >  6:461911/1‑148 (MQ=255)
                                                                                                                                                                             cGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAAttt                                                                                                                                                                                                                           <  8:96866/149‑1 (MQ=255)
                                                                                                                                                                                tGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGaa                                                                                                                                                                                                                        <  7:465616/149‑1 (MQ=255)
                                                                                                                                                                                tGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGa                                                                                                                                                                                                                         <  5:599283/148‑1 (MQ=255)
                                                                                                                                                                                    tAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAg                                                                                                                                                                                                                    <  7:881258/149‑1 (MQ=255)
                                                                                                                                                                                           gttgttGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCCGTGATCTTGGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTCGTTTTTACTTACGTACCATCTGTTACAGAGGTGTGATGGTGTCTTTAGTTAAATCGATGATTCTACTTCAGAAAAGCATTGCGTgg                                                                                >  4:26205/1‑275 (MQ=255)
                                                                                                                                                                                           gttgttGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTACTGGAGAAAAGCTTTGCGTgg                                                                                >  2:392489/1‑275 (MQ=255)
                                                                                                                                                                                           gttgctGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCt                                                                                                                                                                                                              >  7:102544/1‑149 (MQ=255)
                                                                                                                                                                                            ttgttgTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAgc                                                                                                                                                                                                                                      >  8:799186/1‑124 (MQ=255)
                                                                                                                                                                                            ttgttgTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAgc                                                                                                                                                                                                                                      <  7:799186/124‑1 (MQ=255)
                                                                                                                                                                                              gttgtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGCGTGTGTAGTTGTTAATTa                                                                                                                                                           >  2:645849/1‑197 (MQ=255)
                                                                                                                                                                                                   aCGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTg                                                                                                                                                                                                      >  8:379382/1‑149 (MQ=255)
                                                                                                                                                                                                      ccTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGa                                                                                                                                                                                                    <  5:251300/148‑1 (MQ=255)
                                                                                                                                                                                                      ccTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGa                                                                                                                                                                                                    >  5:791971/1‑148 (MQ=255)
                                                                                                                                                                                                      ccTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGa                                                                                                                                                                                                    >  5:288390/1‑148 (MQ=255)
                                                                                                                                                                                                      ccTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCAGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGAt                                                                                                                                                                                                   <  8:566895/149‑1 (MQ=255)
                                                                                                                                                                                                         aTTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGNATTTNCAGGAGCAGTGTGTGTAGTTGTTAATTaag                                                                                                                                                         >  4:743387/1‑186 (MQ=255)
                                                                                                                                                                                                          ttCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGAc                                                                                                                                                                                                <  6:254727/148‑1 (MQ=255)
                                                                                                                                                                                                           tCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGATTATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGAcc                                                                                                                                                                                               >  6:5721/1‑148 (MQ=255)
                                                                                                                                                                                                            cGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCa                                                                                                                                                                                              <  5:739408/148‑1 (MQ=255)
                                                                                                                                                                                                             ggTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTTGGGTGCCTTTACCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGGAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGGGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATc                                                              <  1:207196/275‑1 (MQ=255)
                                                                                                                                                                                                              gTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGt                                                                                                                                                                                            >  6:234777/1‑148 (MQ=255)
                                                                                                                                                                                                              gTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGc                                                                                                                                                                                <  4:805556/160‑1 (MQ=255)
                                                                                                                                                                                                              gTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATc                                                              <  2:528101/274‑1 (MQ=255)
                                                                                                                                                                                                                aaTATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTAtt                                                                                                                                                                                         <  8:600328/149‑1 (MQ=255)
                                                                                                                                                                                                                 atatCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTGCAAGCAGTGATCTTTGCACTGGATGACCGGTATTTTCAGGAGCAGTGTATGTAGTTGTTAATTACGTACACTATGTTACAGAGATGTGAGGGTCTCTTTAGTTAACTCGATGCGGCTCCTGTAGAAAAGCATTGCCTGCCCTCGTACTCGCTACATctcct                                                          >  4:493325/1‑275 (MQ=255)
                                                                                                                                                                                                                     cccGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCa                                                                                                                                                                                     >  6:7956/1‑148 (MQ=255)
                                                                                                                                                                                                                     cccGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAg                                                                                                                                                                                    <  8:247397/149‑1 (MQ=255)
                                                                                                                                                                                                                           cccGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTgg                                                                                                                                                                                                     <  4:359460/126‑1 (MQ=255)
                                                                                                                                                                                                                            ccGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAg                                                                                                                                                                              <  6:741779/148‑1 (MQ=255)
                                                                                                                                                                                                                             cGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATg                                                                                                                                                                                                                                                                       <  1:95625/58‑1 (MQ=255)
                                                                                                                                                                                                                              gACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTc                                                                                                                                                                                      <  8:364763/138‑1 (MQ=255)
                                                                                                                                                                                                                              gACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAgg                                                                                                                                                                                   <  7:612673/141‑1 (MQ=255)
                                                                                                                                                                                                                                cGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAgtgtg                                                                                                                                                                          <  5:234777/148‑1 (MQ=255)
                                                                                                                                                                                                                                 gTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGa                                                                                                                                                                                  >  8:231263/1‑139 (MQ=255)
                                                                                                                                                                                                                                     tGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACg                                                                                                                                                         >  4:128994/1‑160 (MQ=255)
                                                                                                                                                                                                                                      gCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAgtt                                                                                                                                                                   <  7:98595/149‑1 (MQ=255)
                                                                                                                                                                                                                                         ggAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTaca                                                                                                                                                     <  2:128995/160‑1 (MQ=255)
                                                                                                                                                                                                                                          gAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCTCGATGAATTTGAAGCAGTGATCTTTGCACTGTATGACCAGTATTTTCAg                                                                                                                                                                                    >  3:124090/1‑128 (MQ=255)
                                                                                                                                                                                                                                           aaTGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATcc                                                                                                                                                                                                                                                    >  4:601111/1‑63 (MQ=255)
                                                                                                                                                                                                                                            aTGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAAt                                                                                                                                                             >  8:499220/1‑149 (MQ=255)
                                                                                                                                                                                                                                              ggCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCTCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGc                                                                                                                                                                                <  1:124091/128‑1 (MQ=255)
                                                                                                                                                                                                                                               gCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAg                                                                                                                                                                                    >  2:461796/1‑123 (MQ=255)
                                                                                                                                                                                                                                                      tgcctttgccttTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTacac                                                                                                                                                    >  5:845406/1‑148 (MQ=255)
                                                                                                                                                                                                                                                        cctttgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTa                                                                                                                                                                                           >  1:67660/1‑107 (MQ=255)
                                                                                                                                                                                                                                                         ctttgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTg                                                                                                                                                <  7:379382/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           ttgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGt                                                                                                                                               >  5:447480/1‑148 (MQ=255)
                                                                                                                                                                                                                                                           ttgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGt                                                                                                                                               <  6:845406/148‑1 (MQ=255)
                                                                                                                                                                                                                                                            tgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTAtttt                                                                                                                                                                                       <  3:67659/107‑1 (MQ=255)
                                                                                                                                                                                                                                                              cctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACa                                                                                                                                           <  8:880247/149‑1 (MQ=255)
                                                                                                                                                                                                                                                              catttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGTGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCa                                                                                                                                                                                     <  6:264809/105‑1 (MQ=255)
                                                                                                                                                                                                                                                                tttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTa                                                                                                                                                                      >  5:850445/1‑120 (MQ=255)
                                                                                                                                                                                                                                                                tttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACaga                                                                                                                                         >  8:706722/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                tttATCTCTGTAAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACaga                                                                                                                                         >  7:583608/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                tttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTa                                                                                                                                             <  8:158328/145‑1 (MQ=255)
                                                                                                                                                                                                                                                                tttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTa                                                                                                                                             >  7:158328/1‑145 (MQ=255)
                                                                                                                                                                                                                                                                 ttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACaga                                                                                                                                         >  6:528909/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                   aTCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAgaga                                                                                                                                       <  6:21824/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                      tctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAgttgt                                                                                                                                                                 <  2:756646/119‑1 (MQ=255)
                                                                                                                                                                                                                                                                      tctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTg                                                                                                                                                <  8:503593/136‑1 (MQ=255)
                                                                                                                                                                                                                                                                       ctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGAcc                                                                                                                                                                                               >  8:427692/1‑88 (MQ=255)
                                                                                                                                                                                                                                                                       ctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGAcc                                                                                                                                                                                               <  7:427692/88‑1 (MQ=255)
                                                                                                                                                                                                                                                                       ctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGAtgtg                                                                                                                                   >  5:531154/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                         gACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGAt                                                                                                                                 <  5:256765/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                         gACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTTCGTGGCCTCGTAATCGCTATCTCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGCTTTTTTGTTGTATGGAGTTTtaat  >  4:599662/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                 
TTGCAGAAGCAGCCAGCATCGCCCAGGGTAAAGAGCAGCAGAAGCAGGATACGATAAAAAAGCCATAGCTGCCCAGCGGTGGTTATCGGAGTTCGGGGGAGAAAGGGGGGAAAAGGCAAGATGGAAGCGAGAAAAACTCAGGTTGCCACCGATACCGGAACCAGAAATAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAAT  >  NC_000913/1425913‑1426451

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