Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 73,001 A→G G507G (GGT→GGC)  thiP ← thiamine ABC transporter membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091373,0010AG97.3% 231.6 / ‑1.8 75G507G (GGT→GGC) thiPthiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (38/35);  total (39/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 = 1.50e-01

ATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGG  >  NC_000913/72755‑73243
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
aTAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGt                                                                                                                                                                                                                                                  >  1:30549/1‑251 (MQ=255)
   gTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCt                                                                                                                                                                                                                                               <  1:486554/251‑1 (MQ=255)
       gCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTTCTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCt                                                                                                                                                                                                                                           <  2:653924/251‑1 (MQ=255)
             cGCTGGCTGGCGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGAt                                                                                                                                                                                                                                      >  1:67130/1‑250 (MQ=255)
                gggctggcGGCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTGCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGTCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCTGTCCTGGCTGCGATAGGa                                                                                                                                                                                                                                  <  1:156338/250‑1 (MQ=255)
                 ggctggcGTCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAg                                                                                                                                                                                                                                 <  1:73713/251‑1 (MQ=255)
                         tCAGAAAACCGGCGATCAAGTTCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTTCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCGCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAAtt                                                                                                                                                                                                                           <  1:309548/249‑1 (MQ=255)
                            gAAAACCGGCTATCAAGTTCAGCAGGGTACTTG‑ACCCGGCGCCGCTTGGCCCGAGGATCGCTACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGTTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTg                                                                                                                                                                                                                      <  2:26406/251‑1 (MQ=255)
                                aCCGGCGAGCAAGTTACGCAGGGTACTTTTCCCTG‑CGCCGCTTGGCCCGAGGATCTCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAAGTTGCTGGTag                                                                                                                                                                                                                  <  2:93357/251‑1 (MQ=255)
                                      gATCAAGTTCAGAAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCAACTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACAGCGCCGTACTGGCTGCGTTAGGAGCCAATTTGCTGGTAGAGATa                                                                                                                                                                                                             >  1:219478/1‑250 (MQ=255)
                                           agttCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGa                                                                                                                                                                                                                                       >  1:365887/3‑219 (MQ=255)
                                           agttCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGa                                                                                                                                                                                                                                       <  2:365887/217‑1 (MQ=255)
                                            gttCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACACGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGCGATAAAACGGc                                                                                                                                                                                                      >  1:232299/2‑251 (MQ=255)
                                             ttCAGCGGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGGGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGc                                                                                                                                                                                                      >  2:600197/1‑250 (MQ=255)
                                              tCAGCAGGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGc                                                                                                                                                                                                      >  1:628535/1‑249 (MQ=255)
                                                 gcagGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGt                                                                                                                                                                                                 <  1:274904/251‑1 (MQ=255)
                                                 gcagGGTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGt                                                                                                                                                                                                 <  1:425865/251‑1 (MQ=255)
                                                    gggTACTTTTACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCg                                                                                                                                                                                              <  2:30549/251‑1 (MQ=255)
                                                          ttttACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGt                                                                                                                                                                                                                    <  1:344441/223‑1 (MQ=255)
                                                          ttttACCTG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGGGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGt                                                                                                                                                                                                                    >  2:344441/1‑223 (MQ=255)
                                                          ttttACCTG‑CGCCGCTTGGCACGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGAGCGGAAatc                                                                                                                                                                                        >  1:564100/1‑251 (MQ=255)
                                                                    cgccgcTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTAcc                                                                                                                                                                               >  1:142364/1‑251 (MQ=255)
                                                                      ccgcTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTAAGCTAACACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCg                                                                                                                                                                              >  1:429390/1‑250 (MQ=255)
                                                                       cgcTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCTACCGTTACGCTAAAACGCATCGGAAAATGTTGGAAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGAGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATCAAACGGCAGGGTGCGGAAATCATCGTTACCGa                                                                                                                                                                             >  1:609750/1‑250 (MQ=255)
                                                                                                  cgccgcGTTCTACCGGTAAGCTAAAACGCAACGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTGTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTGGCTGGTAGAGATAAAACGGCAGGGTGCGGAAAGCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑Ca                                                                                                                                                 <  2:322304/249‑1 (MQ=255)
                                                                                                   gccgcGTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTCGCAGCAGAATTAACGCGGTGACCGCGACGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGAGGAAATCATCGTTACCGAAC‑AACGCCACCCCACCAAAATCAC‑AAATc                                                                                                                                              >  1:453709/1‑251 (MQ=255)
                                                                                                       cgTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTAGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAAGCAC‑CAATCGAc                                                                                                                                           >  1:413406/1‑250 (MQ=255)
                                                                                                        gTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCACCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGTAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCTGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAg                                                                                                                                         <  1:42736/251‑1 (MQ=255)
                                                                                                        gTTCTACCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCCCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAg                                                                                                                                         <  2:447914/251‑1 (MQ=255)
                                                                                                             aCCGTTAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCa                                                                                                                                       <  1:621366/248‑1 (MQ=255)
                                                                                                                 ttAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCCCCACGCCAAAATCAC‑CAATCCACAGCACGCaggc                                                                                                                                >  2:295334/1‑248 (MQ=255)
                                                                                                                 ttAAGCTAAAACGCATCGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCaggc                                                                                                                                >  1:301157/1‑248 (MQ=255)
                                                                                                                      cTAAAACGCATCGGCCAATGTTGGTAAATCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAACAGAATTAACGCGGTGACCGCGCCGTCCTGGCCGCGAAAGGAGCCAATTTGCTGGTAGAGCTAAAACGGAAGGGTGCGGCAATCATCGTTACCTACCCAAC‑CCACCCCGCCAAAAT‑ACACAAACGACAACACGCA‑GCaaaaa                                                                                                                           >  2:356494/1‑248 (MQ=255)
                                                                                                                              cATCGGCAAATGTTGGTACAGCCAGGTGATATCAGTCCGTTTTAACCTCTCGACCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCCAGCAAAAACCAGCGcc                                                                                                                   >  2:324342/1‑251 (MQ=255)
                                                                                                                                 cGGCAAATGTTTGTAAATCCAGGTGAAATCAGGCAGTTGTAACATCGCGCCCCGGTAGTTTTTCAAGCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAAGTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACAACGCCAAAATCCC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTgc                                                                                                                <  2:653675/251‑1 (MQ=255)
                                                                                                                                 cGGCAAATGTTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTg                                                                                                                 <  1:602166/250‑1 (MQ=255)
                                                                                                                                               aaaGCCAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGttt                                                                                                  <  1:220108/251‑1 (MQ=255)
                                                                                                                                                  gCCAGGTGATATCAGTCAGTTGTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAg                                                                                               <  1:65405/251‑1 (MQ=255)
                                                                                                                                                    cAGGTGATATCAGTCAGTTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAACCAGCGCCTGCGCCAGTGGACGTTTCCggg                                                                                             >  2:577812/1‑251 (MQ=255)
                                                                                                                                                                 tcagtTTTAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTcgac                                                                                >  1:560012/1‑248 (MQ=255)
                                                                                                                                                                        tAACATCTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGCCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGCAATCATCGTTACCGAAC‑CACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTACAGGGCGCGCAGCTCGACCACtttt                                                                         >  2:251620/1‑251 (MQ=255)
                                                                                                                                                                            atctCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCCCCCGTGGACGTTTCCGGGCGCGCAACTCGACCACTTTTacgc                                                                     >  2:334310/4‑248 (MQ=255)
                                                                                                                                                                                cGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCgcg                                                                  <  1:366371/250‑3 (MQ=255)
                                                                                                                                                                                    ccGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCg                                                                                                                                                                                              >  1:177041/1‑124 (MQ=255)
                                                                                                                                                                                    ccGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCg                                                                                                                                                                                              <  2:177041/124‑1 (MQ=255)
                                                                                                                                                                                    ccGGTAGTTTTGAAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCTGTGACCTCGCCTTCCTGGCTTCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGTGCGCGCAGCTCGACCACTTTTAAGCGCGACCAg                                                             <  2:699840/251‑1 (MQ=255)
                                                                                                                                                                                     cGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGcc                                                                                                             <  2:611658/202‑1 (MQ=255)
                                                                                                                                                                                     cGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGAGACCGCGCCGTCCTGGCTGCGATAGGAGCCACTTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGcc                                                                                                             >  1:611658/1‑202 (MQ=255)
                                                                                                                                                                                         aGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCtt                                                         <  1:467399/250‑1 (MQ=255)
                                                                                                                                                                                              ttCAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAAtgc                                                   >  1:617722/1‑249 (MQ=255)
                                                                                                                                                                                                cAATCACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAAtgcc                                                  <  1:577812/250‑4 (MQ=255)
                                                                                                                                                                                                    cACGGTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCg                                             >  1:595239/1‑251 (MQ=255)
                                                                                                                                                                                                      cGGTAAACAGCAGAAAACAGAGGAGCAGCAGAATGAACGCGTTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCTCGACCAGCCTTCAATGCCCAGCGAc                                           <  2:564100/251‑1 (MQ=255)
                                                                                                                                                                                                       ggTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACt                                          <  1:140136/251‑1 (MQ=255)
                                                                                                                                                                                                        gTAAACAGCAGAAAACAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTg                                         >  1:578797/1‑251 (MQ=255)
                                                                                                                                                                                                                  gAAAACAGAGTAGCAGCAGAATTAACGCGATGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAACTCAC‑CAATCGGCAGAACGCAAGCAAAAGCCAGCGCATGCGCCAGCGGAAGATGCAGGACGCGCAGCTCGACCACTTTTAAACGCGACCAGCCTTCAATGCCCAGCGACTAACACAACAt                                >  1:564957/1‑250 (MQ=255)
                                                                                                                                                                                                                    aaaCAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACt                                                                            >  2:187371/1‑204 (MQ=255)
                                                                                                                                                                                                                    aaaCAGAGTAGCAGCAGAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACt                                                                            <  1:187371/204‑1 (MQ=255)
                                                                                                                                                                                                                               cagcagAATTAACGCGGTGACCGCGCCGGCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTGACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAgcgggcgg                  <  2:628535/251‑1 (MQ=255)
                                                                                                                                                                                                                                 gcagAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCCATGCCCAGCg                                             >  1:164445/1‑222 (MQ=255)
                                                                                                                                                                                                                                 gcagAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAgcgggcgatg                >  1:10244/1‑248 (MQ=255)
                                                                                                                                                                                                                                 gcagAATTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTg                <  1:295334/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cagAATTAACGCGGTGACCGCGCAGTCCCGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCCGTGGACGTTTCAGGGCGCGCAGCTAGACCACCTTTAAGCGCGACCAGCCTTCACTGCCCAGCGACTGACACACCATACTGTCGCGGGCGGTg                >  2:289090/1‑250 (MQ=255)
                                                                                                                                                                                                                                     aaTTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGatat            >  1:491849/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aTTAACGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATc           <  1:177651/251‑1 (MQ=255)
                                                                                                                                                                                                                                         aaCGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGTACGCAGGCAAAAGCCAGCGCCTgc                                                                                                                <  2:114247/147‑1 (MQ=255)
                                                                                                                                                                                                                                         aaCGCGGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTgc                                                                                                                >  1:114247/1‑147 (MQ=255)
                                                                                                                                                                                                                                           cgcgGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGGGATATCgcgc       >  1:400685/1‑250 (MQ=255)
                                                                                                                                                                                                                                           cgcgGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAACTCGACCACTTTTAAGCGCGACCAGCCTTCACTGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCgcgcc      >  2:288535/1‑250 (MQ=255)
                                                                                                                                                                                                                                            gcgGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCCGCCCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCAt     >  2:574164/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gcgGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCAt     >  1:167365/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gggTGACCGCGCCGTCCGGGGTGCGATAGGAGCCAATTTGCTGGGAGAGATAAAACGGCAGGGGGCGGAAATCATCGTTACCGAAA‑AACGCCACCCCGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCCCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATc    <  2:17259/250‑1 (MQ=255)
                                                                                                                                                                                                                                             cgGTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCgcgcctc    >  1:263808/1‑248 (MQ=255)
                                                                                                                                                                                                                                               gTGACCGCGCGGTACTGGATGCTATAGGATCCAAATTGCTGGTAGAGATAAACCGGCAGGGTGCGTAAATCAACGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGTACGTTTCAGGGCGCGCAGCTCGACCTCTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACATCATCCTGTAGCGGGCGGTGATATCGCGCATCgg  <  1:243370/251‑1 (MQ=255)
                                                                                                                                                                                                                                               gTGACCGCGCCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCgg  <  1:233736/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
ATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGATCAAATTCAGCAGGGTACTTTTACCCG‑CGCCGCTTGGCCCGAGGATCGCCACCTGCTCGCCGCGTTCCACCGTTAAGCTAAAACGCATCGGCAAATGGTGGTAAAGCCAGGTGATATCAGTCAGTTTTAACATTTCGCCCCGGTAGTTTTTCAATCACGGTAAACAGCAGAAAACAGAGCAGCAGCAGAATTAACGCGGTGACCGCACCGTCCTGGCTGCGATAGGAGCCAATTTGCTGGTAGAGATAAAACGGCAGGGTGCGGAAATCATCGTTACCGAAC‑AACGCCACCACGCCAAAATCAC‑CAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGG  >  NC_000913/72755‑73243

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