Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,762,328 A→G 6.0% F62S (TTC→TCC)  sufC ← SufBCD Fe‑S cluster assembly scaffold protein, ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,762,3280AG 6.0% 225.7 / 3.9 67F62S (TTC→TCC) sufCSufBCD Fe‑S cluster assembly scaffold protein, ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (29/34);  new base G (2/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

CGTTGCGCTTTTTCTCGCCGCCGGAAAAACCAACGTTTACCGAACGGGTTAATAAATCTTCCGGCATCTTCAGGAGAGCGATTTTCTCTTCCATCAAATCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTAATACTTAACATGTTTATTCCTTATCCGACGCTGTGTTCAAGACTGATGGCGAGGAGTTTTTGTGCTTCAACGGCAAATTCCAACGGCAGCTCCGAGAACACGTCTTTGCAGAA  >  NC_000913/1762050‑1762614
                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                              
cGTTGCGCTTTTTCTCGCCGCCGGAAAAACCAACGTTTACCGAACGGGTTAATAAATCTTCCGGCATCTTCAGGAGAGCGATTTTCTCTTCCATCAAATCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGt                                                                                                                                                                                                                                                                                      >  1:435129/1‑289 (MQ=255)
      gcTTTTTCTCGCCGCCGGAAAAACCAACGTTTACCGAACGGTTTAATAAATCTTCCGGCATCGTCAGGAGAGCGATTTTCTCTTCCATCAAATCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTgccgcc                                                                                                                                                                                                                                                                                <  1:117911/289‑1 (MQ=255)
       cTTTTTCTCGCCGCCGGAAAAACCAACGTTTACCGAACGGGTTAATAAATCTTCCGGCATCTTCAGGAGAGCGATTTTCTCTTCCATCAAATCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGAAACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGccc                                                                                                                                                                                                                                                                               >  1:563291/1‑289 (MQ=255)
         ttttCTCGCCGCCGGAAAAACCAACGTTTACCGAACGGGTTAATAAATCTTCCGGCATCTTCAGGAGAGCGATTTTCTCTTCCATCAAATCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGt                                                                                                                                                                                                                                                                             <  2:402508/289‑1 (MQ=255)
          tttCTCGCCGCCGGAAAAACCAACGTTTACCGAACGGGTTAATAAATCTTCCGGCATCTTCAGGAGAGCGATTTTCTCTTCCATCAAATCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTc                                                                                                                                                                                                                                                                            <  2:518946/289‑1 (MQ=255)
          tttCTCGCCGCCGGAAAAACCAACGTTTACCGAACGGGTTAATAAATCTTCCGGCATCTTCAGGAGAGCGATTTTCTCTTCCATCAAATCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTc                                                                                                                                                                                                                                                                            <  2:563291/289‑1 (MQ=255)
           tgctcGCCCCCGGAAAAACCAACGTTTACCGAACGGGTTAATAAATCTTCCGGCATCTTCAGGAGAGCGATTTTCTCTTCCATCAAATCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCa                                                                                                                                                                                                                                                                           <  1:588727/287‑1 (MQ=255)
                                 cGTTTACCGAACGGGTTAATAAATCTTCCGGCATCTTCAGGAGAGCGATTTTCTCTTCCATCAAATCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGc                                                                                                                                                                                                                                                     <  1:499073/289‑1 (MQ=255)
                                             gggTTAATAAATCTTCCGGCATCTTCAGGAGAGCGATTTTCTCTTCCATCAAATCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGa                                                                                                                                                                                                                                         <  1:361168/289‑1 (MQ=255)
                                                  aataaATCTTCCGGCATCTTCAGGAGAGCGATTTTCTCTTCCATCAAATCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAgg                                                                                                                                                                                                                                    <  1:138749/289‑1 (MQ=255)
                                                    taaATCTTCCGGCATCTTCAGGAGAGCGATTTTCTCTTCCATCAAATCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGt                                                                                                                                                                                                                                   >  1:469494/1‑288 (MQ=255)
                                                    taaATCTTCCGGCATCTTCAGGAGAGCGATTTTCTCTTCCATCAAATCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACCCAACCGTGCCGCCCGTCACATCATAATCGTCTCGCCAGGAAAGCGTTGCCGATAAggga                                                                                                                                                                                                                                  >  2:657521/1‑287 (MQ=255)
                                                         cTTCCGGCATCTTCAGGAGAGCGATTTTCTCTTCCATCAAATCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTCCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGACGAGGTACtttt                                                                                                                                                                                                                             >  1:87131/1‑289 (MQ=255)
                                                            ccGGCATCTTCAGGAGAGCGATTTTCTCTTCCATCAAATCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGc                                                                                                                                                                                                                           >  2:158169/1‑288 (MQ=255)
                                                                 agctgcAGGAGAGCGAATTTCTATGCCATTAAATCCTGAAAATCAAAGCGGTGCAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTGTGCAGGAAAAACTGGTTACTGACACCTTGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCGTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCTCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGaa                                                                                                                                                                                                                      <  2:317968/283‑1 (MQ=255)
                                                                         gagaGCGATTTTCTCTTCCATCAAATCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTgg                                                                                                                                                                                                              >  1:357099/1‑288 (MQ=255)
                                                                          agagCGATTTTCTCTTCCATCAAATCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTgg                                                                                                                                                                                                              <  2:435129/287‑1 (MQ=255)
                                                                            agCGATTTTCTCTTCCATCAAATCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACt                                                                                                                                                                                                                                >  1:30085/1‑267 (MQ=255)
                                                                            agCGATTTTCTCTTCCATCAAATCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACt                                                                                                                                                                                                                                <  2:30085/267‑1 (MQ=255)
                                                                                   ttCTCTTCCATCAAATCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATg                                                                                                                                                                                                    >  2:212381/1‑288 (MQ=255)
                                                                                                  tCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAg                                                                                                                                                                                                                                                  >  2:102605/1‑227 (MQ=255)
                                                                                                  tCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAg                                                                                                                                                                                                                                                  <  1:102605/227‑1 (MQ=255)
                                                                                                       aaaaTCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGGACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGa                                                                                                                                                                               >  2:220337/1‑289 (MQ=255)
                                                                                                            cAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTc                                                                                                                                                                          <  1:158169/289‑1 (MQ=255)
                                                                                                                          gTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCcg                                                                                                                                                             <  2:357099/288‑1 (MQ=255)
                                                                                                                            ttCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTTGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCcgcg                                                                                                                                                           <  1:493881/288‑1 (MQ=255)
                                                                                                                            ttCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCcgcg                                                                                                                                                           <  2:469494/288‑1 (MQ=255)
                                                                                                                            ttCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCcgcg                                                                                                                                                           <  1:212381/288‑1 (MQ=255)
                                                                                                                             tCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATCCTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACCTCTTAATCGTCTCGCCCGGCTAGATTTGCCGATAAGTTACTTTTGCCCGAACCGTTTCGCCCCATAACGGCGTGAAATTCGCCGGGATGAACGTCAAGGCTTAATCCgcgcc                                                                                                                                                         >  2:182533/1‑288 (MQ=255)
                                                                                                                                                    cGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTg                                                                                                                                  <  1:174647/289‑1 (MQ=255)
                                                                                                                                                          aaGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACgtgt                                                                                                                            <  1:503441/289‑1 (MQ=255)
                                                                                                                                                              gCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGGACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTaa                                                                                                                          <  1:220337/287‑1 (MQ=255)
                                                                                                                                                              gCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTaaa                                                                                                                         >  1:494452/1‑288 (MQ=255)
                                                                                                                                                                            aaaCTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAACTCTTTAATACTTaat                                                                                                          >  1:482114/1‑288 (MQ=255)
                                                                                                                                                                                  tttccttacacCTTGCATCTCCACCGGATACTGGAAGGCCATTAAGATGCCTTCGCCCGCCCGATCTTCCGTCGACATCTCAAGCAAAGCTTTTCCTTTTAGCTCAACCGTGCCCCCCGTCACTTCATCATCTTCTCTCCCGGCAAGCGTTGCCGATAAGCTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGCGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTTACGTGTAAATCTTTAATACTTAACATGttt                                                                                                    <  2:16039/282‑1 (MQ=255)
                                                                                                                                                                                            ccTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTAATACTTAACATGTTTATTCCTTAt                                                                                           >  1:638284/1‑288 (MQ=255)
                                                                                                                                                                                                  aTCGCCACGGGAAACTGGAGGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTGTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTAATACTTAACATGTTTATTCCTTATCCGACg                                                                                     <  2:482114/288‑1 (MQ=255)
                                                                                                                                                                                                   tctcCACCGGATACTGGAAGGGCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTAATACTTAACATGTTTATTCCTTATCCGACGCt                                                                                   <  1:313667/289‑1 (MQ=255)
                                                                                                                                                                                                   tctcCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTAATACTTAACATGTTTATTCCTTATCCGACGCt                                                                                   <  1:626141/289‑1 (MQ=255)
                                                                                                                                                                                                     tcCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTa                                                                                                                   >  1:43592/1‑255 (MQ=255)
                                                                                                                                                                                                     tcCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTa                                                                                                                   <  2:43592/255‑1 (MQ=255)
                                                                                                                                                                                                                      ggCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTAATACTTAACATGTTTATTCCTTATCCGACGCTGTGTTCAAGACTGATGGCg                                                                >  2:526601/1‑289 (MQ=255)
                                                                                                                                                                                                                      ggCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTAATACTTAACATGTTTATTCCTTATCCGACGCTGTGTTCAAGACTGATGGCg                                                                >  2:72707/1‑289 (MQ=255)
                                                                                                                                                                                                                         cATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTAATACTTAACATGTTTATTCCTTATCCGACGCTGTGTTCAAGACTGATGGCgag                                                              <  1:526601/288‑1 (MQ=255)
                                                                                                                                                                                                                           tAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTAATACTTAACATGTTTATTCCTTATCCGACGCTGTGTTCAAGACTGATGGCgaggag                                                           >  1:409344/1‑289 (MQ=255)
                                                                                                                                                                                                                                  gCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGGACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAAt                                                                                                                                                                >  2:190157/1‑181 (MQ=255)
                                                                                                                                                                                                                                  gCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGGACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAAt                                                                                                                                                                <  1:190157/181‑1 (MQ=255)
                                                                                                                                                                                                                                  gCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCcgcg                                                                                                                                                           <  2:53160/186‑1 (MQ=255)
                                                                                                                                                                                                                                  gCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCcgcg                                                                                                                                                           >  1:53160/1‑186 (MQ=255)
                                                                                                                                                                                                                                  gCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTaa                                                                                                                          >  2:82210/1‑219 (MQ=255)
                                                                                                                                                                                                                                  gCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTaa                                                                                                                          <  1:82210/219‑1 (MQ=255)
                                                                                                                                                                                                                                  gCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTAATACTTAACATGTTTATTCCTTATCCGACGCTGTGTTCAAGACTGATGGCGAGGAGTTTTtgt                                                    >  2:279562/1‑289 (MQ=255)
                                                                                                                                                                                                                                       cGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTg                                                                                                                                  <  1:285916/206‑1 (MQ=255)
                                                                                                                                                                                                                                       cGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTg                                                                                                                                  >  2:285916/1‑206 (MQ=255)
                                                                                                                                                                                                                                       cGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTAATACTTAACATGTTTATTCCTTATCCGACGCTGTGTTCAAGACTGATGGCGAGGAGTTTTTGTGCtt                                                >  1:592953/1‑288 (MQ=255)
                                                                                                                                                                                                                                                    ttCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTAATACTTAACATGTTTATTCCTTATCCGACGCTGTGTTCAAGACTGATGGCGAGGAGTTTTTGTGCTTCAACGGCAAATTcc                                  <  2:409344/289‑1 (MQ=255)
                                                                                                                                                                                                                                                        ggCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCttt                                                                                                                    >  1:558216/1‑203 (MQ=255)
                                                                                                                                                                                                                                                        ggCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCttt                                                                                                                    <  2:558216/203‑1 (MQ=255)
                                                                                                                                                                                                                                                           gACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTAATACTTAACATGTTTATTCCTTATCCGACGCTGTGTTCAAGACTGATGGCGAGGAGTTTTTGTGCTTCAACGGCAAATTCCAACGGc                            <  2:638284/288‑1 (MQ=255)
                                                                                                                                                                                                                                                            aCAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTAATACTTAACATGTTTATTCCTTATCCGACGCTGTGTTCAAGACTGATGGCGAGGAGTTTTTGTGCTTCAACGGCAAATTCCAACGGCa                           >  2:654725/1‑288 (MQ=255)
                                                                                                                                                                                                                                                              aGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTAATACTTAACATGTTTATTCCTTATCCGACGCTGTGTTCAAGACTGATGGCGAGGAGTTTTTGTGCTTCAACGGCAAATTCCAACGGCAGc                         >  2:434046/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                cgcAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTAATACTTAACATGTTTATTCCTTATCCGACGCTGTGTTCAAGACTGATGGCGAGGAGTTTTTGTGCTTCAACGGCAAATTCCAACGGCAGCTc                       >  2:432233/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                 gcaagcaaATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTAATACTTAACATGTTTATTCCTTATCCGACGCTGTGTTCAAGACTGATGGCGAGGAGTTTTTGTGCTTCAACGGCAAATTCCAACGGCAGCTcc                      <  1:72707/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aagcaaATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTAATACTTAACATGTTTATTCCTTATCCGACGCTGTGTTCAAGACTGATGGCGAGGAGTTTTTGTGCTTCAACGGCAAATTCCAACGGCAGCTCCgag                   <  2:494452/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                          tCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTAATACTTAACATGTTTATTCCTTATCCGACGCTGTGTTCAAGACTGATGGCGAGGAGTTTTTGTGCTTCAACGGCAAATTCCAACGGCAGCTCCGAGAACACg             >  1:535760/1‑288 (MQ=255)
                                                                                                                                                                                                                                                                             ttGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGTTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTAATACTTAACATGTTTCTTCCTTATCCGACGCTGTGTTCAAGCCTGATGGCGAGGAGTTTTTGTGCTTCAACGGCAACTTCCAACGGCAGCCCCGAGGACACGTCtt         >  2:42012/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                                    tGAACTCAACCGTGCCTCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTAATACTTAACATGTTTATTCCTTATCCGACGCTGTGTTCAAGACTGATGGCGAGGAGTTTTTGTGCTTCAACGGCAAATTCCAACGGCAGCTCCGAGAACACGTCTTTGCAGaa  <  1:612652/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                              
CGTTGCGCTTTTTCTCGCCGCCGGAAAAACCAACGTTTACCGAACGGGTTAATAAATCTTCCGGCATCTTCAGGAGAGCGATTTTCTCTTCCATCAAATCCTGAAAATCAAAGCGGTCGAGCGTTTCCTGGCCGCGATAGCTGCGCACCGCATTAAGTGCCGTTTGCAGGAAAAACTGGTTACTGACACCTGGAATCTCCACCGGATACTGGAAGGCCATAAAGATGCCTTCGCCCGCGCGATCTTCCGGCGACAGCGCAAGCAAATCTTTGCCTTTGAACTCAACCGTGCCGCCCGTCACTTCATAATCTTCTCGCCCGGCAAGCGTTGCCGATAAGGTACTTTTGCCCGAACCGTTTGGCCCCATAATGGCGTGAACTTCGCCGGGATGAACGTCGAGGCTTAATCCGCGCAGGATAGCTTTATCTTCCACGCTGACGTGTAAATCTTTAATACTTAACATGTTTATTCCTTATCCGACGCTGTGTTCAAGACTGATGGCGAGGAGTTTTTGTGCTTCAACGGCAAATTCCAACGGCAGCTCCGAGAACACGTCTTTGCAGAA  >  NC_000913/1762050‑1762614

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