Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,870,464 C→T V423M (GTG→ATG)  dgoT ← D‑galactonate:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,870,4640CT98.7% 216.8 / NA 75V423M (GTG→ATG) dgoTD‑galactonate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (36/37);  minor base G (0/1);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

AGCTCGCTGGAGGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGTTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGCAAAACCGTTACCGAAGAATGCCAGCGCCATCAGGCAC  >  NC_000913/3870217‑3870708
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
aGCTCGCTGGATGGGATTAGCTATGCGCTGGATAGCACAGCCATTCTGCGTTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATCCCGCTTACTGGCGCGCCACGCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCTCACCGGCTTTAAGAAAGTGGATTCTGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACACCGCCCATcag                                                                                                                                                                                                                                                   <  2:447912/251‑1 (MQ=255)
    cGCTGGAGGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGTTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACAGCGCCCATCAGCAgg                                                                                                                                                                                                                                               <  1:418240/251‑1 (MQ=255)
     gCTGGAGGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGTTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGTTGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAgg                                                                                                                                                                                                                                               >  1:122105/1‑250 (MQ=255)
     gCTGGAGGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGTTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGGGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAgg                                                                                                                                                                                                                                               >  2:220805/1‑250 (MQ=255)
        ggAGGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGTTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGatat                                                                                                                                                                                                                                           <  1:120512/251‑1 (MQ=255)
        ggAGGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGTTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGACTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGAGCTTCACATCGCCCATCAGCAGGttat                                                                                                                                                                                                                                           >  2:25043/1‑251 (MQ=255)
         gAGGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGTTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGatat                                                                                                                                                                                                                                           >  1:50681/1‑250 (MQ=255)
              gATTAGCTATGCGCTGGATAGCACAGACATTCTGCGTTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAag                                                                                                                                                                                                                                        >  1:527982/1‑248 (MQ=255)
                 tAGCTATGCGCTGGATAGCACAGACATTCTGCGTTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGAAGGATATAAGAGAgcg                                                                                                                                                                                                                                  >  2:7183/1‑251 (MQ=255)
                   gcgaTGCGCTGGATAGCACAGACATGCTGCGTTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACAGCGCCCATCAGCAGGATATAAGAGAgcgcg                                                                                                                                                                                                                                <  2:508161/248‑1 (MQ=255)
                      aTGCGCTGTATAGCACAGACTTGCGTCGTTCGGTGTGACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACGCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGGGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACCTCGCCCATCAGCAGGATATAAGAGAGCGCGCCa                                                                                                                                                                                                                             <  2:303226/251‑1 (MQ=255)
                               aTAGCACATACCTTCTGCGTTTGGTGTTACTGTACCAGCCAGAGAAGCTCAATTCCGCTTACTGGCGCGCCGCCCGCAAACTGCGTTTAAAACTCTGCGGCAATATAGTAAAAAGCCGATGCCTTACCGCACCGGCTATAAGAAAGTGGATTATGCGCAGTGGGGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAAcgcg                                                                                                                                                                                                                    <  2:162225/251‑1 (MQ=255)
                                 aGCACAGCCATTCTGCGTTCGGTGTTAGTGAAGCAGCCAGAGAAGCTCAATACCGATTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATAGAAGAGAGCGCGCCAATCAACGCGAc                                                                                                                                                                                                                  <  2:386003/251‑1 (MQ=255)
                                       gACATTCTGCGTTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGAGGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGACGga                                                                                                                                                                                                            <  2:122105/251‑1 (MQ=255)
                                          aTTCTGCTTTCGGTGTTACTGAAGCAGCCAAATAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGat                                                                                                                                                                                                         <  1:498303/251‑1 (MQ=255)
                                            tCTGCGTTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATTCCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCAGCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGata                                                                                                                                                                                                        <  2:51932/250‑1 (MQ=255)
                                            tCTGCGTTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGata                                                                                                                                                                                                        <  1:324753/250‑1 (MQ=255)
                                                               aaGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGaa                                                                                                                                                                                    <  1:87324/251‑1 (MQ=255)
                                                               aaGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGa                                                                                                                                                                                     >  1:191715/1‑250 (MQ=255)
                                                                   agcCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACCCTGCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGAGATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCg                                                                                                                                                                                <  2:310656/251‑1 (MQ=255)
                                                                      cAGAGAAGTTCAATACCGCTTACTGGAGCGCCACCCGCAAACTGCGTGTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAATTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTaa                                                                                                                                                                             <  2:433330/251‑1 (MQ=255)
                                                                          gaAGCTCAATACCGCTTAATGGCGCTCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGACGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCACCGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCt                                                                                                                                                                         >  2:93679/1‑251 (MQ=255)
                                                                                    aCCGCTTACTGGCGCGCCACCCGCAATCTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTTACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTTTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTTCGCCAGGt                                                                                                                                                                <  1:372794/250‑1 (MQ=255)
                                                                                    aCCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGt                                                                                                                                                                <  1:561689/250‑1 (MQ=255)
                                                                                      cGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGCAACCCTCCGCCAGGAAcc                                                                                                                                                             >  2:191407/1‑251 (MQ=255)
                                                                                          tACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCcca                                                                                                                                                          <  1:204520/250‑1 (MQ=255)
                                                                                                cgcgcCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAg                                                                                                                                                   <  1:491226/251‑1 (MQ=255)
                                                                                                 gcgcCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGc                                                                                                                                                  >  1:103406/1‑251 (MQ=255)
                                                                                                  cgcCACCCGCAAAATGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGGTATAAACCAGTGCCGGTGCGAAAACGTAACCCTGCGCCAGGTACCCCACCACCAGCg                                                                                                                                                 >  2:158475/1‑251 (MQ=255)
                                                                                                   gcCACCCGCAAACTGCGTTTCACACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCCCATCGCCCATCAGCAGGATATAAGATAGCGCGCCAATCAACGCGACGCCGGAGATATAAACCAGTGCGGGTGCGAAACGGTAACCCAGCGCCAGGGACCCCAACACCAGccg                                                                                                                                                >  2:422048/1‑249 (MQ=255)
                                                                                                      aCCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAAc                                                                                                                                             <  1:7183/251‑1 (MQ=255)
                                                                                                         cGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGt                                                                                                                                          >  1:69929/1‑251 (MQ=255)
                                                                                                              actccGTTTAAAACTCTCCGGCAAGAGAGTAAAAAGCCGATGCCTGCCCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGTACATATTTCCCAGTTGAAGGATTAGCCAACTCGCCTCCCATCGCCCATCAGCAGGATATAAGAGAGCGCCCCAATCAACGCGACGTCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATgc                                                                                                                                     <  1:317460/247‑1 (MQ=255)
                                                                                                              aCTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATgc                                                                                                                                     >  1:70200/1‑251 (MQ=255)
                                                                                                                  cGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGAGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGACCGGTGATgccgcc                                                                                                                                 >  1:35813/1‑251 (MQ=255)
                                                                                                                   gTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGccc                                                                                                                                >  1:500954/1‑251 (MQ=255)
                                                                                                                     ttAAAACTCTTCGGCCAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGAGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCCATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGGGATGCCGccc                                                                                                                                >  2:529044/1‑249 (MQ=255)
                                                                                                                          aCTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGaccac                                                                                                                         <  2:500954/251‑1 (MQ=255)
                                                                                                                           ctctTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGACTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATAAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGGGATATAAACCAGTGCAGGTGCGCAACCGTAACCCTGCGCCAGTTACCCCACCACCAGCAGCACGGTGATACCGCCAAGaccccc                                                                                                                        >  2:461140/1‑248 (MQ=255)
                                                                                                                                   aaaGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGGGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGt                                                                                                                <  2:224509/250‑1 (MQ=255)
                                                                                                                                             aaaaGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACAcgcc                                                                                                      >  1:241028/1‑251 (MQ=255)
                                                                                                                                                   cGATGCCGGACCGCACCGGCTTAAAGAAAGTTGATTATGCGCAGTGGTGTGGCAAATATTCCCAAGTTGAAGGATTAGCCAACGCGCTTCCCATCGCCCATCAGCAGGATAAAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCTCCGGCGAAGTTAAACACGCCGCCGGt                                                                                                <  2:548743/251‑1 (MQ=255)
                                                                                                                                                    taTGCCTGACCGCAGCGGCTTTAAGAAAGTGGATTATGCGCAGCGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGAGCTACACATCGCCCATAAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGCCGGAGATATAAACCCAAGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGCACGGTGATGCCGCCAAGACCACCGGCCCAGTTGACCACGCCGCCggct                                                                                               >  2:55869/2‑249 (MQ=255)
                                                                                                                                                    gATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGtt                                                                                               >  1:466527/1‑251 (MQ=255)
                                                                                                                                                        ccTGTCCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCTACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGTCGAAGTTAAACACGCCGCCGGTTAAAc                                                                                           <  1:377235/251‑1 (MQ=255)
                                                                                                                                                                 aCCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTATGGCCCATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACACTGCGCCAGGTACCCCACCACCAGc                                                                                                                                                  >  2:79819/1‑187 (MQ=255)
                                                                                                                                                                 aCCGGCTTGAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATAGTCCACAGTTGAAGGATTAGCCAACGCGCGTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGc                                                                                                                                                  <  1:79819/187‑1 (MQ=255)
                                                                                                                                                                    ggCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGCAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGAGGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGATAAACCAATCAGgcgc                                                                                >  1:398069/1‑250 (MQ=255)
                                                                                                                                                                    ggCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATCCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGgcgc                                                                                >  1:245886/1‑250 (MQ=255)
                                                                                                                                                                      cTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATc                                                                             >  1:307652/1‑251 (MQ=255)
                                                                                                                                                                         tAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCACCGGTTAAACCAATCCGGCGCATCggg                                                                          >  1:64420/1‑250 (MQ=255)
                                                                                                                                                                            gAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGcc                                                                       >  1:111370/1‑251 (MQ=255)
                                                                                                                                                                              aaGTGGATTATGCGCGGTGGGTTGTCACATATTCCGCAGGTGAAGGATTAGCCAACGCGCTTCACGTCGCCCATCAGCAAGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCCGCGCCAGGTACCCCGCCACCAGCGGAACGGTGATGCCGCCCAGACCGCCGGCGAAGTTAAACACGCCGCCGGGTAAACCAATCATGCGTAGCGGTGCCag                                                                     <  2:35813/251‑1 (MQ=255)
                                                                                                                                                                              aaGTGGATTATGCGCAGTGGTGTGGCCCATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCTGTGCCag                                                                     <  2:318740/251‑1 (MQ=255)
                                                                                                                                                                                gTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGtt                                                                                               <  1:411449/223‑1 (MQ=255)
                                                                                                                                                                                   gATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAgaga                                                                  >  1:108230/1‑249 (MQ=255)
                                                                                                                                                                                       aTGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGGGCCAGAGAAGAAAc                                                            >  1:293780/1‑251 (MQ=255)
                                                                                                                                                                                           gcAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGTCGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGACGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGc                                                        >  2:222801/1‑251 (MQ=255)
                                                                                                                                                                                            cAGTGGTTTGGCCCATATTCCACAGTTGAAGGTTTCGCCAACGCGCTTCACATCGCCCATCAGCCGGATATAAGAGAGCGCGCCAATCAACGAGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCCCCACCCGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAGACCAATCAGGCGCATCGGAGCCAGAGAAAAAACCAg                                                         >  2:339456/1‑249 (MQ=255)
                                                                                                                                                                                            cAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAACCCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAACCCAGCg                                                       >  2:219712/1‑251 (MQ=255)
                                                                                                                                                                                            cAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACACGCTTCACATAAACCATCAGCAGGATATAAGACAGCGCGCCAATCAACGCGACGGCGGAGATATAAACAAGTGCAGGTGCGAACCCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGACCGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGCGAAGACACCAGCg                                                       >  1:233326/1‑251 (MQ=255)
                                                                                                                                                                                                     ggCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGACGAAACCAGCGACCAGGTa                                               >  1:43960/1‑250 (MQ=255)
                                                                                                                                                                                                       cacaTATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATc                                            <  1:175158/251‑1 (MQ=255)
                                                                                                                                                                                                               ccACAGTTGAAGGATTAGCCAACGCGCGTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGCaa                                     <  1:430669/250‑1 (MQ=255)
                                                                                                                                                                                                                 acaGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGCAAAAc                                  <  1:222801/251‑1 (MQ=255)
                                                                                                                                                                                                                        aaGGATTAGCCCAAGCGCTTCACCGCGCCCATCAGCAGGATAAAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGGGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGCAAAACCGTTAcc                           <  2:541120/251‑1 (MQ=255)
                                                                                                                                                                                                                           gcttAGCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGCGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGCAAAACCGTTACCga                         >  1:360144/3‑250 (MQ=255)
                                                                                                                                                                                                                            aTTAGCCAACGCGCTTCACATCGCCCATCAGCAGGATAGAAGAGAGCGCGCCAATCAACGCGGCGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTCCCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGCAAAACCGTTACCgaag                       <  2:318582/251‑1 (MQ=255)
                                                                                                                                                                                                                                gCCAACGCGCTTCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGCAAAACCGTTACCGAAGAATg                   <  1:529503/251‑1 (MQ=255)
                                                                                                                                                                                                                                           tCACATCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGCAAAACCGTTACCGAAGAATGCCAGCGCCa          <  1:25043/249‑1 (MQ=255)
                                                                                                                                                                                                                                            cacaTCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGTCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGCAAAACCGTTACCGAAGAATGCCAGCGCCATCa       <  1:525112/251‑1 (MQ=255)
                                                                                                                                                                                                                                             acaTCGCCCATCAGCAGGATATAAGAAGGCGCGCCAATCAACGCTACGGCTGAGATATAAACCAGTGCAGGTGCGAAACCGGAACCCGGCGCCAGGTACCCCACCGCCAGCGGAACGGTGATGCCGCCCAGGCCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGCCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGCAAAACCGTTACCGAAGAATGCCATCGCCATCAg      <  2:241028/251‑1 (MQ=255)
                                                                                                                                                                                                                                               aacGCCCACCAGCAGGATATAAGAGAGCGCGCCAAGCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCCCCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGCAAAACCGTTACCGAAGAATGCCAGCGCCATCAGGc    <  2:103406/249‑1 (MQ=255)
                                                                                                                                                                                                                                               aTCGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGAATCGGTGCCAGAGAAGAGACCAGCTACCAGGTAATCGAAGCAAAACCGTTACCGAAGAATGCCAGCGCCATCAGGc    >  1:757/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 cGCCCATCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGCAAAACCGTTACCGAAGAATGCCAGCGCCATCAGGcac  <  1:220805/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
AGCTCGCTGGAGGGGATTAGCTATGCGCTGGATAGCACAGACATTCTGCGTTCGGTGTTACTGAAGCAGCCAGAGAAGCTCAATACCGCTTACTGGCGCGCCACCCGCAAACTGCGTTTAAAACTCTTCGGCAAGATAGTAAAAAGCCGATGCCTGACCGCACCGGCTTTAAGAAAGTGGATTATGCGCAGTGGTGTGGCACATATTCCACAGTTGAAGGATTAGCCAACGCGCTTCACATCGCCCACCAGCAGGATATAAGAGAGCGCGCCAATCAACGCGACGGCGGAGATATAAACCAGTGCAGGTGCGAAACCGTAACCCTGCGCCAGGTACCCCACCACCAGCGGAACGGTGATGCCGCCCAGACCACCGGCGAAGTTAAACACGCCGCCGGTTAAACCAATCAGGCGCATCGGTGCCAGAGAAGAAACCAGCGACCAGGTAATCGAAGCAAAACCGTTACCGAAGAATGCCAGCGCCATCAGGCAC  >  NC_000913/3870217‑3870708

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