Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,622,348:1 +A 5.1% coding (115/639 nt) purN → phosphoribosylglycinamide formyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,622,3481.A 5.1% 1068.9 / 53.4 255coding (115/639 nt)purNphosphoribosylglycinamide formyltransferase 1
Reads supporting (aligned to +/- strand):  ref base . (124/118);  new base A (6/7);  total (130/125)
Fisher's exact test for biased strand distribution p-value = 7.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-03

Only 100 of 259 total aligned reads displayed.

CACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGAAAATCGGTATCATCAAAGCCTCTGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTTTATGCGCATTCTCAGCCCGGCGTTTGTCTCCCACTATGCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCG  >  NC_000913/2622084‑2622618
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                              
cACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTACTCAATGCCAACGGTGAAAACGCGTGGAAAATCGGTATCATCAAAGCCTCTGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGAGTTCg                                                                                                                                                                                                                                                                       >  2:113255/1‑274 (MQ=255)
                                                      cAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGAAAATCGGTATCATCAAAGCCTCTGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTCTTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACCCATACGCTCATCCCCAgcac                                                                                                                                                                                                                >  4:357449/1‑273 (MQ=255)
                                                             cTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGAAAATCGGTATCATCAAAGCCTCTGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGa                                                                                                                                                                                                          <  2:94348/274‑1 (MQ=255)
                                                                                aCGGTGAAAACGCGTGGAAAATCGGTATCATCAAAGCCTCTGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGAcc                                                                                                                                                                                      >  2:27934/1‑275 (MQ=255)
                                                                                     gAAAACGCGTGGAAAATCGGTATCATCAAAGCCTCTGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGa                                                                                                                                                                                  <  2:484663/274‑1 (MQ=255)
                                                                                                             atcaAAGCCTCTGATTCCGAAAAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTa                                                                                                                                                          <  2:332900/274‑1 (MQ=255)
                                                                                                             atcaAAGCCTCTGATCCCTAACAATGCGTGGTTATCGAATATAGAATATTGTGGGGCTTATGTCCGGAAACGGAAGGAAGTAACAGGTAATTATTGACGCCTGTAAAAACAACAAAATTAAAGGCACCGTACGGGCAGGTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTGGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGAGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTa                                                                                                                                                          <  4:234663/274‑1 (MQ=255)
                                                                                                                     ctctGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑a                                                                                                                                                                                                                                                                               >  8:277693/1‑149 (MQ=255)
                                                                                                                      tctGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑Ac                                                                                                                                                                                                                                                                              >  6:800469/1‑149 (MQ=255)
                                                                                                                       ctGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGAtgt                                                                                                                                              >  1:159492/1‑276 (MQ=255)
                                                                                                                        tGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTTCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCCGTCTTAAGCAATCAGCCAG‑Aca                                                                                                                                                                                                                                                                             >  4:552880/1‑147 (MQ=255)
                                                                                                                         gATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑Acgcg                                                                                                                                                                                                                                                                           >  8:1002829/1‑149 (MQ=255)
                                                                                                                           ttCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGtt                                                                                                                                                                                                                                                                         >  5:1202395/1‑149 (MQ=255)
                                                                                                                            tCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTc                                                                                                                                                                                                                                                                        <  6:961083/149‑1 (MQ=255)
                                                                                                                            tCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTc                                                                                                                                                                                                                                                                        <  8:1482493/149‑1 (MQ=255)
                                                                                                                            tCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCa                                                                                                                                                      <  2:541559/263‑1 (MQ=255)
                                                                                                                              cGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCgg                                                                                                                                                                                                                                                                      >  7:1251353/1‑149 (MQ=255)
                                                                                                                               gAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGc                                                                                                                                                                                                                                                                     >  8:1107569/1‑149 (MQ=255)
                                                                                                                                  caacGTGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCG‑TCGGCCtt                                                                                                                                                                                                                                                                  <  7:371796/149‑1 (MQ=255)
                                                                                                                                  caacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCtt                                                                                                                                                                                                                                                                  <  7:1107569/149‑1 (MQ=255)
                                                                                                                                   aacGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCG‑TCGGCCtta                                                                                                                                                                                                                                                                 <  7:275647/149‑2 (MQ=255)
                                                                                                                                      gcgTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTg                                                                                                                                                                                                                                                                 >  7:653542/1‑146 (MQ=255)
                                                                                                                                          ggTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGACGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACgcgc                                                                                                                                                                                                                                                          <  7:1453024/149‑1 (MQ=255)
                                                                                                                                              aTCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCAt                                                                                                                                                                        >  2:37010/1‑227 (MQ=255)
                                                                                                                                                   ataatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGc                                                                                                                                                                                                                                                 >  7:222184/1‑149 (MQ=255)
                                                                                                                                                     aatGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCgg                                                                                                                                                                                                                                               <  6:1046181/149‑1 (MQ=255)
                                                                                                                                                      atGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCGGGCTGGTTTTATGCGCATtct                                                                                                                >  2:625039/1‑275 (MQ=255)
                                                                                                                                                       tGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGt                                                                                                                                                                                                                                             >  6:1532821/1‑149 (MQ=255)
                                                                                                                                                       tGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGt                                                                                                                                                                                                                                             >  5:833344/1‑149 (MQ=255)
                                                                                                                                                        gAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACTCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTa                                                                                                                                                                                                                                            <  7:877654/149‑1 (MQ=255)
                                                                                                                                                           tatTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCG‑TCGGCCttat                                                                                                                                                                                                                                                                <  5:309204/126‑3 (MQ=255)
                                                                                                                                                           tatTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCgg                                                                                                                                                                                                                                               <  6:1641138/143‑1 (MQ=255)
                                                                                                                                                           tatTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTg                                                                                                                                                                                                                                         >  7:435309/1‑149 (MQ=255)
                                                                                                                                                              tgtgGTGCNTATTTCCGNNAACGGAAGTAATTTACAGGCAATNATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGc                                                                                                                                                                                                                                        <  3:64881/147‑1 (MQ=255)
                                                                                                                                                               gtggtgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAAc                                                                                                                                                                                                                                     <  6:1202395/149‑1 (MQ=255)
                                                                                                                                                                  gtgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAAcg                                                                                                                                                                                                                                                             <  4:655936/122‑1 (MQ=255)
                                                                                                                                                                  gtgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCa                                                                                                                                                                                                                                  >  6:1504643/1‑149 (MQ=255)
                                                                                                                                                                   tgCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCAt                                                                                                                                                                                                                                 >  8:704770/1‑149 (MQ=255)
                                                                                                                                                                      ttATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCt                                                                                                                                                                                                                                                                   <  8:454175/112‑1 (MQ=255)
                                                                                                                                                                       tATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTg                                                                                                                                                                                                  <  4:186582/176‑1 (MQ=255)
                                                                                                                                                                          ttCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCCACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATTTACGCACCCGATGTGTTCGGTCTGGCTGGTTTTATGCGCATTTTCAGCCCGGCGTTTGtccc                                                                                              >  4:433231/1‑271 (MQ=255)
                                                                                                                                                                            ccGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATc                                                                                                                                                                                                                        <  8:577716/149‑1 (MQ=255)
                                                                                                                                                                              ggCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGc                                                                                                                                                                                                                      >  5:202092/1‑149 (MQ=255)
                                                                                                                                                                              ggCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTTTATGCGCATTCTCAGCCCGGCGTTTGtc                                                                                                <  2:433235/267‑1 (MQ=255)
                                                                                                                                                                               gCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGcc                                                                                                                                                                                                                     >  6:369709/1‑149 (MQ=255)
                                                                                                                                                                                cAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATc                                                                                                                                                                                                                        >  2:689566/1‑145 (MQ=255)
                                                                                                                                                                                cAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCa                                                                                                                                                                                                                    <  8:1251353/149‑1 (MQ=255)
                                                                                                                                                                                 aaCGGAAGTAATTTACAGGCAATTATTGACGCCTGT‑AAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAgc                                                                                                                                                                                                                  >  6:996770/1‑149 (MQ=255)
                                                                                                                                                                                   cGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATa                                                                                                                                                                                                                                <  2:212459/134‑1 (MQ=255)
                                                                                                                                                                                      aaGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGt                                                                                                                                                                                                              <  5:1528444/149‑1 (MQ=255)
                                                                                                                                                                                      aaGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTTTATGCGCATTCTCAGCCCGGCGTTTGTCTCCCACTATGCCGGGc                                                                                <  2:657193/275‑1 (MQ=255)
                                                                                                                                                                                         tAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTg                                                                                                                                                                                                           >  8:950132/1‑149 (MQ=255)
                                                                                                                                                                                          aaTTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGa                                                                                                                                                                                                          >  6:297827/1‑149 (MQ=255)
                                                                                                                                                                                           aTTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGAc                                                                                                                                                                                                         >  5:784633/1‑149 (MQ=255)
                                                                                                                                                                                            tttACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACa                                                                                                                                                                                                        <  5:1622087/149‑1 (MQ=255)
                                                                                                                                                                                            tttACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACa                                                                                                                                                                                                        <  6:1363177/149‑1 (MQ=255)
                                                                                                                                                                                                   gCAATTATTGACGCCTGTAAACCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCGAACGCG‑TCGGCCttatt                                                                                                                                                                                                                                                               >  6:58568/1‑84 (MQ=38)
                                                                                                                                                                                                   gCAATTATTGAAGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGa                                                                                                                                                                                                          >  7:958947/1‑140 (MQ=255)
                                                                                                                                                                                                    cAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGaa                                                                                                                                                                                                >  6:1613086/1‑149 (MQ=255)
                                                                                                                                                                                                             aCGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGtt                                                                                                                                                                                                                                                                         >  1:44082/1‑67 (MQ=255)
                                                                                                                                                                                                                ccTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCgg                                                                                                                                                                                    <  8:605901/149‑1 (MQ=255)
                                                                                                                                                                                                                ccTGAAAAACCAACAAATTTAAAGGCACATTACGGGCAGTTTACAGTAATAAGGCAG‑ACGCGTTCGGCCTTGCACGCTGCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACACTCGTGAAGCCTATGCCCgg                                                                                                                                                                                    <  6:532650/149‑1 (MQ=255)
                                                                                                                                                                                                                 cTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCtt                                                                                                                                                                                                                                                                  >  5:494494/1‑70 (MQ=255)
                                                                                                                                                                                                                  tGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGa                                                                                                                                                                                  >  8:421587/1‑149 (MQ=255)
                                                                                                                                                                                                                  tGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACg                                                                                                                                                        <  1:185860/175‑1 (MQ=255)
                                                                                                                                                                                                                       aaCCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATAGCCAGCGCGTTTGACAGTCGTGAAGCTTATGACCGGGAGTTGa                                                                                                                                                                             <  8:39637/149‑1 (MQ=255)
                                                                                                                                                                                                                           aacaaAATTAAGGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCa                                                                                                                                                                         >  7:1352464/1‑149 (MQ=255)
                                                                                                                                                                                                                            acaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATAc                                                                                                                                                                                                                               >  6:1301407/1‑94 (MQ=255)
                                                                                                                                                                                                                            acaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGa                                                                                                                                                                                  <  7:856916/139‑1 (MQ=255)
                                                                                                                                                                                                                            acaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGt                                                                                                                                                                                >  6:869250/1‑141 (MQ=255)
                                                                                                                                                                                                                            acaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGt                                                                                                                                                                                <  5:869250/141‑1 (MQ=255)
                                                                                                                                                                                                                            acaaAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCAt                                                                                                                                                                        >  8:1325388/1‑149 (MQ=255)
                                                                                                                                                                                                                              aaaaTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGa                                                                                                                                                                      <  7:500762/149‑1 (MQ=255)
                                                                                                                                                                                                                              aaaaTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGa                                                                                                                                                                      >  5:621488/1‑149 (MQ=255)
                                                                                                                                                                                                                               aaaTTAAAGGCACCGTACGGGCAGTTTTCATCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGt                                                                                                                                                                                                   <  8:1160664/119‑1 (MQ=255)
                                                                                                                                                                                                                               aaaTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGt                                                                                                                                                                                                   >  7:1160664/1‑119 (MQ=255)
                                                                                                                                                                                                                                aaTTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATa                                                                                                                                                                                                                                <  4:749726/89‑1 (MQ=255)
                                                                                                                                                                                                                                  ttAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATc                                                                                                                                                                  >  7:719542/1‑149 (MQ=255)
                                                                                                                                                                                                                                  ttAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATc                                                                                                                                                                  <  8:1640225/149‑1 (MQ=255)
                                                                                                                                                                                                                                   caagggtgccGTACGGGTAGTTTTTAGCAATAGNCCCG‑GCGCGTTCGGCCTTGACCGCGNCCGCCCAGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCNTGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTTTATGCGCATTCTCAGCCCGGCGTTTGTCTCCCACTATGCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTaca                               <  3:163085/271‑1 (MQ=255)
                                                                                                                                                                                                                                       ggCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAg                                                                                                                                                                                 >  5:433599/1‑129 (MQ=255)
                                                                                                                                                                                                                                       ggCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAg                                                                                                                                                                                 <  6:433599/129‑1 (MQ=255)
                                                                                                                                                                                                                                       ggCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGtt                                                                                                                                                                               >  8:1432830/1‑131 (MQ=255)
                                                                                                                                                                                                                                       ggCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACAt                                                                                                                                                             >  6:664237/1‑149 (MQ=255)
                                                                                                                                                                                                                                              tACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCa                                                                                                                                                      >  8:1585324/1‑149 (MQ=255)
                                                                                                                                                                                                                                               aCGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCAc                                                                                                                                                     <  5:1349763/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    cAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGAt                                                                                                                                                <  5:59160/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    cAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGAt                                                                                                                                                <  7:1629710/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       ttttCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGAtgtg                                                                                                                                             <  6:855492/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            aGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGAtt                                                                                                                                                                           <  8:870714/114‑1 (MQ=255)
                                                                                                                                                                                                                                                              cAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgc                                                                                                                                      <  8:1083533/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               aaTAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgct                                                                                                                                     >  5:1011076/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                aTAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctg                                                                                                                                    >  6:762798/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 tAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctgg                                                                                                                                   >  7:1631496/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 tAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctgg                                                                                                                                   <  5:996770/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aGGCCGAACGCG‑TCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctggct                                                                                                                                 >  8:31090/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   aGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTgctggct                                                                                                                                 >  7:1070983/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    ggCCC‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGCATTGCAACGCATACGCTCCTCGCCCGCGCCATTGACAGCCGTGAAACCTATGACCGGGAGTTGATTCATGACATCGACATGCACGCACCCGATGTGGTCGTGCGCGCTTGCTTTATGCGCATGCTCAGCACGGCGTCTGTCTCCCACCATACCGAGCGTTTGATGCACAGTCCCCCTTCTCGGCTCCCGACATCACCCGGTTTCCATACCCATCCTCATGCGCTGGAGATTGGca  >  4:517341/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                       cG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGtt                                                                                                                             >  5:901946/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        g‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGttt                                                                                                                            >  5:1051998/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                              
CACCTTCAACTGCGGCGTCGGGATGATTATTGCCCTGCCTGCTCCGGAAGTGGACAAAGCCCTCGCCCTGCTCAATGCCAACGGTGAAAACGCGTGGAAAATCGGTATCATCAAAGCCTCTGATTCCGAACAACGCGTGGTTATCGAATAATGAATATTGTGGTGCTTATTTCCGGCAACGGAAGTAATTTACAGGCAATTATTGACGCCTGTAAAACCAACAAAATTAAAGGCACCGTACGGGCAGTTTTCAGCAATAAGGCCG‑ACGCGTTCGGCCTTGAACGCGCCCGCCAGGCGGGTATTGCAACGCATACGCTCATCGCCAGCGCGTTTGACAGTCGTGAAGCCTATGACCGGGAGTTGATTCATGAAATCGACATGTACGCACCCGATGTGGTCGTGCTGGCTGGTTTTATGCGCATTCTCAGCCCGGCGTTTGTCTCCCACTATGCCGGGCGTTTGCTGAACATTCACCCTTCTCTGCTGCCGAAATATCCCGGATTACACACCCATCGTCAGGCGCTGGAAAATGGCG  >  NC_000913/2622084‑2622618

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑