Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,879,061 A→T 15.2% 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,0610AT15.2% 173.8 / 5.3 76M353K (ATG→AAG) gyrBDNA gyrase, subunit B
Reads supporting (aligned to +/- strand):  ref base A (30/33);  new base T (2/10);  total (32/44)
Fisher's exact test for biased strand distribution p-value = 5.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTACCGCCGTCAC  >  NC_000913/3878839‑3879292
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                        
gAGCCCCCTCCCACCAGGTGTGATCCGGAAGGGTCGGGATCAGGTCCCGGGCTGTCTTCCAGTTTCCCCGGCGGCCCCTTTAGGCCGCGGGAACCTTTCCGGCGGGTATTTTCCCGCGCGCGACCCCCCTCTTCGGGGGCCCGGGCAGCATCGATAATTTTGCCAACCACGTTTTTCGCTTCGTTTGGGTTTGCCAGCAGATTTTCTGCCAGAAGTTCGTTCATCTGCTGTTCAAccg                                                                                                                                                                                                                          <  2:171776/238‑1 (MQ=255)
         tCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCCGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGGCATTTCATCCCCGCGCCGGCCCGCGCCCCGGGCCCGGGCAGTATAGCTGATTTTGGCCAAACCTCGCTTCGCGTCGGTTGCGGTTTTTTGCAGGCACTATGCCACCGCTGCGTTCTTCTTCTGCTCTAccgccggcatca                                                                                                                                                                                                                 >  2:146631/1‑232 (MQ=255)
           cgccgggaaCAGTCGGGAACGCGCCGGACGGCGTCCCGGGCATCCTGCAGGTTGCCCGGGCCGCCCGCCAAGGGCGAGCGCCACTTGCCGGGGGGTCTTTTCCGGCGCGGGCCGCGCCGCTTCACGGGCACGGGCACCATCGATATTTTTGCCACCCCCGTTTTTCGCGCGGGTTGGGTTTTCAGGCAGTTATTCCGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCt                                                                                                                                                                                                              <  2:1367909/231‑1 (MQ=255)
             ccaGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGGCATTTTACCCGCGCGGCGCGCCCCCTCACGGGGACGGGGGGGACCCGGAAATTTTCCAACCACCATTTTCGCGGCGGTTGGGTGTTTTAGCCGGGATTCTGCCAGCCGTTCGTTTATCTTCTGTTCATccgcccccttctcccc                                                                                                                                                                                                             >  2:596705/1‑227 (MQ=255)
                  tACAGTTCGGAAAGCGCCGGATCGCGTTCCGGGCATTCTGCCAGTTGGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGGGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGTTTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGTTGTTCAACCGCCGATTTCACCTCagaa                                                                                                                                                                                                         <  1:308699/237‑1 (MQ=255)
                   gcAGTTCGGGAACCGCGGGTGCGCGTTCCGGGCATTCGGCCAGTTTCCCCGGAGGGCCCGCTAAGTCGGGCGCTCCTTTACGGGGGGTCATTTCCCGCGCGCGACGCGCCGCTTCCCGGGCACGGGCAGCATCGATAATTTTGCCAACCACGTTTTTCGCGTGGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCTTCTGCTGTTCAACCGCCGATTTCACCTCAgaagaa                                                                                                                                                                                                      <  2:11523/238‑1 (MQ=255)
                     gcttggaaacgcggcgGTTCGCGTCCCGGGCATTCTCCGATTTGCCCCGCCAGGCCCCTAAAGTCGGGCGCCCCTGTGGGGCGGGTCATTCCCCGCGCGCGCCGCGCCGCTTCACGGGCACGGGCGGCATCGATAATTTTCCCAACCACGTTTTTCGTGTCGGTTGGGTTTTCCAGCGGGTTTTCTCCAAGCAGTTCGTTCATCTGTTGTTAAACCGCCGATTTCACCTCAGAAGAAAc                                                                                                                                                                                                    <  2:636172/229‑1 (MQ=255)
                           gAACGCGCCGGACCGCGTCCCTGGCAGTCGGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCCCTTTTACGGGGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTCCCAACCACGATTTCCGCGTCGGTTGGGTTTTCCAGCAGGTTTTCTGCCAGCAGTTCGTTCTTCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGtt                                                                                                                                                                                               <  2:396778/238‑1 (MQ=255)
                           gAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTTTCGCGCGCGCCGCGCCGCCTCCCGGCCCCGGGCAGCATCGGCAATTTTTCCACACCCGACTTTTTCGTCGGGGTGGGTTTCCCGCCGGGATTTTGTCTCCAGCTACTTTTTTTTCTGTTCCAACGCCGCCTTCTCCTCCCTAGAAAACAGttt                                                                                                                                                                                              >  2:1336818/1‑239 (MQ=255)
                                        cgcgTTCCGGGCAGTCTGCGAGTTCGCCCGGCGGGCCGGCTAAGCGGAGCGCCCCTTTCCGGCGGGTCTTTTCACGGGCGCGCCGCGCCGCTTCCCGGGCGCGGGGAGCATCGATAATTTTGCCAACCACGATTTTCGGGTGGGTTGGTTTTCCAGCAGTTATTCTGCCAGCATTTCTTTCTTCTGTTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTgg                                                                                                                                                                                 <  2:625087/238‑1 (MQ=255)
                                          gggtccttgcgagccgtccggttGCCCGGGAGGGCCCGGTAAGCCGAGCGCCCCTTTACGGGGGTTAATTTCACGCGCGCGCCGGCCCGCTTCCCGGGCAGGGGCAGCTTCGATAATTTTGCCACCCCGGTTTTTCGCGGCGTTTGGTTTTTCCAGCGGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTgg                                                                                                                                                                                 <  2:1277686/217‑1 (MQ=255)
                                          cgTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCACTTTTCGCGCGCGCCGCCCCGCTTCACCGGGACGGGGGGCAACGATAATTTTGCCAACCACCACTTTTGCGGCGGTTGGGGTTTTCCGCAGGGGTTTTTCCCCACGGTCTTTCATCTGCTGTTCAAACCCCCCGTTTTCCTCACAAGAAACCCACTTGTCttttttttgggt                                                                                                                                                                               >  2:939279/1‑230 (MQ=255)
                                            tcccgggCAGTTCGCCAGTTTGCCCGGCCGGCCCGCGAAGCCGAGCGCCCCTTTCGGGGGGGTATTTTCACGCGCGCGCCGCGCGGTTCCCCGGGAGCGGGCAGCATCGATAATTTTGCAAACCACGTTTTCCGGGTGGGTTGGTTTTTCCAGCAGGTATTCTCCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGggagg                                                                                                                                                                             <  2:628268/234‑1 (MQ=255)
                                                                               gCTAGGTCGAGCGCACCTTTACGGCGGGTCTTTTCCCGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGAc                                                                                                                                           <  1:492065/239‑1 (MQ=255)
                                                                                       cgggccctcttttCCGGGGGGTCATTTCCCGCGCGCGCGGCGCCGCTGCCGGGGCACGGGCAGCCTCGATATTTTTGCCACCCCCTTTTTGCGTTTTGGTTGGTTTTTCCGGCAGGTATTCGCCCACCAGTTCGTTCATCTGTTTTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATc                                                                                                                                    <  2:628023/229‑1 (MQ=255)
                                                                                       cgcgcgcccttttggggggggTTATTTCACGCGCGCGACGCGCGGCTCCAGGGGCAGGGGCACCATCGATATTTTTCCCAACCCCGTTTTTCGCGTCGGTGGGGTTTTCCGGCAGGTATTCTCCCAGCATTTCTTTCATCTGTTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATc                                                                                                                                    <  2:288926/221‑1 (MQ=255)
                                                                                                  tgagggggggTCATTGCCCGCGCGCGCGGCGCCGCTTCACGGGCAGGGGCAGCATCGATAATTTTCCCCCCCACGTTTTCCGCGTCGGTTGGGTTTCCCAGCAGGTTTTCTGCCAGCAGTTCGTTCATCTGTTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCAcg                                                                                                                         <  2:900402/232‑1 (MQ=255)
                                                                                                      ggcggGTCATTCCCCGCGCGCGCCGCGCCGCTTCCCGGGCACGGGCAGCATCGATAATTTTCCCAACCACGATTTCCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCa                                                                                                                     <  1:146091/239‑1 (MQ=255)
                                                                                                       gcggGTCTTTTCACCCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGTTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCa                                                                                                                     <  1:1326172/238‑1 (MQ=255)
                                                                                                       gcggGTCATGTCAGGCGCCGGCCGCGCGGCTACAGGGGCACGGACAGCGTGAATATTTTTCCAACCCACGTTTTTCGCGTCGTTTGGTTTGTCCAGCAGGTTTCCGCCCAGCAGTTCGTTCATCTGTTTTCCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCa                                                                                                                     <  2:1297212/238‑1 (MQ=255)
                                                                                                                     gcgcgggCCGCCCCGCTTCACGGGCCGGGGCAGAGTCGTATATTTTCCCAACCACTATTTTCGCGCTGTTTGGGTTTCCAAGCAGGTATTCTGCAAGCATTTGTTTAATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGgcg                                                                                                      <  2:742935/239‑1 (MQ=255)
                                                                                                                      cgggcgACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGGTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGCAACCAGTTTGTCTTTGGTCTGGGCGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGACTTCACGCGCATCCTAACCGgcggcg                                                                                                      >  2:872418/4‑238 (MQ=255)
                                                                                                                      cgcgcgACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGATAATTTCGGGTCCGGCACTTTTACGGAAACGACCGCAATCAGGCCTTCACGCGCATAGTCACCGGTTgc                                                                                                       >  2:488126/1‑237 (MQ=255)
                                                                                                                        cgcgACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGGGAATTTCGGGTCCGGCACTTTCACGGAAACGCCCGCAATCTGGCCTTCAGGCGCTTCGTCCCCGCTGGCCCtc                                                                                                   >  2:61391/1‑238 (MQ=255)
                                                                                                                          cgACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTg                                                                                                   >  1:322136/1‑237 (MQ=255)
                                                                                                                                         cgggcacgggcaGATTGGATAATTTTCCAACCCACGATTTTCGCGTGGGTTGTGTTTCCCAGCGGGTTTTCGGCCAGCATTTCTTTCTTCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCt                                                                                  <  1:336593/239‑1 (MQ=255)
                                                                                                                                          gggcacgggaaGCATCGATAATTTTGCCACCCACGATTTTCGCGTCGGTGGGTTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCt                                                                                  <  1:1198522/238‑1 (MQ=255)
                                                                                                                                              acgggcaGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAg                                                                              >  2:422835/1‑238 (MQ=255)
                                                                                                                                              acgggcaGCATCGATAATTTTCCCAACCCCGATTTTCGCGTCGTTTGGGTTTTCCAGCAGTTTTTCGGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGc                                                                             <  1:440348/239‑1 (MQ=255)
                                                                                                                                              acggccaGCATCTATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTTTTCTGCCAGCATTTCGTTCATCTGCTGTTAACCCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGc                                                                             <  2:1096902/239‑1 (MQ=255)
                                                                                                                                               cgggcaGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAg                                                                              >  1:1064874/1‑237 (MQ=255)
                                                                                                                                               cgggcaGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAg                                                                              >  2:68176/1‑237 (MQ=255)
                                                                                                                                               cgggcaGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAg                                                                              >  1:783619/1‑237 (MQ=255)
                                                                                                                                                   gagtatgaaTAATGTTGCTATCCACGATTTTCGGTGGGGTTGGGTTTTCCATGGGTGATTCCGTGCACAGTTTGTACATCGGATGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCcttc                                                                         <  1:445884/230‑1 (MQ=255)
                                                                                                                                                    agcaTCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCcttctt                                                                       >  2:1214009/1‑239 (MQ=255)
                                                                                                                                                      caTCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTTGCCTTCTTTg                                                                     >  2:245373/1‑239 (MQ=255)
                                                                                                                                                      caTCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTg                                                                     >  2:893336/1‑239 (MQ=255)
                                                                                                                                                      caTCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTg                                                                     >  2:508849/1‑239 (MQ=255)
                                                                                                                                                       aTCGATAATTTTGCCAACCGCGATTTTCGCGTCGTTTGGTTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCGGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCttt                                                                      <  1:1359505/237‑1 (MQ=255)
                                                                                                                                                       aTCGATAATTTGGCCACCCACGATTTGCGCGTCGGTTGGTTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCTTCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGt                                                                    <  1:532057/239‑1 (MQ=255)
                                                                                                                                                        tagATAATTTTCCCACCCACGATTTTCGCGCCGGTTGTGTTTCCCAGCGGGTATTCTGCCAGCATTTCGTTCATCGGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCcttct                                                                        <  2:946703/232‑1 (MQ=255)
                                                                                                                                                           aTAATTTTGCCAACCACGCTTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTTACCGCAGTTCGTTCATCTGCTGTTCAACCGCCGCTTTCACCTCAGAAGGAACCAGTTTGTCTTTGGTCTGGGAGCGGCATTTCGGGTCCGGCACTTTCACGGCCTCGACCGCAATCCGGCATTCACGCGCATCGTCACCGGGGGCGTTGACTTTCGCTTTTTTGCTGCAGCATTCTTTGTCCAt                                                                >  2:747158/1‑239 (MQ=255)
                                                                                                                                                               ttttGCCAACCACTATTTGCGGGTGGTTTGGGTTTTCCGGAAGGTATTCTGCCAGCAGTTCGTGCATCTGTTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAg                                                            <  1:740951/239‑1 (MQ=255)
                                                                                                                                                                tttgtccAACGACTGTTTTCGCGTCGGTTGTGTTTGCCAGAAGGTATTCTCCCAGCATTTCTTTCAGCTGTTGTTAAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAg                                                            <  1:488126/234‑1 (MQ=255)
                                                                                                                                                                 ttGCCAACCACATTTTTCGGGCCGTTGGGTTTTTCAGCCTGGTTTTCTGCCACGGTTTGTTTTACCTCTGTTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGc                                                          <  2:488704/239‑1 (MQ=255)
                                                                                                                                                                  tGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGt                                                              >  1:946703/1‑234 (MQ=255)
                                                                                                                                                                     ccAACCCCGTTTTCCGCGTGGTTTGGGTTTTCCAGCAGGTATGCGGCGAGCGGTTCGTTCATCTGCTTTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGtt                                                       <  2:410838/239‑1 (MQ=255)
                                                                                                                                                                     ccAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCCGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGCAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTATTTGTCCATGTAGGCGt                                                        >  1:1210082/1‑238 (MQ=255)
                                                                                                                                                                     ccAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATCCTGCCAGCAGTTCGTTCATCTGTTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGtt                                                       <  1:928743/239‑1 (MQ=255)
                                                                                                                                                                     caaaCCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGTTATTCTGCCAGCAGTTCTTTCATCTGCTTTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGtt                                                       <  1:122586/237‑1 (MQ=255)
                                                                                                                                                                      cccccccgtttttTGCGGCGGGTGGGTTTTTCCAGCAGGTTTTCCGCCGGCTGTTCCTTCCTTTGCTTTCCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGtt                                                       <  2:1228690/230‑1 (MQ=255)
                                                                                                                                                                      cAACCACTATTTTCGCGTGGGTTGGGTTTTCCAGCAGGTTTTCGGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGtt                                                       <  1:960852/238‑1 (MQ=255)
                                                                                                                                                                       acccACGTTTTGCCCGGCGGTTGGGTTTTCCAGCAGTTTTCCACCAAGATGTTCGTTAATCGGCTTTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCa                                                     <  1:331492/237‑1 (MQ=255)
                                                                                                                                                                                    gcggcgGTGGGTTTTCCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCGCCTGTTCAGCCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTGTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGGCCATGTAGGCGTTCAGGGTACGGGTCa                                         <  1:1116401/238‑1 (MQ=255)
                                                                                                                                                                                    cgttggtttGTGTTTTCCAGCAGTTATCCCGCCAGGAGTCCGTCCATCTGGTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCAt                                        <  1:376224/232‑1 (MQ=255)
                                                                                                                                                                                            tGGGTTTTCCAGCAGGTTTTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCa                                 <  1:596705/238‑1 (MQ=255)
                                                                                                                                                                                              ggTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACgg                              >  1:1313841/1‑239 (MQ=255)
                                                                                                                                                                                                 ttttCCAGCGGTTATTCTGCCAGCTGTTCGTTTATCTGCTGTACAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGa                             <  1:897731/238‑1 (MQ=255)
                                                                                                                                                                                                 tttcCCAGCAGTTTTTCGGCCGGTAGTTCGTACATTGGCTGATCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGa                             <  1:1064205/238‑1 (MQ=255)
                                                                                                                                                                                                  ttcccAGCGGGTTTTCTCCCCGCAGTTCTTTCTTCTCTTGTTCAACCGCCGATTTCCCCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAg                           <  1:188007/236‑1 (MQ=255)
                                                                                                                                                                                                   ttcccagaggttATCCTGCCGGCAGTGCGTTCTTCTGCGGTTAACCCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAg                           <  1:1357132/228‑1 (MQ=255)
                                                                                                                                                                                                           gggTTTTCTGCCAGCCGTTCTTTCTTCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGcc                     <  1:1244337/236‑1 (MQ=255)
                                                                                                                                                                                                            ggTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAgg                  >  2:842048/1‑239 (MQ=255)
                                                                                                                                                                                                             gTATTCTGCCAGCAGTTCTTTCATCTGCTTTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAgg                  <  1:1057828/238‑1 (MQ=255)
                                                                                                                                                                                                               aTTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGa               >  2:971927/1‑239 (MQ=255)
                                                                                                                                                                                                               aTTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGa               >  2:127132/1‑239 (MQ=255)
                                                                                                                                                                                                                   tGCCAGAATTTCGTTCTTCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAc           <  1:1185375/239‑1 (MQ=255)
                                                                                                                                                                                                                     ccAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGCAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAcc          >  1:1303228/1‑238 (MQ=255)
                                                                                                                                                                                                                       agcagTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAccgcc       >  2:426679/1‑239 (MQ=255)
                                                                                                                                                                                                                       agcagTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGGGAGTAccgc        >  2:791272/1‑238 (MQ=255)
                                                                                                                                                                                                                       agcagTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGCAGCCTTCTTTGTCCATGTAGGCGTTCATGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAcagc        >  2:1357594/1‑238 (MQ=255)
                                                                                                                                                                                                                        gcagTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAccgccg      >  1:206676/1‑239 (MQ=255)
                                                                                                                                                                                                                        gcagTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAccgccg      >  1:77958/1‑239 (MQ=255)
                                                                                                                                                                                                                        gcagTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAccgccg      >  1:784689/1‑239 (MQ=255)
                                                                                                                                                                                                                         cagTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTAccgccg      >  1:536294/1‑238 (MQ=255)
                                                                                                                                                                                                                            ttcgttcATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCGGCCAGGTGAGTACCGCCGTccc  >  1:962044/1‑237 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                        
GAGTCCCCTTCCACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACTTTCACGGAAACGACCGCAATCAGGCCTTCACGCGCATCGTCACCGGTGGCGCTGACTTTGGCTTTTTTGCTGTAGCCTTCTTTGTCCATGTAGGCGTTCAGGGTACGGGTCATCGCCGCACGGAAGCCTGCCAGGTGAGTACCGCCGTCAC  >  NC_000913/3878839‑3879292

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: