Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,879,062 A→T 15.7% M353K (ATG→AAG)  gyrB ← DNA gyrase, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,879,0620AT15.7% 169.1 / 5.6 74M353K (ATG→AAG) gyrBDNA gyrase, subunit B
Reads supporting (aligned to +/- strand):  ref base A (30/31);  new base T (2/10);  total (32/42)
Fisher's exact test for biased strand distribution p-value = 5.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAA‑GTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTT‑GCCAACCACGATTTTCGCGTCGGTTGGG‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTACCGCCGTCAC  >  NC_000913/3878840‑3879293
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                        
gAGCCCCCTCCCACCAGGTGTGATCCGGAAGGGTCGGGATCAGGTCCCGGGCTGTCTTCCAGTTTCCCCGGCGGCCCCTTTAG‑GCCGCGGGAACCTTTCCGGCGGGTATTTTCCCGCGCGCGACCCCCCTCTTCGGGGGCCCGGGCAGCATCGATAATTTT‑GCCAACCACGTTTTTCGCTTCGTTTGGG‑TTTGCCAG‑CAGATTTTCTGCCAGAAGTTCGTTCATCTGCTG‑TTCAAccg                                                                                                                                                                                                                          <  2:171776/238‑1 (MQ=255)
         tCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCCGTCTGCCAGTTTGCCCGGCAGGCCCGCTAA‑GTCGAGCGCACCTTTACGGCGGGGCATTTCATCCCCGCGCCGGCCCGCGCCCCGGGCCCGGGCAGTATAGCTGATTTTGGCCAAAC‑CTCGCTTCGCGTCGGTTGCG‑GTTTTTTG‑CAGGCACTATGCCACCGCTGCGTTCTTCTTCTG‑CTCTAccgccggcatca                                                                                                                                                                                                                 >  2:146631/1‑232 (MQ=255)
           cgccgggaaCAGTCGGGAACGCGCCGGACGGCGTCCCGGGCATCCTG‑CAGGTTGCCCGGGCCGCCCGCCAAGGGCGAGCGCCACTTGCCGGGGGGTCTTTTCCGGCGCGGGCCGCGCCGCTTCACGGGCACGGGCACCATCGATATTTTT‑GCCACCCCCGTTTTTCGCGCGGGTTGGG‑TTTTCAGG‑CAGTTATTCCGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCt                                                                                                                                                                                                              <  2:1367909/231‑1 (MQ=255)
             ccaGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAA‑GTCGAGCGCACCTTTACGGCGGGGCATTTTACCCGCGCGGCGCGCCCCCTCACGGGGACGGGGGGGACCCGGAAATTT‑TCCAACCACCATTTTCGCGGCGGTTGGG‑TGTTTTAG‑CCGGGATTCTGCCAGCCGTTCGTTTATCTTCTG‑TTCATccgcccccttctcccc                                                                                                                                                                                                             >  2:596705/1‑227 (MQ=255)
                  tACAGTTCGGAAAGCGCCGGATCGCGTTCCGGGCATTCTGCCAGTTGGCCCGGCAGGCCCGCTAA‑GTCGAGCGCACCTTTACGGGGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTT‑GCCAACCACGTTTTTCGCGTCGGTTGGG‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCATCTGTTG‑TTCAACCGCCGATTTCACCTCagaa                                                                                                                                                                                                         <  1:308699/237‑1 (MQ=255)
                   gcAGTTCGGGAACCGCGGGTGCGCGTTCCGGGCATTCGGCCAGTTTCCCCGGAGGGCCCGCTAA‑GTCGGGCGCTCCTTTACGGGGGGTCATTTCCCGCGCGCGACGCGCCGCTTCCCGGGCACGGGCAGCATCGATAATTTT‑GCCAACCACGTTTTTCGCGTGGGTTGGG‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCTTCTGCTG‑TTCAACCGCCGATTTCACCTCAgaagaa                                                                                                                                                                                                      <  2:11523/238‑1 (MQ=255)
                     gcttggaaacgcggcgGTTCGCGTCCCGGGCATTCTCCGATTTGCCCCGCCAGGCCCCTAAA‑GTCGGGCGCCCCTGTGGGGCGGGTCATTCCCCGCGCGCGCCGCGCCGCTTCACGGGCACGGGCGGCATCGATAATTTT‑CCCAACCACGTTTTTCGTGTCGGTTGGG‑TTTTCCAG‑CGGGTTTTCTCCAAGCAGTTCGTTCATCTGTTG‑TTAAACCGCCGATTTCACCTCAGAAGAAAc                                                                                                                                                                                                    <  2:636172/229‑1 (MQ=255)
                           gAACGCGCCGGACCGCGTCCCTGGCAGTCGGCCAGTTTGCCCGGCAGGCCCGCTAA‑GTCGAGCGCCCTTTTACGGGGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTT‑CCCAACCACGATTTCCGCGTCGGTTGGG‑TTTTCCAG‑CAGGTTTTCTGCCAGCAGTTCGTTCTTCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGtt                                                                                                                                                                                               <  2:396778/238‑1 (MQ=255)
                           gAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAA‑GTCGAGCGCACCTTTACGGCGGGTCATTTTTCGCGCGCGCCGCGCCGCCTCCCGGCCCCGGGCAGCATCGGCAATTTT‑TCCACACCCGACTTTTTCGTCGGGGTGG‑GTTTCCCG‑CCGGGATTTTGTCTCCAGCTACTTTTTTTTCTG‑TTCCAACGCCGCCTTCTCCTCCCTAGAAAACAGttt                                                                                                                                                                                              >  2:1336818/1‑239 (MQ=255)
                                        cgcgTTCCGGGCAGTCTGCGAGTTCGCCCGGCGGGCCGGCTAA‑GCGGAGCGCCCCTTTCCGGCGGGTCTTTTCACGGGCGCGCCGCGCCGCTTCCCGGGCGCGGGGAGCATCGATAATTTT‑GCCAACCACGATTTTCGGGTGGGTT‑GG‑TTTTCCAG‑CAGTTATTCTGCCAGCATTTCTTTCTTCTGTTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTgg                                                                                                                                                                                 <  2:625087/238‑1 (MQ=255)
                                          cgTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAA‑GTCGAGCGCACCTTTACGGCGGGTCACTTTTCGCGCGCGCCGCCCCGCTTCACCGGGACGGGGGGCAACGATAATTTT‑GCCAACCACCACTTTTGCGGCGGTTGGG‑GTTTTCCG‑CAGGGGTTTTTCCCCACGGTCTTTCATCTGCTG‑TTCAAACCCCCCGTTTTCCTCACAAGAAACCCACTTGTCttttttttgggt                                                                                                                                                                               >  2:939279/1‑230 (MQ=255)
                                            tcccgggCAGTTCGCCAGTTTGCCCGGCCGGCCCGCGAA‑GCCGAGCGCCCCTTTCGGGGGGGTATTTTCACGCGCGCGCCGCGCGGTTCCCCGGGAGCGGGCAGCATCGATAATTTT‑GCAAACCACGTTTTCCGGGTGGGTTGGT‑TTTTCCAG‑CAGGTATTCTCCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGggagg                                                                                                                                                                             <  2:628268/234‑1 (MQ=255)
                                                                              gCTAG‑GTCGAGCGCACCTTTACGGCGGGTCTTTTCCCGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTT‑GCCAACCACGATTTTCGCGTCGGTTGGG‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGAc                                                                                                                                           <  1:492065/239‑1 (MQ=255)
                                                                                      cgggccctcttttCCGGGGGGTCATTTCCCGCGCGCGCGGCGCCGCTGCCGGGGCACGGGCAGCCTCGATATTTTT‑GCCACCCCCTTTTTGCGTTTTGGTTGGT‑TTTTCCGG‑CAGGTATTCGCCCACCAGTTCGTTCATCTGTTT‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATc                                                                                                                                    <  2:628023/229‑1 (MQ=255)
                                                                                                 tgagggggggTCATTGCCCGCGCGCGCGGCGCCGCTTCACGGGCAGGGGCAGCATCGATAATTTT‑CCCCCCCACGTTTTCCGCGTCGGTTGGG‑TTTCCCAG‑CAGGTTTTCTGCCAGCAGTTCGTTCATCTGTTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCAcg                                                                                                                         <  2:900402/232‑1 (MQ=255)
                                                                                                     ggcggGTCATTCCCCGCGCGCGCCGCGCCGCTTCCCGGGCACGGGCAGCATCGATAATTTT‑CCCAACCACGATTTCCGCGTCGGTTGGG‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCa                                                                                                                     <  1:146091/239‑1 (MQ=255)
                                                                                                      gcggGTCTTTTCACCCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTT‑GCCAACCACGATTTTCGCGTCGGTTGGG‑TTTTCCAG‑CAGTTATTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCa                                                                                                                     <  1:1326172/238‑1 (MQ=255)
                                                                                                      gcggGTCATGTCAGGCGCCGGCCGCGCGGCTACAGGGGCACGGACAGCGTGAATATTTTT‑CCAACCCACGTTTTTCGCGTCGTTTGGT‑TTGTCCAG‑CAGGTTTCCGCCCAGCAGTTCGTTCATCTGTTT‑TCCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCa                                                                                                                     <  2:1297212/238‑1 (MQ=255)
                                                                                                                    gcgcgggCCGCCCCGCTTCACGGGCCGGGGCAGAGTCGTATATTTT‑CCCAACCACTATTTTCGCGCTGTTTGGG‑TTTCCAAG‑CAGGTATTCTGCAAGCATTTGTTTAATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGgcg                                                                                                      <  2:742935/239‑1 (MQ=255)
                                                                                                                     cgggcgACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTT‑GCCAACCACGATTTTCGCGTCGGGTGGG‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGCAACCAGTTTGTCTTTGGTCTGGGCGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGACTTCACGCGCATCCTAACCGgcggcg                                                                                                      >  2:872418/4‑238 (MQ=255)
                                                                                                                     cgcgcgACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTT‑GCCAACCACGATTTTCGCGTCGGTTGGG‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGATAATTTCGGGTCCGGCACTTTTACGGAAACGACCGCAATCAGGCCTTCACGCGCATAGTCACCGGTTgc                                                                                                       >  2:488126/1‑237 (MQ=255)
                                                                                                                       cgcgACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTT‑GCCAACCACGATTTTCGCGTCGGTTGGG‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGGGAATTTCGGGTCCGGCACTTTCACGGAAACGCCCGCAATCTGGCCTTCAGGCGCTTCGTCCCCGCTGGCCCtc                                                                                                   >  2:61391/1‑238 (MQ=255)
                                                                                                                         cgACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTT‑GCCAACCACGATTTTCGCGTCGGTTGGG‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTg                                                                                                   >  1:322136/1‑237 (MQ=255)
                                                                                                                                        cgggcacgggcaGATTGGATAATTTT‑CCAACCCACGATTTTCGCGTGGGTTGTG‑TTTCCCAG‑CGGGTTTTCGGCCAGCATTTCTTTCTTCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCt                                                                                  <  1:336593/239‑1 (MQ=255)
                                                                                                                                         gggcacgggaaGCATCGATAATTTT‑GCCACCCACGATTTTCGCGTCGGTGGGT‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCt                                                                                  <  1:1198522/238‑1 (MQ=255)
                                                                                                                                             acgggcaGCATCGATAATTTT‑GCCAACCACGATTTTCGCGTCGGTTGGG‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAg                                                                              >  2:422835/1‑238 (MQ=255)
                                                                                                                                             acgggcaGCATCGATAATTTT‑CCCAACCCCGATTTTCGCGTCGTTTGGG‑TTTTCCAG‑CAGTTTTTCGGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGc                                                                             <  1:440348/239‑1 (MQ=255)
                                                                                                                                             acggccaGCATCTATAATTTT‑GCCAACCACGATTTTCGCGTCGGTTGGG‑TTTTCCAG‑CAGGTTTTCTGCCAGCATTTCGTTCATCTGCTG‑TTAACCCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGc                                                                             <  2:1096902/239‑1 (MQ=255)
                                                                                                                                              cgggcaGCATCGATAATTTT‑GCCAACCACGATTTTCGCGTCGGTTGGG‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAg                                                                              >  1:783619/1‑237 (MQ=255)
                                                                                                                                              cgggcaGCATCGATAATTTT‑GCCAACCACGATTTTCGCGTCGGTTGGG‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAg                                                                              >  2:68176/1‑237 (MQ=255)
                                                                                                                                              cgggcaGCATCGATAATTTT‑GCCAACCACGATTTTCGCGTCGGTTGGG‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAg                                                                              >  1:1064874/1‑237 (MQ=255)
                                                                                                                                                  gagtatgaaTAATGTT‑GCTATCCACGATTTTCGGTGGGGTTGGG‑TTTTCCAT‑GGGTGATTCCGTGCACAGTTTGTACATCGGATG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCcttc                                                                         <  1:445884/230‑1 (MQ=255)
                                                                                                                                                   agcaTCGATAATTTT‑GCCAACCACGATTTTCGCGTCGGTTGGG‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCcttctt                                                                       >  2:1214009/1‑239 (MQ=255)
                                                                                                                                                     caTCGATAATTTT‑GCCAACCACGATTTTCGCGTCGGTTGGG‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTTGCCTTCTTTg                                                                     >  2:245373/1‑239 (MQ=255)
                                                                                                                                                     caTCGATAATTTT‑GCCAACCACGATTTTCGCGTCGGTTGGG‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTg                                                                     >  2:893336/1‑239 (MQ=255)
                                                                                                                                                     caTCGATAATTTT‑GCCAACCACGATTTTCGCGTCGGTTGGG‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTg                                                                     >  2:508849/1‑239 (MQ=255)
                                                                                                                                                      aTCGATAATTTT‑GCCAACCGCGATTTTCGCGTCGTTTGGT‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCATCGGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCttt                                                                      <  1:1359505/237‑1 (MQ=255)
                                                                                                                                                      aTCGATAATTTG‑GCCACCCACGATTTGCGCGTCGGTTGGT‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCTTCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGt                                                                    <  1:532057/239‑1 (MQ=255)
                                                                                                                                                       tagATAATTTT‑CCCACCCACGATTTTCGCGCCGGTTGTG‑TTTCCCAG‑CGGGTATTCTGCCAGCATTTCGTTCATCGGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCcttct                                                                        <  2:946703/232‑1 (MQ=255)
                                                                                                                                                          aTAATTTT‑GCCAACCACGCTTTTCGCGTCGGTTGGG‑TTTTCCAG‑CAGGTATTCTTACCGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGCTTTCACCTCAGAAGGAACCAGTTTGTCTTTGGTCTGGGAGCGGCATTTCGGGTCCGGCACTTTCACGGCCTCGACCGCAATCCGGCATTCACGCGCATCGTCACCGGGGGCGTTGACTTTCGCTTTTTTGCTGCAGCATTCTTTGTCCAt                                                                >  2:747158/1‑239 (MQ=255)
                                                                                                                                                              tttt‑GCCAACCACTATTTGCGGGTGGTTTGGG‑TTTTCCGG‑AAGGTATTCTGCCAGCAGTTCGTGCATCTGTTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAg                                                            <  1:740951/239‑1 (MQ=255)
                                                                                                                                                               tttgtccAACGACTGTTTTCGCGTCGGTTGTG‑TTTGCCAG‑AAGGTATTCTCCCAGCATTTCTTTCAGCTGTTG‑TTAAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAg                                                            <  1:488126/234‑1 (MQ=255)
                                                                                                                                                                tt‑GCCAACCACATTTTTCG‑GGCCGTTGGG‑TTTTTCAGCCTGGTTTTCTGCCA‑CGGTTTGTTTTACCTCTGTTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGc                                                          <  2:488704/239‑1 (MQ=255)
                                                                                                                                                                 t‑GCCAACCACGATTTTCGCGTCGGTTGGG‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGt                                                              >  1:946703/1‑234 (MQ=255)
                                                                                                                                                                    ccAACCCCGTTTTCCGCGTGGTTTGGG‑TTTTCCAG‑CAGGTATGCGGCGAGCGGTTCGTTCATCTGCTT‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGtt                                                       <  2:410838/239‑1 (MQ=255)
                                                                                                                                                                    ccAACCACGATTTTCGCGTCGGTTGGG‑TTTTCCAG‑CCGGTATTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGCAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTATTTGTCCATGTAGGCGt                                                        >  1:1210082/1‑238 (MQ=255)
                                                                                                                                                                    ccAACCACGATTTTCGCGTCGGTTGGG‑TTTTCCAG‑CAGGTATCCTGCCAGCAGTTCGTTCATCTGTTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGtt                                                       <  1:928743/239‑1 (MQ=255)
                                                                                                                                                                    caaaCCACGATTTTCGCGTCGGTTGGG‑TTTTCCAG‑CAGTTATTCTGCCAGCAGTTCTTTCATCTGCTT‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGtt                                                       <  1:122586/237‑1 (MQ=255)
                                                                                                                                                                     cccccccgtttttTGCGGCGGGTGGGTTTTTCCAG‑CAGGTTTTCCGCCGGCTGTTCCTTCCTTTGCTT‑TCCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGtt                                                       <  2:1228690/230‑1 (MQ=255)
                                                                                                                                                                     cAACCACTATTTTCGCGTGGGTTGGG‑TTTTCCAG‑CAGGTTTTCGGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGtt                                                       <  1:960852/238‑1 (MQ=255)
                                                                                                                                                                      acccACGTTTTGCCCGGCGGTTGGG‑TTTTCCAG‑CAGTTTTCCACCAAGATGTTCGTTAATCGGCTT‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCa                                                     <  1:331492/237‑1 (MQ=255)
                                                                                                                                                                                   gcggcgGTGGGT‑TTTCCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCATCGCCTG‑TTCAGCCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTGTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGGCCATGTAGGCGTTCAGGGTACGGGTCa                                         <  1:1116401/238‑1 (MQ=255)
                                                                                                                                                                                   cgttggtttGTG‑TTTTCCAG‑CAGTTATCCCGCCAGGAGTCCGTCCATCTGGTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCAt                                        <  1:376224/232‑1 (MQ=255)
                                                                                                                                                                                           tGGG‑TTTTCCAG‑CAGGTTTTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCa                                 <  1:596705/238‑1 (MQ=255)
                                                                                                                                                                                             gg‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACgg                              >  1:1313841/1‑239 (MQ=255)
                                                                                                                                                                                                ttttCCAG‑CGGTTATTCTGCCAGCTGTTCGTTTATCTGCTG‑TACAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGa                             <  1:897731/238‑1 (MQ=255)
                                                                                                                                                                                                tttcCCAG‑CAGTTTTTCGGCCGGTAGTTCGTACATTGGCTG‑ATCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGa                             <  1:1064205/238‑1 (MQ=255)
                                                                                                                                                                                                 ttcccAG‑CGGGTTTTCTCCCCGCAGTTCTTTCTTCTCTTG‑TTCAACCGCCGATTTCCCCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAg                           <  1:188007/236‑1 (MQ=255)
                                                                                                                                                                                                  ttcccagaggttATCCTGCCGGCAGTGCGTTCTTCTGCGG‑TTAACCCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAg                           <  1:1357132/228‑1 (MQ=255)
                                                                                                                                                                                                          gggTTTTCTGCCAGCCGTTCTTTCTTCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGcc                     <  1:1244337/236‑1 (MQ=255)
                                                                                                                                                                                                           ggTATTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAgg                  >  2:842048/1‑239 (MQ=255)
                                                                                                                                                                                                            gTATTCTGCCAGCAGTTCTTTCATCTGCTT‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAgg                  <  1:1057828/238‑1 (MQ=255)
                                                                                                                                                                                                              aTTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGa               >  2:971927/1‑239 (MQ=255)
                                                                                                                                                                                                              aTTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGa               >  2:127132/1‑239 (MQ=255)
                                                                                                                                                                                                                  tGCCAGAATTTCGTTCTTCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAc           <  1:1185375/239‑1 (MQ=255)
                                                                                                                                                                                                                    ccAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGCAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAcc          >  1:1303228/1‑238 (MQ=255)
                                                                                                                                                                                                                      agcagTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAccgcc       >  2:426679/1‑239 (MQ=255)
                                                                                                                                                                                                                      agcagTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGGGAGTAccgc        >  2:791272/1‑238 (MQ=255)
                                                                                                                                                                                                                      agcagTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGCAGCCTTCTTTGTCCATGTAGGCGTTCATGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAcagc        >  2:1357594/1‑238 (MQ=255)
                                                                                                                                                                                                                       gcagTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAccgccg      >  1:784689/1‑239 (MQ=255)
                                                                                                                                                                                                                       gcagTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAccgccg      >  1:77958/1‑239 (MQ=255)
                                                                                                                                                                                                                       gcagTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAccgccg      >  1:206676/1‑239 (MQ=255)
                                                                                                                                                                                                                        cagTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAccgccg      >  1:536294/1‑238 (MQ=255)
                                                                                                                                                                                                                           ttcgttcATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCGGCCAGGTGAGTACCGCCGTccc  >  1:962044/1‑237 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                        
GAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAA‑GTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTT‑GCCAACCACGATTTTCGCGTCGGTTGGG‑TTTTCCAG‑CAGGTATTCTGCCAGCAGTTCGTTCATCTGCTG‑TTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTACCGCCGTCAC  >  NC_000913/3878840‑3879293

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑