Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,621 A→G 10.4% intergenic (+87/‑109) hrpB → / → mrcB putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,6210AG10.4% 244.3 / 17.5 87intergenic (+87/‑109)hrpB/mrcBputative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base A (43/35);  new base G (4/5);  total (47/40)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.62e-03

TGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAG  >  NC_000913/164356‑164895
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                     
tGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTAACCCTCACCCTAACCCTCTCCCTCAAggggaag                                                                                                                                                                                                                                                                               >  2:442918/1‑270 (MQ=255)
   tGGAGTTGCTTTAACCGTCCCAAAGGCCATTACAAATCACACGAGATTGACGCGACGGCTGGAAAGGGGCGTACCGTGGGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCGGGCCGGACGGCCCGGCAAATACTGCACCGACGCGACGGAAGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGTGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgggt                                                                                                                                                                                                                                                                            <  4:446251/274‑2 (MQ=255)
         tGCTTCCACCGTCCCAAAGGCCATTACAGATCACACCAGTGTTGAGCGACTTCTGGAAAGTGGCGTACCGTGAGGTGCAAAAAGGGATGAAAGGGCGTTATCACAAACACGTCTGGCCGGACGACCCGGCAAATACGGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCGTCTCCCTGTGGGGGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTTCACCCTCACCCTTTCCCTCTCCCTCAAGGGAGAGGGGACCGAt                                                                                                                                                                                                                                                                      <  4:259678/274‑1 (MQ=255)
                                      aTCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATAAGCGCGCCCGGTCACACTCACCCTCACCCACTCCCTCAAGTGAGAGGGGACCGATCGGGACAAAATGTTGAGAGATCtcttcctc                                                                                                                                                                                                                                      >  4:112615/1‑275 (MQ=255)
                                                                     gAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCCACAAATTTGAGAGATCTCTTCTCCTGTCTTGTCACCGAAGAAcctgtctcttat                                                                                                                                                                                                       >  4:133498/1‑263 (MQ=255)
                                                                      aGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGCGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCt                                                                                                                                                                                                     >  2:23442/1‑275 (MQ=255)
                                                                         gTACCGTGAGGTGCAAAAAGGGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGt                                                                                                                                                                                                                              <  2:133499/248‑1 (MQ=255)
                                                                               tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTg                                                                                                                                                                                                                                                      >  2:513762/1‑219 (MQ=255)
                                                                                         aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGACCTCAAATATTGAGAGATATCTTCTTCTGTCTGGTAACAGACTAACAAAAAATCCGGCTTTTGCGCCTGAATATTGCg                                                                                                                                                                                  >  4:398795/1‑275 (MQ=255)
                                                                                              gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCGCCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaa                                                                                                                                                                             >  2:209396/1‑275 (MQ=255)
                                                                                                tGANNGGGCGTTATCCCAAACATGGCTGGCCGGACGACCCGGCAAATACTGCANCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgagag                                                                                                                                                                                                                                                 <  4:513757/206‑1 (MQ=255)
                                                                                                    agtgCGTTAGCCCAAGCAGGTCTGGCCGGACGGCCAGGAAAATAAGGGACCGAGGCGAAGGATGAAAAAGTAGTCGTAAGGTTGGTTTTCTCCCTCTCCCGGGGGGGGAGGGCCGGCGTGAGGGCATCAGCGCGCACGTTCACTCTCACCCTAACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAGGAAAAATTTTGAGAGATATCTTCTTCTGTCTTGTACCAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAGAAGc                                                                                                                                                                        <  4:23442/271‑1 (MQ=255)
                                                                                                            aTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCAACCTAACCCCCTCCCTCAAGGGAGAGGGGACCGAACGAGCCAAAAGTTTGAGAGATATCTTCTTCTGTCTTGTTAAAGAAGAACAGACAACCGGGCTTTTGCGCCTGAATATTGAGGAAAAAAAACATGGCCcgg                                                                                                                                                               >  4:85601/1‑275 (MQ=255)
                                                                                                                     aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcg                                                                                                                                                       <  2:112616/274‑1 (MQ=255)
                                                                                                                      tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                                                                                                                                                                                      <  8:315993/149‑1 (MQ=255)
                                                                                                                        tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGCCGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCACCTCCCTGTGGGCGATGGCCTGGGTGAGGGCATCAGCGCGCCAGCTCATCCTCACCCTAACCCTCCCCCTAAAGGGAGAGGGGGCCGACCGAGCACAAACTTGATACGCTCTCTTCTTCTGTCGTTTACCAGACACAACGAAAATAGTGTTTGTGAGCCTGCATATTGCGGGGAAAAAGCAGGACCGGTAATGGccaggac                                                                                                                                                   >  4:193408/1‑270 (MQ=255)
                                                                                                                          tGGCCGGACTACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTAGGAAGGGTGGTTCTCTCCCTCTCCCTCTGGTAGAGAGCCGGGGAGAGGGCAGCAGCGCGAACGTTAACCCTTACCCTAACGCTACCTCTCAAGGGAGAGGGGACCCATGCAGCACAACATTTGAGAGTTATCTTCTTCTGTCCTGTAACAGACGAACAGACAATCCGCGTTTTCCGCTCGAATCGGGCGGTAAAAAAGCACGCCCGGGACGTACCCCGAGcc                                                                                                                                                 >  4:302312/1‑275 (MQ=255)
                                                                                                                              cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                                                                                                                                                                                              <  6:364/149‑1 (MQ=255)
                                                                                                                              cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                                                                                                                                                                                               >  7:150318/1‑148 (MQ=255)
                                                                                                                              cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                                                                                                                                                                                               <  8:150318/148‑1 (MQ=255)
                                                                                                                                    acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGc                                                                                                                                       <  2:85602/275‑1 (MQ=255)
                                                                                                                                                    aCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAtttt                                                                                                                                                                                                                                                       >  7:191583/1‑149 (MQ=255)
                                                                                                                                                        acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCACGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTcc                                                                                                                   >  2:81427/1‑275 (MQ=255)
                                                                                                                                                           cgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGCGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATTTTGCGGAGAAAAAGCATTGCCGGGAATGACCGCGAGACAAGTGGACGCAAAGGGAAACCGACCCGTCCggg                                                                                                                >  4:215484/1‑274 (MQ=255)
                                                                                                                                                           cgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGat                                                                                                                                                                                                                                               <  6:419271/149‑1 (MQ=255)
                                                                                                                                                           cgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGGGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGGGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGt                                                                                                                <  2:302316/275‑1 (MQ=255)
                                                                                                                                                            gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                                                                                                                                                                                              >  7:307252/1‑119 (MQ=255)
                                                                                                                                                            gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                                                                                                                                                                                              <  8:307252/119‑1 (MQ=255)
                                                                                                                                                            gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTc                                                                                                               <  2:416390/275‑1 (MQ=255)
                                                                                                                                                                          tATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCtt                                                                                                                                                                                                                                >  7:467874/1‑149 (MQ=255)
                                                                                                                                                                           aTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTgg                                                                                                                                           >  2:474755/1‑232 (MQ=255)
                                                                                                                                                                               gTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                                                                                                                                                                                               >  8:375870/1‑99 (MQ=255)
                                                                                                                                                                               gTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                                                                                                                                                                                               <  7:375870/99‑1 (MQ=255)
                                                                                                                                                                               gTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                                                                                                                                                                                           >  5:140671/1‑149 (MQ=255)
                                                                                                                                                                                tAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCcgtc                                                                                            <  2:193410/274‑1 (MQ=255)
                                                                                                                                                                                  aGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACaga                                                                                                                                                                                                                        <  6:353147/149‑1 (MQ=255)
                                                                                                                                                                                     ttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                                                                                                                                                     <  5:47830/118‑1 (MQ=255)
                                                                                                                                                                                     ttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTga                                                                                                                                                                                                                                                     >  6:47830/1‑118 (MQ=255)
                                                                                                                                                                                     ttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaaga                                                                                                                                                                                                                    >  5:270941/1‑149 (MQ=255)
                                                                                                                                                                                       ggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAAc                                                                                                                                                                                                                  <  8:191583/149‑1 (MQ=255)
                                                                                                                                                                                            tctccttcccctCGTGGGAGAGGGGCGGCGGGAGGTCGTCCGCGCGCGCGTTCACCCGAACCCTAACCCTATCCCTCAAGGGGGAGGGGACCGATCGAGCACAAATTTTGAGAGGTATCTTCTTCTGTCTTGTAAAAGAAGAACAGAAAATCGGGCTTGTGCGCCTGAATAGTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGTAAAGGGAAACCGACGCGTCCGGTCAATCAAAAGGGAAGCCGTCGTCGTTACGAAGa                                                                               <  4:156042/275‑1 (MQ=255)
                                                                                                                                                                                                   ctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAAc                                                                                                                                                                                                                           >  2:523507/1‑129 (MQ=255)
                                                                                                                                                                                                       ctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACagaa                                                                                                                                                                                                                      <  4:523502/129‑1 (MQ=255)
                                                                                                                                                                                                       ctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCtttt                                                                                                                                                                                                  >  5:202333/1‑149 (MQ=255)
                                                                                                                                                                                                        tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                                                                                                                                                                                          >  6:390959/1‑94 (MQ=255)
                                                                                                                                                                                                        tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAAt                                                                                                                                                                                                                                                          <  5:390959/94‑1 (MQ=255)
                                                                                                                                                                                                         gtgGGGGAGGGACGGGGCTAGGGGTTCAGCGCGCACGTCCAACGTCCACCTAACCCATGCCCTAAATGGAGAGGGGAACGATCGAGAACAAATTTTGAGAGTTATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGTTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTac                                                                  <  4:81426/275‑1 (MQ=255)
                                                                                                                                                                                                         gtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTgc                                                                                                                                                                                                <  6:270941/149‑1 (MQ=255)
                                                                                                                                                                                                          tgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTCCGAAGATGACgataatga                                                                   >  2:219723/1‑271 (MQ=255)
                                                                                                                                                                                                           gggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTacga                                                                <  2:249337/275‑1 (MQ=255)
                                                                                                                                                                                                              agagGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGCAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCCAGCCAATTGGAGGCAAAGGGAAACCGACGCGTCCGGTCAACCAAAAGGTCAGCCGGCGTCGTTACGAAGATGACGATGAttcagaaga                                                             >  2:61443/1‑268 (MQ=255)
                                                                                                                                                                                                               gagGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                                                                                                                                                                                                    <  6:294698/93‑1 (MQ=255)
                                                                                                                                                                                                               gagGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTgag                                                                                                                                                                                                                                                    >  5:294698/1‑93 (MQ=255)
                                                                                                                                                                                                                agGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCAGGTCAAACAAAAGGTAAGCAGTCGTCGTTACGAAGACGACGATACTTacgacagtt                                                           >  2:85314/1‑271 (MQ=255)
                                                                                                                                                                                                                   gCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTatg                                                        <  2:215487/275‑1 (MQ=255)
                                                                                                                                                                                                                              ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCt                                                                                                                                                                                                     >  4:105758/1‑123 (MQ=255)
                                                                                                                                                                                                                               gCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa                                                                                                                                                                          <  6:140671/149‑1 (MQ=255)
                                                                                                                                                                                                                                  tCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                                                                                                                                                                                        >  8:466149/1‑116 (MQ=255)
                                                                                                                                                                                                                                  tCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCgg                                                                                                                                                                                                        <  7:466149/116‑1 (MQ=255)
                                                                                                                                                                                                                                  tCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTg                                                                                                                                                                                                 <  2:105759/123‑1 (MQ=255)
                                                                                                                                                                                                                                   cAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCg                                                                                                                                                                                                                                                                        >  4:406163/1‑54 (MQ=255)
                                                                                                                                                                                                                                    aGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATg                                                                                                                                                                     <  8:467874/149‑1 (MQ=255)
                                                                                                                                                                                                                                       gcgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCg                                                                                                                                                                                                                                                                    <  2:406168/54‑1 (MQ=255)
                                                                                                                                                                                                                                           aCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGCAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGa                                                                                                                                                              >  7:260745/1‑149 (MQ=255)
                                                                                                                                                                                                                                           aCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGa                                                                                                                                                              <  5:234477/149‑1 (MQ=255)
                                                                                                                                                                                                                                              ttCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg                                                                                                                                                                         >  5:74316/1‑135 (MQ=255)
                                                                                                                                                                                                                                              ttCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg                                                                                                                                                                         <  6:74316/135‑1 (MQ=255)
                                                                                                                                                                                                                                              ttCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCGTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATg                              <  2:161412/274‑1 (MQ=255)
                                                                                                                                                                                                                                                   cctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATAACcgc                                                                                                                                                      >  5:16406/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         cctAACCCTCTCCCTCAAGGGAGAGGGTACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGt                  <  2:20756/275‑1 (MQ=255)
                                                                                                                                                                                                                                                          ctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                                                                                                                                               >  5:264178/1‑149 (MQ=255)
                                                                                                                                                                                                                                                             aCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGAGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGAtt                                                                    >  1:53795/1‑221 (MQ=255)
                                                                                                                                                                                                                                                              cccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                                                                                                                                                    <  7:236097/140‑1 (MQ=255)
                                                                                                                                                                                                                                                              cccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                                                                                                                                                    >  8:236097/1‑140 (MQ=255)
                                                                                                                                                                                                                                                               ccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                                                                                                                                                                                                 <  6:182159/78‑1 (MQ=255)
                                                                                                                                                                                                                                                               ccTCTCCCTTGAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACa                                                                                                                                                                                                                 >  5:182159/1‑78 (MQ=39)
                                                                                                                                                                                                                                                               ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGCGCCAATTGGa                                                                                                                                          <  8:260745/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa                                                                                                                                          >  7:356398/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               ccTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATACTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa                                                                                                                                          >  8:458277/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGa                                                                                                                                                              >  7:469806/1‑127 (MQ=255)
                                                                                                                                                                                                                                                                 tctcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGa                                                                                                                                                              <  8:469806/127‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGCGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTACg                                                                                                                  >  1:12496/1‑169 (MQ=255)
                                                                                                                                                                                                                                                                   tcCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTAAGACGATTATGATGAATATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAg               >  1:63575/1‑268 (MQ=255)
                                                                                                                                                                                                                                                                   ccccaaaaGGGAGGGGGGACCGAGAGNNCAAAAATGTTGGCAGATTGTTTCTTCTGTCNTGTAACAGAAGAAACGAAAATCGGGCTTTGGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCNAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGggca           <  3:63575/266‑1 (MQ=255)
                                                                                                                                                                                                                                                                     ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCaaa                                                                                                                                    >  6:192461/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     ccTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCaaa                                                                                                                                    >  7:23391/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       tCAAGGGAGAGGGGACCGATCGAGCACAACTTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCTGGCTTTTGCGCCTGAATATGGCGGAGAAAAAGCATGGCCGGGACTGACGGCGAGACAATTAGAAGCAAAGGGAGACAGACGCATACGGTAACACAAAAGGTACGCGGACTGCGTTACGAAAATTACGAAGATTCCGACGTGTATGCTACCTATATCGATGAAAAACGGATTCCGAGCAGAGGAAAgagataagcgag    >  4:498795/1‑264 (MQ=255)
                                                                                                                                                                                                                                                                       tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg                                                                                                                                  >  8:419722/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       tCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGggcaaaggcaa    <  2:398800/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aaGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCAAATTGGACGCAAAAGGAAACCGACGCGACAGGTCCAACAAAAGGTAAGCCGTCGTCGTTACGAAGNAGACGATGATTACGACGAATATGATGACGAAGAGGCTGAAGAATAGATGCCACGAAAAGGGAAGggcaaagagaaaa  >  3:111148/1‑269 (MQ=255)
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                     
TGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCACAAATTTTGAGAGATATCTTCTTCTGTCTTGTAACAGAAGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAAAG  >  NC_000913/164356‑164895

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