Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,080,362 A→C 15.4% V77G (GTG→GGG)  yeeA ← putative transporter, FUSC family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,080,3620AC15.4% 150.9 / 4.3 80V77G (GTG→GGG) yeeAputative transporter, FUSC family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (41/23);  new base C (2/13);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 4.46e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCAGTAAAGAGCCGAGGATCACATCGCCGCTTCGCCATAACGCCGTATCAATTTCACCTGTCGGGGAACCCACAACAATTGCCAGCGTCACCCCAATCAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGATAGTAAACAGGCGGATAATGAGAAAAGTGAGCAGGAATGCCAGCGCGACCCGAGTACCATGCACAATGCGGTAATGGCGGTATACCCGGATTTCAAACGGGCTTAATGACTTATCGGCACGCACGCAGAAA  >  NC_000913/2080127‑2080598
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              
tGCCAGTAAAGAGCCGAGGATCACATCGCCGCTTCGCCATAACGCCGTATCAATTTCACCTGTCGGGGAACCCACAACAATTGCCAGCGTCACCCCAATCAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCCCCGCGCACCAGACTAACATCACCGGTAACGAGATTAACTCCAGTTGCAGAGCGAAAAAACCTAAAAACAAAACAAACAcca                                                                                                                                                                                                                                             >  2:808921/1‑238 (MQ=255)
           agCCGAGGATCACATCGCCGCTTCGCCATAACGCCGTATCAATTTCACCTGTCGGGGAACCCACAACAATTGCCAGCGTCACCCCAATCAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCCGACTAACATCAGCGGTAACGCGATTAACTCCAGTTGCAGAGCGCTAAGACCCAAAACCCAACCCAAAAccctccacttacaa                                                                                                                                                                                                                                  >  2:78838/1‑227 (MQ=255)
            gCCGAGGATCACATCGCCGCTTCGCCATAACGCCGTATCAATTTCACCTGTCGGGGAACCCACAACAATTGCCAGCGTCACCCCAATCAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCACCGGTAACGAGATTAACTCCAGTTGCCAAGCGATAAGACCTAAAATCCAACCCAACAcccaaccccccaaa                                                                                                                                                                                                                                  >  2:87328/1‑226 (MQ=255)
            gCCGAGGAACACAGCGCCGTTTCCCAATAACGCCGTATCAATTTCACCTGCCGGGGAACCCACAACAATTGCCAGCGTCACCCCAATCAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACCCCGTACCGCCAAt                                                                                                                                                                                                                                   <  1:637279/237‑1 (MQ=255)
                aGGATCACATCGCCGCTTCGCCATAACGCCGTATCAATTTCACCTGTCGGGGAACCCACAACAATTGCCAGCGTCACCCCAATCAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAATTTAACTCCAGTTGCAGACCGAAAAGACCCAAAACCGAACCAAACACCGAACGCCCAAAAcccac                                                                                                                                                                                                                             >  2:908750/1‑235 (MQ=255)
                aGGATCACATCGCCGCTTCGCCATAACGCCGTATCAATTTCACCTGTCGGGGAACCCACAACAATTGCCAGCGTCACCCCAATCAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGGAACGAGATTAACTCCAGTTGCAGAGCGCTAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTc                                                                                                                                                                                                                             >  1:300875/1‑239 (MQ=255)
                       caTCGCCGCTTCGCCATAACGCCGTATCAATTTCACCTGTCGGGGAACCCACAACAATTGCCAGCGTCACCCCAATCAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCACCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAGGACCTAAACTCGAACCCAACACCGTACCGCCAATACGCTCCAGGgcg                                                                                                                                                                                                                      >  1:145541/1‑239 (MQ=255)
                         tCGCCGTTTCCCCATAACCCCGTTTCAATTTCACCTGGCGGGGAACCCACAACGATTGCCAGCGTCCCCCCAATCAATAACCCTTGATCCGGTTTCTTGCCCGGCGCCGGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTACCATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACGGTACCGCCAATACGCTCAAAGGcgcg                                                                                                                                                                                                                    <  1:109034/239‑1 (MQ=255)
                          cgccgcTTCGCCATAACGCCGTATCAATTTCACCTGTCGGGGAACCCACAACAATTGCCAGCGTCACCCCAATCAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTACAGAGATTAACTCCAGTTCCAAAGCGAAAAACCCTAAAATCGAACCAAAAACCGCACCGCCAATACGCCCAAAGGagcgg                                                                                                                                                                                                                   >  2:923033/1‑238 (MQ=255)
                          cgccgcTTCGCCATAACGCCGTATCAATTTCACCTGTCGGGGAACCCACAACAATTGCCAGCGTCACCCCAATCAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAAACCTAAAATCGAACCCAACAACGAACCGCCAAAAAGCCCAAAGgcgcaa                                                                                                                                                                                                                   >  2:329783/1‑237 (MQ=255)
                              ccgctgccacaccgcctttacAATTCCGCCTGCGGGGGAACCCAAAACAATTCCCAGCGCCCCCCCAAACAAAAAACCTTTAACCGGTTTCTTCCCCAGCGCCACCCACCCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACCACAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGAAAAGCCCTAAAATCGAACCCAACCCCGTACCGCCAATACGCTCAAAGGCGCGAGGGa                                                                                                                                                                                                               <  2:675598/219‑1 (MQ=255)
                                  cGCCATAACGCCGTATCAATTTCACCTGTCGGGGAACCCACAACAATTGCCAGCGTCACCCCAATCAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAAAACCGCACCACCCAAACGCTCAAAAGCGCGCGGGACaaa                                                                                                                                                                                                           >  2:727047/1‑238 (MQ=255)
                                   gccacacccccctaccAATTTCCACTGTCGGGGAACCCCAAAAAATTCCCAGCGTCACCCCAATCAATAACCCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCCCCAGACTAACACCAGCGGTAACGAGATTAACCCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACg                                                                                                                                                                                                          <  2:1011231/224‑1 (MQ=255)
                                    tataagagccagaTCTATTCCCCCTCTCGGGGAACCCACACCAGTTCCCCGGGCCCCCCCAATCAATAAACCTTGATCCGGCTTCTTGCCCAGCCCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGAAACGAGATTAACCCCGGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGAACCGCCAATACGCTGAAAGGCGCGAGGGACAACGt                                                                                                                                                                                                         <  2:24290/227‑1 (MQ=255)
                                              gTATCAATTTCACCTGTCGGGAAACCCACAACAATTGCCAGCGTCACCCCAATCAATAACCCTTGATCCGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGaa                                                                                                                                                                                               <  1:912726/239‑1 (MQ=255)
                                               tATCAATTTCACCTGTCGGGGAACCCACAACAATTGCCAGCGTCACCCCAATCAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGa                                                                                                                                                                                                >  1:1011231/1‑237 (MQ=255)
                                               tATCAATTTCACCTGTCGGGGAACCCACAACAATTGCCAGCGTCACCCCAATCAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAAc                                                                                                                                                                                              >  1:274060/1‑239 (MQ=255)
                                                    ttttCCCCTTTGGGGTAACCCACAACAATGCCAAGCGTCACCCCAAAAAATAAACCTTGATCCGGTTTCTTGCCCAGCGCCAGCCAATCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCATAGCGATAAGACCTAAAATCGAACCCAACACCGTGCCGTCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGAt                                                                                                                                                                                         <  2:302641/238‑1 (MQ=255)
                                                     tttccccgcgcGGGGAACCCACAACAATTGCCAGCGTCCCCCCAATCAATAAACCTGGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGAtt                                                                                                                                                                                        <  1:234087/229‑1 (MQ=255)
                                                         cccTGTCGGGGAACCCACAACAATTGCCAGCGTCACCCCAATCAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGcc                                                                                                                                                                                    <  1:329783/238‑1 (MQ=255)
                                                         aCCTGTCGGGGAACCCACAACAATTGCCAGCGTCACCCCAATCAATAACCCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCCAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAAAACCGTACCGCCAAAACCCCCAAAGCCCCGAGGGCACACGTTACCCCAAAACAAAATTGGc                                                                                                                                                                                     >  2:594894/1‑238 (MQ=255)
                                                           cTGTCGGGGAACCCACAACAATTGCCAGCGTCACCCCAATCAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGaa                                                                                                                                                                                               >  1:24290/1‑226 (MQ=255)
                                                               gggggAACCCAAAACAATTGCCAGCGCCCCCCCAACCAATAACCCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCCCCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACCCCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATa                                                                                                                                                                              <  1:993677/237‑1 (MQ=255)
                                                                ggggCACCCACAACAATTGCCGAGTCCTCCCGAAACAAAAAAACGGAGTACGGTGTTGGTCCCAGCGCCAGCCAACCGCAAAGGAACATGCCTGCGGCGGACCAGACAAACATCAGCGAGACCAAGATGAACTCCGGTTGCAGAGCGATAAGCCCTAAAATCGAACCCAACCCCGTACCGCCAATACGCTCAAAGGCGCGAGGGAGAACGTTACCCCAGAACGAGATTGGCCCCATagt                                                                                                                                                                            <  2:154913/239‑3 (MQ=255)
                                                                gggaAACCCCCAACAATTCCCACGTCCACCCCAAACAAAAAACCTTGATCCGGTTTCTTCCACAGCCCCAGCCAACCGAAAAGAAACAGGGCCCCCGCGCCCCAGACTACCATCAGCGGAAACGAGATTAACCCCAGTTGAAGACCGATAAGACCTAAAATCGACCCCAACCCCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAAt                                                                                                                                                                            <  2:693673/239‑1 (MQ=255)
                                                                 ggaaaCCCACAAAACGTGCCAGCGTCCCCCCAGTCTAGAACCCTGGGTACATTGTCTTCCCCAGCGCAGCCCAACTGCAAAGGAACATGGCCGCCGCACACCAGACTAACATCACCGGTAACTAGATGAACTCTGTTTGCAGAACAGTAAGTCCTAAAGGCGAACCCACCTCCGTACCGCCAATACGCCCAAAGGCGCGAGCGACAACGTTACCCCAGAACGAGATTGGCCCCATAAt                                                                                                                                                                            <  2:1066040/235‑1 (MQ=255)
                                                                   gcaCCCTCAACAATTCCCATCGTACACCCAAACAAAAACCCTTTAACCGGTTCCTTGCCCGGCACCAGCCAGCCGAAAAGGAACATGGCCGCCGCGCACCGGACAAACCTCAGGGGTAACGAGATTAACCCCAGGTGCTGAGCGATAAGACCAAAAACCGAACCAACCCCCGGCCCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCATAATcac                                                                                                                                                                         <  2:300875/236‑1 (MQ=255)
                                                                     aCCCACAACAATTGCCAGCGTCACCCAAAACAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATcacca                                                                                                                                                                       <  1:1011563/239‑1 (MQ=255)
                                                                        cacaACAATTGCCAGCGTCACCCCAATCAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCaccacc                                                                                                                                                                     >  1:2995/1‑238 (MQ=255)
                                                                        cacaACAATTGCCAGCGTCACCCCAATCAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCAccacca                                                                                                                                                                    >  1:605335/1‑239 (MQ=255)
                                                                        cacaACAATTGCCAGCGTCACCCCAATCAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAAATTGGCCCCATAATCAcaacca                                                                                                                                                                    >  1:447246/1‑239 (MQ=255)
                                                                                ttGCCAGCGTCACCCCAATCAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTACCATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTCAAATCGAACCCAACACCGTACCGCCAATACGCTCTAAGGCGCGAGGGACAACGTTCACCCAGACCGAGATTGGCCCCATAATCACCCCCATGGTGAcc                                                                                                                                                            >  1:943501/1‑239 (MQ=255)
                                                                                 tGCCAGCGACACCCCAAACAATAAACCTTGATACGGGTTCTTGCCCAGCGCCAGCCAACAGCAAAGGAACATGTCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCCACCCCGTACCGCCAATACGCACAAAGGCGCTAGGGCACACATTACCCCAGAACGCGATTGGCCCCCAAATCACCACCACGGAGACCa                                                                                                                                                           >  2:726457/1‑239 (MQ=255)
                                                                                     aGCGTCACCCCAATCAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCAAAAACGCCCAAAGGCGCGAGGGCCAACGTGACCCCAGAACGCAATTGGCCCCACAATACCCACCATGGTGACCAGCg                                                                                                                                                        >  2:672937/1‑238 (MQ=255)
                                                                                       cGTCACCCCAATCAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCCCCATGGTGACCAGCGGcc                                                                                                                                                     >  1:551014/1‑239 (MQ=255)
                                                                                       cGTCACCCCAATCAATAAACCTTGATACGATTTCTTGCCCAGCGCCAGCCCACCGCAAAGGAACATGGCCGCCGAGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGAAGAGCGCTAAAACCTAAAATATAACCCAACCCCGTACCGCCAATACGCTCAAAGGCGCGAGGGCCAACGTTACCCCAACACAGAATCGGCCCCATACTCACCACCATCGGGCCCAACGGcc                                                                                                                                                     >  2:708445/1‑239 (MQ=255)
                                                                                           ccccccCATCAAAAAACCTTGACTCGGTTTCCTGCCCACCGCCACCCACCCGCAAAGGAACATGGCCCCCGCCCCACAGGCTAACAGAAGCGGTAACGAGATTAACTCCTGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAgg                                                                                                                                                  <  2:932062/237‑1 (MQ=255)
                                                                                             ccccAATCAAAAACCCTGTGACGGTTTTCTTACCCAGGCCCCGCCAACCAAAAATGAAAAGGGCCGCCCCGCACCAGACAAAAAGCAGGGGTAACGAAATGAACTCCAGTGGCAGAGGGATAAGACCTAAAATGAAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTg                                                                                                                                                <  2:435536/239‑1 (MQ=255)
                                                                                                 ttgcggttggctggcgctggggTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGCACACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCAAGGAGACCACCGGCCAGGTGCt                                                                                                                                              >  1:913019/21‑238 (MQ=255)
                                                                                                   tcaatAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCACAGGGCCGAGGGACAACGTTACCCCCGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGCGCTTTcc                                                                                                                                          >  1:426759/1‑239 (MQ=255)
                                                                                                    caatAACCCTTGATCCGGTTTCTTGCCAAGCGCCAGCCAACGGCAAAGGAACATGGCCGGCGCGCACAAGACTACACACCGCGGTAACGAGACTAACCCCAGTTGCAGCGCGATAAGACCTAAAATCGAACCCAACCCCGAACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCg                                                                                                                                         <  2:87734/239‑1 (MQ=255)
                                                                                                    caaaAACCCTTGAGTCGGTTCCTTGCCCAGCGCCAGCCAACCGCAAAGGAACAGGCCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGCCCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCg                                                                                                                                         <  1:594894/239‑1 (MQ=255)
                                                                                                     aaaaaaCCTTTGGTACGTTTTCTGTCCCACGGCCAGCCAACCGCAAAGGACCATGGCCGCCCCGCACCCGACTAACAACAGAGGTAACAAGATTAACTCCAGTTCCGAAGGAAGAAGACCTAAAATCGAACCCAACACCGAACCGCCAATACGCCCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGTCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTcct                                                                                                                                         <  2:54077/236‑2 (MQ=255)
                                                                                                       tAACCCTTGATCCGGTTTCTTGCCCAGCGCCAGCCACCCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGa                                                                                                                                      <  1:472952/239‑1 (MQ=255)
                                                                                                       tAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCCTAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCgggct                                                                                                                                     >  1:287610/1‑237 (MQ=255)
                                                                                                         aaccgtggATCCGGTTTCTTGCCCAGCGCCAGCCAACCGAAAAGGAACATGGCCGCCGCGCGCAAGACAAACTACGGCGGTAACGAGATTAACCCCAGTTGGAGAGCGATAAGCCCTAAAACCGACCCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCGCCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGAt                                                                                                                                     <  1:347647/232‑1 (MQ=255)
                                                                                                                gATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGATAGTAAAc                                                                                                                              >  1:580918/1‑239 (MQ=255)
                                                                                                                gATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAACAACCACCATGGTGACAAACGGCCAAGAGCTGTCCGGGAAAGAAAAc                                                                                                                              >  2:542586/1‑239 (MQ=255)
                                                                                                                 actcGGTTTCTGGCCCCACGCCAGCCAACCGAAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGATAGTAAACa                                                                                                                             <  1:923033/236‑1 (MQ=255)
                                                                                                                 aTACGGTTTCTTGCCCAGGGCCAGCCAACCGAAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAAGCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGATAGTAAACa                                                                                                                             <  1:370048/239‑1 (MQ=255)
                                                                                                                        tttCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCAAGAACGAGATTGGCCCCATAATAACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGTTAGTAAACAGGCGGa                                                                                                                       >  1:230969/1‑239 (MQ=255)
                                                                                                                          tCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCCCCAGACTAACATCAGCGGTAACGAGATTAACCCCAGTTGCAGGGCGATAAGACCTAAAATCGAACCCAACCCCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCTATAATCACCACCATGGTGACCAGCGGc                                                                                                                                                      <  1:37294/207‑1 (MQ=255)
                                                                                                                              gcccgggcccagccaACCCAAAAGGAACCGGGCCCCCCCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGACCCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGATAGTAAACAGGCGGATAATga                                                                                                                 <  1:808921/231‑1 (MQ=255)
                                                                                                                                  aGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATcacca                                                                                                                                                                       >  1:41331/1‑181 (MQ=255)
                                                                                                                                    cgcCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCgg                                                                                                                                        >  1:535931/1‑210 (MQ=255)
                                                                                                                                    cgcCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGGGCTTTCCGGGAAAGTAAACAGGCGGATAATGAGAAAAg                                                                                                           >  1:456030/1‑239 (MQ=255)
                                                                                                                                              ccGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCACTACGCCCAAAGGAGCGAGGGAAAACGTTACCCCAGAACGAGATTAGCCCCATAATCACCACCACGGTGACCAGCGGCCAGGGGCTTTCCGGGTGAAAAAACAAGCGGAACATGAAGAAAGAGGGGAGGaa                                                                                                 >  2:782761/1‑239 (MQ=255)
                                                                                                                                                        catgcgccgcccccccacAGATTAACATGATCGGTAACGAGATTAACCCCTTTTGCAGTCGGATAAAGCATAAAATCAAACCCAACACCAGGCCGCACATACGCTCAAAGGCGCGAGGGAAAACTTTACCCCATAACGAGATTGGCCCCATAATCACCACCATGATGACCAGCGGCCAGGTGCTTTCCGGGATAGTAAACAGGCGGATAATGAGAAAAGTGAGCAGGAATGCCAgcgc                                                                                        <  2:943501/221‑1 (MQ=255)
                                                                                                                                                         cctGGCCCCCGCACACCAAACTAACATCAGCTGTAACGAGATTAACTCCAGTTGCCGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGATAGTAAACAGGCGGATAATGAGAAAAGTGAGCAGGAATGCCAGCGCGa                                                                                      <  1:391180/237‑1 (MQ=255)
                                                                                                                                                               cgccgcGCACCAGCCTAACACCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAACCGAACCCAACCCCGTACCGGCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGATAGTAAACAGGCGGATAATGAGAAAAGTGAGCAGGAATGCCAGCGCGACCCGa                                                                                 <  1:908750/238‑1 (MQ=255)
                                                                                                                                                               cgccgcGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGCCCTAAAATCGAACCCAACCCCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGATAGTAAACAGGCGGATAATGAGAAAAGTGAGCAGGAATGCCAGCGCGACCCGa                                                                                 <  1:727047/238‑1 (MQ=255)
                                                                                                                                                                 ccgcGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGTTGACCAGCGGCCAGGTGCTTTCCGGGATAGTAAACAGGCGGATAATGGGAAAAGTGAGGAGGAATGCCAGCGCGACCCGAGt                                                                               >  1:1046999/1‑238 (MQ=255)
                                                                                                                                                                 ccgcGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGATAGTAAACAAGCGGATAATGCGAAAAGTGAGCAGGAATGCCAGCGCGACCCGAGt                                                                               >  1:917424/1‑238 (MQ=255)
                                                                                                                                                                           gtcTAACATCAGCGGCAACGAATATAACCCCAGTGGCAGGGTGATAAGCACAAAAACCGCCCCCACACCTGACCGCCAATACTCCAAAAGGGGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGATAGTAAACAGGCGGATAATGAGAAAAGTGAGCAGGAATGCCAGCGCGACCCGAGTACCATGCACAAt                                                                   <  2:287610/237‑1 (MQ=255)
                                                                                                                                                                                  cATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGGGCTTTCCGGGATAGTAAACAGGCGGATAATGAGAAAAGTGAGCAGGAATGCCAGCGCGACCCGAGTACCATGCCCAATGCGGt                                                              >  1:1044668/1‑239 (MQ=255)
                                                                                                                                                                                  cATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGCCCAGGTGCTTTCCGGGATCGTAAACAGGCGGATAATGAGAAAAGGGAGCAGGAATGCCAGCGCGACCCCGGGACCAGGCACaaaccggg                                                              >  2:781309/1‑233 (MQ=255)
                                                                                                                                                                                       gCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGa                                                                                                                                      >  2:951097/1‑162 (MQ=255)
                                                                                                                                                                                       gCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGa                                                                                                                                      <  1:951097/162‑1 (MQ=255)
                                                                                                                                                                                       gCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGATAGTAAACAGGCGGATAATGAGAAAAGTGAGCAGGAATGCCAGCCCGACCCGAGTACCATGCACAATGCGGGAATgg                                                         >  1:857238/1‑239 (MQ=255)
                                                                                                                                                                                       gCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACCTTACCCCAGACCGAGATTGGCCCCCTAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGATACTAAACAGGCGGGTAATGGGAAAAGTGAGCAGGAATGCCAACACCACCCGAGTACCATGCCCaaagagggaaag                                                          >  2:280631/1‑228 (MQ=255)
                                                                                                                                                                                         tggaaCGAGATTAACTCAAGTTGCAGAGCGATAAGCCCTAAAATCGAACCCAACCACTGCCCGCCAACACGGTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGATAGTAAACAGGCGGATAATGAGAAAAGTGAGCAGGAATGCCAGCGCGACCCGAGTACCATGCACAATGCGGTAATg                                                          <  1:726457/233‑1 (MQ=255)
                                                                                                                                                                                             aCGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTc                                                                                                                                           >  1:551739/1‑151 (MQ=255)
                                                                                                                                                                                                     aaCCCCAGTTGCAGGGCGATAACACCTAAATCCGACCCCAACCCCGTACCGCAAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGATAGTAAACAGGCGGATAATGAGAAAAGTGAGCAGGAATGCCAGCGCGACCCGAGTACCATGCACAATGCGGTAATGGCGGTATACCCGGAt                                          <  1:781309/239‑1 (MQ=255)
                                                                                                                                                                                                        tCCAGTTCCAGAGGAATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGATAGTAAACAGGCGGATAATGAGAAAAGTGAGCAGGAATGCCAGCGCGACCCGAGTACCATGca                                                                       <  1:594701/208‑1 (MQ=255)
                                                                                                                                                                                                            gTTGCAGAGCGCTTAGACCTAAAGTCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCTCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGCGACCAGCGGCCCGGAGCTTTCCGGACTAGTAAAAGGGCGGATAATGAGAAAAGTAAGCGGGCATGCCAACAGGAACACAGGACCATGCAAAACGCGGCACTGgggggacacaagggtgtacaga                                    >  2:807410/1‑221 (MQ=255)
                                                                                                                                                                                                                                   agcGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGATAGTAAACAGGCGGATAATGAGAAAAGTGAGCAGGAATGCCAGCGCGACCCGAGTACCATGCACAATGCGGTAATGGCGGTATACCCGGATTTCAAACGGGCTTAATGACTTATCGgcacg            <  1:87328/237‑1 (MQ=255)
                                                                                                                                                                                                                                   aTCGAACCCAACCCCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGATAGTAAACAGGCGGATAATGAGAAAAGTGAGCAGGAATGCCAGCGCGACCCGAGTACCATGCACAATGCGGTAATGGCGGTATACCCGGATTTCAAACGGGCTTAATGACTTATCGgcacg            <  1:1003146/239‑1 (MQ=255)
                                                                                                                                                                                                                                       aaCCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGATAGTAAACAGGCGGATAATGAGAAAAGTGAGCAGGAATGCCAGCGCGACCCGAGTACCATGCACAATGCGGTAATGGCGGTATACACGGATTTCAAACGGGCTTATACACTTATCGgcacgcac         >  1:312005/1‑238 (MQ=255)
                                                                                                                                                                                                                                       aaCCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGATAGTAAACAGGCGGATAATGAGAAAAGTGAGCAGGAATGCCAGCGCGACCCGAGTACCATGCACAATGCGGTAATGGCGGAATAACCGGATTTCAAACGGGCTTAATGACTTTTCGgcaccacc         >  1:750706/1‑234 (MQ=255)
                                                                                                                                                                                                                                       aaCCCAACACCGCACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGATAGTAAACAGGCGGATAATGAGAAAAGGGAGCAGGAATGCCAGCCCGACCCGAGTTCCCTGCCCAATGCGGCAATGGCGGAACACCCCGACTCTAAAAGGGGCTTAATGAcacttcggcgaaaa          >  2:117051/1‑225 (MQ=255)
                                                                                                                                                                                                                                              cacCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGAGCTTTCCGGGATAGTAAACAGGCGGATAATGAGAACAGTGAGCAGGAATGCCAGCGCGACCCGAGTACCATGCACAATGCGGCAATGGCGGGATACCACGATTTTAAACCGGCTTAAGGACCTATCGGCACGCACGCCGaaa  >  1:239077/1‑238 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              
TGCCAGTAAAGAGCCGAGGATCACATCGCCGCTTCGCCATAACGCCGTATCAATTTCACCTGTCGGGGAACCCACAACAATTGCCAGCGTCACCCCAATCAATAAACCTTGATACGGTTTCTTGCCCAGCGCCAGCCAACCGCAAAGGAACATGGCCGCCGCGCACCAGACTAACATCAGCGGTAACGAGATTAACTCCAGTTGCAGAGCGATAAGACCTAAAATCGAACCCAACACCGTACCGCCAATACGCTCAAAGGCGCGAGGGACAACGTTACCCCAGAACGAGATTGGCCCCATAATCACCACCATGGTGACCAGCGGCCAGGTGCTTTCCGGGATAGTAAACAGGCGGATAATGAGAAAAGTGAGCAGGAATGCCAGCGCGACCCGAGTACCATGCACAATGCGGTAATGGCGGTATACCCGGATTTCAAACGGGCTTAATGACTTATCGGCACGCACGCAGAAA  >  NC_000913/2080127‑2080598

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