Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,608,229 A→G 14.7% M439V (ATG→GTG)  hyfF → hydrogenase 4, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,608,2290AG14.7% 144.4 / 2.0 79M439V (ATG→GTG) hyfFhydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (29/32);  new base G (9/4);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCAT  >  NC_000913/2607999‑2608466
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   
tGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTgcg                                                                                                                                                                                                                                             <  2:104193/239‑1 (MQ=255)
tGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTgcg                                                                                                                                                                                                                                             <  2:1205680/239‑1 (MQ=255)
tGCTAAAGCACGGCGCGCGCGATCTCCACGTGGTCTGTGGGATGCTCAAAATCATGCCCTTGGCCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTgcg                                                                                                                                                                                                                                             <  1:656549/239‑1 (MQ=255)
  cTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTgcggg                                                                                                                                                                                                                                           >  2:714331/1‑237 (MQ=255)
   tCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCgg                                                                                                                                                                                                                                          <  2:328249/239‑1 (MQ=255)
    cAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCggg                                                                                                                                                                                                                                         <  2:928399/239‑1 (MQ=255)
            gCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCCCGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGGCTGGTGCGGGTGGGCGCGGGGGGGTTTATg                                                                                                                                                                                                                                 >  1:467738/1‑239 (MQ=255)
            gCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGGCCGGTACGGATGGCTGGGCGGGGGTTAATg                                                                                                                                                                                                                                 >  2:250552/1‑239 (MQ=255)
                                gtcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGGTGGCTGCGCGGGTGTTTATGGCGAAACCGGCGCAGGGCGt                                                                                                                                                                                                          >  2:539782/1‑239 (MQ=255)
                                gtcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGCTGGCTGGGCGGGGGTTAATGGCGCAACCGCCGCCGGCCGt                                                                                                                                                                                                          >  2:971344/1‑239 (MQ=255)
                                gtcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCCGCTTCTCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCGCCGCACAGACGGCCGCGCGGGTGGTagtcggggagacgccgccgagcgg                                                                                                                                                                                                             >  1:208262/1‑216 (MQ=255)
                                 tcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGGTGGCTGCGCGGGGGTTAATGGCGAAACCGCCGCAGGCCGtt                                                                                                                                                                                                         >  2:398736/1‑239 (MQ=255)
                                          tGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGGCCGGTACGGGTGGCTGCGGCGGGGTTGATGGCGAAACCGCCGCCGGGCGTCACACGGGGt                                                                                                                                                                                                >  2:257889/1‑239 (MQ=255)
                                           gCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGGACGGGTGGCGGCGCGGGGGTTTATGGCGGAACCGCCCCCGGGCGGTAACCgggggg                                                                                                                                                                                               >  1:61257/1‑237 (MQ=255)
                                            cTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGCCTGGCACGTAATCACCTGCTGCTTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGGCCGGGAGGGCTCGGCGGGGTTTATGGGGCAACCGCCCCCGGTCCGTAACCCgggggg                                                                                                                                                                                               >  1:1204440/1‑237 (MQ=255)
                                            cTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGGACCGGTGGCTGCGCGGGGGTTAATGGCGGAAACGCCGCGGGCCGTTTAACCgggggg                                                                                                                                                                                               >  1:813643/1‑237 (MQ=255)
                                                    gctGCCATTTCGCGCCCGTGGTGTTGGGGGGCGGGGGCTGGCGCCTGACGGGGATCCCCCCCTGAAAATTTTTCCTAGGATATTTGATGCCCTTTACCGCCGGCCTGGCAGGTAACCCCCGGCGGATTACCGTCCTGCTGTTTTTGCGGTTATCGCGGGTGCTGGGGGGCCTGTTAGGGAGGGCTGCGCGGGTGTTAAGGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGc                                                                                                                                                                                      <  1:66894/237‑1 (MQ=255)
                                                      tGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGGCTGGGACGGATGGCTGGGCGGGGGGTTATGGCGGAACCCCCGCAGGCCGTTTACCGGGGGGAATCTCgctgg                                                                                                                                                                                    >  1:1208061/1‑235 (MQ=255)
                                                      tGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCCGGTAAGGGTGGCTGCGCGGGGGTTTATGGCGGAACCCGCGCCGGGCGGTTAACCgggcggtgctcgccgg                                                                                                                                                                                     >  1:395561/1‑226 (MQ=255)
                                                       gCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGATCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGGGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGGGATCTCGGCTGGt                                                                                                                                                                                  >  2:436020/1‑239 (MQ=255)
                                                           tttACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCCGGCCGTTAACCGGGGTGATCTCGGCTGGTTgg                                                                                                                                                                               >  2:626860/1‑237 (MQ=255)
                                                           ttcaCCGCCGTGCTGTGTGGGGGCGGTGCGCTGGGGCGGGCAGGTAGGCCCCCCTTCACCATTTTTCTTAGGGAATTTATGCCCTTTCCCGCCGGACTGGCACGTAATCACCTGCTGATTAGCGTCCTGCTGTTTTTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGAGGGCTGCTGGGGTGTTAAGGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGac                                                                                                                                                                              <  1:626860/236‑1 (MQ=255)
                                                                cgccgTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGGATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCGTGGTACGGTTGGCTGCGCGCGTGTTTATGGTGGAACCACCGCAGGCCGTTTACCGGGGTGGTGTCGGGTGGTTGtccacc                                                                                                                                                                          >  1:647571/1‑238 (MQ=255)
                                                                      tGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGGTGGCTGCGCGGGTGTTAATGGCGGAACCGCCGCAGGCCGTTAACCGGGGGGATCTCGGCTGGGTGACCACCCCGCCa                                                                                                                                                                   >  2:70663/1‑239 (MQ=255)
                                                                                                        ggATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTAAGGATGGTTGCGCGGGGGTTTATGGCGGAACCGCCGCAGGCCGGTAACCGGGGGGGTCTCGGCCGGGTGCCCACCACCGCCAATGGGAATCTACTGGTCCTTATGCTGGTGATg                                                                                                                                  >  1:912921/1‑239 (MQ=255)
                                                                                                         gATGCCGCCCTTCAACATTTTTCTAAGCGAATTTATGCCCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGTTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGa                                                                                                                                    <  2:688214/236‑1 (MQ=255)
                                                                                                              cgccCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTCATGGCGGAAACCCCCCCGGCCGGTAACCGGGGGGATCTCGGCCGGTTTGCCACCCCGCCCATGGGTAGTCTCGTGGTCATGTTTCTGGCGACGGGAAc                                                                                                                             >  1:525615/1‑238 (MQ=255)
                                                                                                              cgccCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAAc                                                                                                                             >  2:645274/1‑238 (MQ=255)
                                                                                                                   ttCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGGCGTTAAACGGGGGGGTTTAGGGTGGTTGAACACCCCGCCAATGGGGATTCTGCTGGTCATGGTGCTGGGGGTTGGAACGCaaat                                                                                                                       >  1:172456/1‑239 (MQ=255)
                                                                                                                   ttCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCCGGTACGGATGGCTGCGCGGGGGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGGTCTCGGCTGGTTGCCCACCTCCCCCATGCTGGTTCTGCGGGGCATTGTGCTGGCGGAGGGGACGccaat                                                                                                                       >  1:837186/1‑235 (MQ=255)
                                                                                                                             tttttAGCGATTTTATGCCAATTCCCGCCGGACTGGCCCTTAATCACCTGCTAATTTTCGTCTTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGCCCTGGTACGGATGGCGGCGCGGGTGTTAATGGCGAAGCCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCt                                                                                                             <  2:45105/236‑1 (MQ=255)
                                                                                                                                  aGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGGTGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGGGATTCTGCTGGTCATGATGCTGGCGATGGGGACGCATATTCCACAACCTGtca                                                                                                         >  2:34285/1‑238 (MQ=255)
                                                                                                                                  aGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGGGTTAATGGCGAAACCGCCGCAGGCCGCTAATCGGGGGGATTCCGGGTGGGGGGACAACTCGCACAAGGGGTCTCAGCTGGTGCTGTTGCTGGCGGTGGGAACGCATATTTcacagaacgtaa                                                                                                         >  1:326483/1‑230 (MQ=255)
                                                                                                                                                   ttACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGGGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGGGATCTCGGCTGGGTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGGTGGGAACGCCTATTCCACAACCTGTCATCAGGGACCGtgcgggc                                                                                        >  1:1032927/1‑239 (MQ=255)
                                                                                                                                                    tACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGGGTTAATGGCGGAAACGCCGCAGGCCGGTAACCGGGGGGGTCTCGGGCGGGTGGCCCCCACCCCCATGGGGATGTTGGTGGTCATGGTTGTGGCGGTTGGGAAGCGTATTCCACCACCTTGCAACAGGGTCCTggcgggg                                                                                        >  1:353189/1‑237 (MQ=255)
                                                                                                                                                     cccgccgGACTGGCACGTAATCACCTGCGGATTATCGTCCTGCTGTTATTGCTGTTAACGCGGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGgcgc                                                                                     <  2:98772/238‑1 (MQ=255)
                                                                                                                                                     ccccccgGACTGGCACGTAATCCCCTGCTTATTATCGTCCTGCTTTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTCGGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGgcgc                                                                                     <  2:215196/235‑1 (MQ=255)
                                                                                                                                                       ctccgggtgggCTCGTAATCACCGTCTTATTAGCGTTGTGGTTTGTTTGTGTTTAACGTTGGTGCTGCTGGGCCTGGTAGGTAGGGGCGCGCGGGTGTTAATGGCAAAACCGCCGCAGGCCGTTAATCGGGGTGATCTCGGCTGGTTGCCCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCtt                                                                                   <  1:953264/230‑1 (MQ=255)
                                                                                                                                                       acctggAGTGTCACGAAAGCACCTGCTGATTATCGTCCTGTTGTTATTGCTGTTAAGGCTGCAGCGGGCGGGCCTGGTACGAGGGGCGGGGCGGGTGTTAATGGCGAACCCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCtt                                                                                   <  1:539782/235‑1 (MQ=255)
                                                                                                                                                        gcgggACGGGCACGTACTCCCGTGTTGATTTTCGCCCGGTTTTTTTTGCTGTGAACGCTGGTGCTGGCGGGCCTGGTCCGGTTGGCTGCGCGGGTTTAAATGGCGAACCCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTc                                                                                  <  1:436020/236‑1 (MQ=255)
                                                                                                                                                        gccgGACTGGCTCGTAATCACCTGCTAGTTATCGTCGTGCTTTTATTGCTGTTAGCGCTGGGGCTGGCGGGCCTGGTCGGGATGGCTGCGCGGTTGTTAATGGCGAAGCCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTc                                                                                  <  2:8815/239‑1 (MQ=255)
                                                                                                                                                          cggactgcacctcatcctccTGCTGTTTTTCGTCTTGTTGTTATGGTTGTTAACGCGGGTGCTGGCGGGCCTGGTAGGGATGGCGGCGCGGGTTTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCa                                                                                <  1:351440/221‑1 (MQ=255)
                                                                                                                                                          cgacgtggacagTAATCACTGGTTGATTTTCTGTCTGGTTTTATTGCTTTTAGCTCGGGTTCTGGCGCCCCTGGTCCGGGTGGCTGCGCGGGTGTTAATGGGGAACCCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCa                                                                                <  1:70663/228‑1 (MQ=255)
                                                                                                                                                           ggACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTTATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGGTTCTGCGGGGCCTGAGGCTGGCGATTGGAACCCACATTCCACAACCCGTCATCAGGAAACTGGGGGGCGGCTTcc                                                                                 >  1:879963/1‑238 (MQ=255)
                                                                                                                                                           ggACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGGGATCTCGGCTGGTTGGCCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTcc                                                                                 >  2:1037625/1‑237 (MQ=255)
                                                                                                                                                           gcatggcaacgtacccccgtgctgATTAGGGTCCTGTTGTTGTTGCTGTTACGGTGGGTGCGGGCGGGCGGGGTGCGGATGGCGGGGTGGGTTTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCa                                                                                <  1:95018/219‑1 (MQ=255)
                                                                                                                                                                ggCAGGTAACCCCCTGCTGTTTGTCCTCGTGCTGTTTTTGCTGTTACCGCTGTTGCGGGGGGGCCTGGTACGGATGGCGGCGCGGGTGTTAAGGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACtata                                                                           <  1:250552/238‑1 (MQ=255)
                                                                                                                                                                 gCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGCTGGCTGCGCGGGGGTTAAAGGCGAAACCGCCGCCGGGCGTTAAACGGGGGGGTCTCGGGTGGTTGACCACCTCGCCCATGGGGAATCTGGTGGTCTTGATTCGGGCGCAGGGAAAGCAACTTACCAAACCCGGCATCGCGAACCTGGCGGGTGCGTCCCCtctca                                                                          >  1:837307/1‑234 (MQ=255)
                                                                                                                                                                  cACGTAATCACCTGCTGATTATCGTCCGGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGt                                                                         <  2:1070859/238‑1 (MQ=255)
                                                                                                                                                                  cACGTAAGCACTTGTTGATTACCGTCTTGCTGTTATTGCTGTTAAGGCTGGGGCTGGCGGGCCTGGTACGAATGGCTGCGCGGGTGTTAATGGGGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGac                                                                                                                                                                              <  2:402042/139‑1 (MQ=255)
                                                                                                                                                                          cACCTTGTGGTTATCTTCCTGTTGTAATTGGTGTAAGGCCTGGTGCTGGCGGGCCGGGTACGGTTGCGTGCGCGGGTGTAAAGGGCAAACCCGCCGCAGGCCGTTAACCGGGGTGATCGCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAgg                                                                <  1:714331/239‑1 (MQ=255)
                                                                                                                                                                                ctgCTTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGGGTTAATGGCGACACCGCCGCAGGCCGTTTACCGGGGTGATCTCGGCTGGGTTACCACCTCCCCAATGGCGATTCTGCTGGGCATACGGCTGGCGATTGGAACGCCTATTTCCCCACCTGCCCTCAGGTGCCCGCCCGGGGCTTCCCCTCTACTCCTCTCCGGGaacc                                                           >  1:791854/1‑235 (MQ=255)
                                                                                                                                                                                ctgATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGGGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGGTGACCACCCCGCCCATGGTGATTCTGCTGGTCAATCTGCTGGGGATGGGAAACCCAATTTCACAACCTGTCATCCGGGTACTGGCGGGGGCGTCCCCATTAGTCGTCTCACGGACGc                                                           >  1:721590/1‑238 (MQ=255)
                                                                                                                                                                                   aTTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGa                                                       <  2:61257/239‑1 (MQ=255)
                                                                                                                                                                                   aTTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGGACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGGGATCTCCGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGGCATGATGCTGGCGGTGGGGACGCATATTCCACCAACCGTCATCAGGGTCCTGGCGGGCGCGTCCACTATATTCCtctcccggccgcccg                                                        >  1:443413/1‑227 (MQ=255)
                                                                                                                                                                                      aTCGTCCTGCTGTTATTGCTGTTACCGCTGGTGGTGGCGGGGCTGGTACGGATGGCTGCGCGGTTGTTAAGGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATct                                                    <  2:879963/239‑1 (MQ=255)
                                                                                                                                                                                        cGTCCTGCTGTTTTGGCTGTTAACGCTGGTGCTGGCGGGCCGGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATctgc                                                  <  2:813643/239‑1 (MQ=255)
                                                                                                                                                                                             ttgtgtgttttgtGTTAACGCTGGTGCTGGGGGGCCTGGTACGGGTGGCGGCGGGGGTTTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGca                                             <  2:572578/227‑1 (MQ=255)
                                                                                                                                                                                             tgttgTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTTTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGca                                             <  2:1006854/236‑1 (MQ=255)
                                                                                                                                                                                             tgctgTGATTGCGGTTAACGCGGGTGCTGGCGGGCCTGGTACGGATGGCGGCGCGGGTGTAAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGca                                             <  2:534792/239‑1 (MQ=255)
                                                                                                                                                                                             gtcggTTTTTGCTGTGACCTCTGGGGCTGGCGCGCTTGGTACGGGTGTCGGCCCGGGTGTATAGGTCAAAACCGCCGCAGGCCGTTAACCGGGTTGATCACGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGTTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCATGATCCTTGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGca                                             <  1:523471/235‑1 (MQ=255)
                                                                                                                                                                                                ttttttggCTGTAAAGGCGGGGGCTGGGGGGCCTGTTAGGGAGGGCGGGGCGGGTGTGAAGGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGcaca                                           <  1:971344/231‑1 (MQ=255)
                                                                                                                                                                                                tgTTTTTGTTGTTACCGTTGGTGCTGGGGGGCCTGGTAGGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGcaca                                           <  2:173819/238‑1 (MQ=255)
                                                                                                                                                                                                tgTTATTGCTGTTAACGCTGGTGCGGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGtcatca                                                                                                      <  2:1147717/180‑1 (MQ=255)
                                                                                                                                                                                                  ttATTGCTGTTAACGCTGGTGCTGGCGGGGCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAAACGGGGTGGTCTCGGCTGGTTGACCAACCCGCCAATGGTGAATCTGCTGGGCATGATGGTGGCGATGGGAACGCATATTCCACAACCTGTCTCCAGGGTGCTGGTGGGCGCGTCCACTAAAATCCGCTCCGGGACGccccgtctccccgccccacc                                        >  1:911525/1‑220 (MQ=255)
                                                                                                                                                                                                   tATTGCTGTTACCGCTGGTGCTGGGGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACg                                        <  2:1139306/238‑1 (MQ=255)
                                                                                                                                                                                                    aTTGCTGTTAACGCTGGTGCTGGCGGGCCCTGTACGGATGGCTCCGCGGGTGTTAATGGCGGAACCGCCGCAGGCCGGTAACCCGGGGGGTCTCGGCTGGGTGACCACCCCGCCAATGGGGATTCTGCCGGGCCTGATGCCGGCGGTGGGCACCCATATTTCACAACCTGTCATCAGGATCCCTGCCGGCGCCTCCCCACTTACACTCTCAGGGACGCAccacctttctgtcccccct                                       >  1:688459/1‑220 (MQ=255)
                                                                                                                                                                                                     ttGCTGTTACCGCTGGTGCTGGCGGGCTTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAg                                     <  2:443413/239‑1 (MQ=255)
                                                                                                                                                                                                         tGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTCACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCAc                                  >  2:550758/1‑238 (MQ=255)
                                                                                                                                                                                                          gTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGGGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGCGCTGCCTGCACAACGTAGcc                                   >  1:1171744/1‑235 (MQ=255)
                                                                                                                                                                                                            tAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTgcg                                                                                                                                                                                                                                             <  1:887675/38‑1 (MQ=255)
                                                                                                                                                                                                            tAACGCTGGTGCTGGCGGGCCCGGTACGGATGGCTgcg                                                                                                                                                                                                                                             >  2:887675/1‑38 (MQ=255)
                                                                                                                                                                                                                  tGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGCCGCACGATCTGCCTGCACAACGGAGCACCTGGCATGAt                        >  2:384090/1‑239 (MQ=255)
                                                                                                                                                                                                                     tGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGGGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCCCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAAACCGTCATCAGGATCCTGGCGGGCGCCTCCCCCATAGTCCTCTCCAGGACGCACCACCTGCCTGCCCAACCTCACACCCGGCATGTtttt                     >  1:980744/1‑239 (MQ=255)
                                                                                                                                                                                                                          gCGGGCCTGGTAGGGATGCCGGCGGGGGTGTAAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACt                                                                              <  1:528700/177‑1 (MQ=255)
                                                                                                                                                                                                                                cgggTAGGGAGGGCTGCGGGGGTGTAAAGGGCGAACCCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAgg          <  1:645274/237‑1 (MQ=255)
                                                                                                                                                                                                                                    tACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCctct                                                                    >  2:528700/1‑177 (MQ=255)
                                                                                                                                                                                                                                        ggATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGt                                                                         <  2:1195699/169‑1 (MQ=255)
                                                                                                                                                                                                                                        ggATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGt                                                                         >  1:1195699/1‑169 (MQ=255)
                                                                                                                                                                                                                                          aTGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCAt  <  2:172456/238‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   
TGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCAT  >  NC_000913/2607999‑2608466

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