Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,205,882 C→T W17* (TGG→TAG)  purD ← phosphoribosylamine‑‑glycine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,205,8820CT97.3% 242.0 / ‑5.8 74W17* (TGG→TAG) purDphosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (46/26);  minor base A (0/2);  total (46/28)
Fisher's exact test for biased strand distribution p-value = 1.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAGCGGCGCTTCCGGGCCGACGATGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTT  >  NC_000913/4205638‑4206127
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
cGCACCTGCGGTGGGGCCGGAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAGCGGCGCTTCCGGGCCGACGATGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCa                                                                                                                                                                                                                                                 <  2:595107/251‑1 (MQ=255)
 gCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAGCGGCGCTTCCGGGCCGACGATGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGcc                                                                                                                                                                                                                                                  >  2:610433/1‑249 (MQ=255)
 gCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAGAGGCGATTCCGGGCCGACGATGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGCGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGcc                                                                                                                                                                                                                                                  >  2:582736/1‑249 (MQ=255)
   aCCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAGCGGCGCTTCCGGGCCGACGCTGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGAGCCGCTTTCTAGGCCAgcg                                                                                                                                                                                                                                              >  1:645099/1‑251 (MQ=255)
         gtttGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCCGCGGCGCTTCCGGGCCGACGATGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCt                                                                                                                                                                                                                                        <  2:252182/249‑1 (MQ=255)
         ggTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCTCCAGCGGCGCTTCCGGGCCGACGATGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTGCTAGGCCAGCGCGTGCt                                                                                                                                                                                                                                        <  2:342470/251‑1 (MQ=255)
         ggTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAGCGGCGCTTCCGGGCCGACGATGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCt                                                                                                                                                                                                                                        >  2:193033/1‑251 (MQ=255)
             gggCCGAAGTTTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTGCACCAGCGGCGCTTCCGGGCCGACGATGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCGTTCTAGGCCAGCGCGTGCTcgc                                                                                                                                                                                                                                     <  1:582736/250‑1 (MQ=255)
               gCCGAAGATTTTCAGCCCGTCGGCGGGGAAGGTATCGACCCCGCCTTTCACCAGAGGCGCTTCCGGGCCGACGATGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGcgcc                                                                                                                                                                                                                                  <  2:156229/251‑1 (MQ=255)
               gCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCCCGCCTTTCACCAGCGGCGCTTCCGGGCCGACGATGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGcgcc                                                                                                                                                                                                                                  <  1:166128/251‑1 (MQ=255)
               gCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAGCGGCGCTTCCGGGCCGACGATGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGcgcc                                                                                                                                                                                                                                  <  1:477684/251‑1 (MQ=255)
                ccGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAGCGGCGCTTCCGGGCCGACGATGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGcgccc                                                                                                                                                                                                                                 >  2:545769/1‑251 (MQ=255)
                   acgATTTTCAGCCCGGCTGCGCTGAAGGTATCGACCACGCCTTTCACCAGCGGCGCTTCCGGGCCGCCGAGGGTCAGATCAATCTTTTCGTTTTGGGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTGGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTGCTATGCCAGCGCGTGCTCGcgcccgcc                                                                                                                                                                                                                              <  2:459255/249‑1 (MQ=255)
                   aaGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAGCGGCGCTTCCGGGCCGACGATGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGTGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGcgcccgcc                                                                                                                                                                                                                              >  1:408950/1‑251 (MQ=255)
                                cggcggCGCGGAAGGTATCGACCACGCCTTTCACCAGCGGCGCTTCCGGGCCGACGATGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAACACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACt                                                                                                                                                                                                                 >  2:244965/1‑251 (MQ=255)
                                cggcggCGCGGAAGGTATCGACCACGCCTTTCACCAGCGGCGCTTCCGGGCCGACGATGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGTCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACt                                                                                                                                                                                                                 >  2:37438/1‑251 (MQ=255)
                                 ggcggcGCGGAAGGTATCGACCACGCCTTTCACCAGCGGCGCTTCCGGGCCGACGATGGGCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCAAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTAGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACt                                                                                                                                                                                                                 >  2:487135/1‑250 (MQ=255)
                                   gggcgcgagaGGTATCAACCACTCTTTCCACCGGGGGCCCTTCCGGGCCGACGATGGTCAGATCATTCTTTTCTTTTTGTCCAAAATCAAACAGCGCTGGGATATCGGTCCCGCCAAGAGCACCGTTTTTCAGCGCGGTTTCCCGTGCAGTCCCTGCTTTCCCCGGAGCAACAAAAACAGTCTCATCCAGCGGCTCCTGGGTCTCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAAt                                                                                                                                                                                                              <  1:654166/242‑1 (MQ=255)
                                           aatgtcgcGACCACCCTTTTCACCAGCGGGCTCTCCGGCCCGAGGATGTTCGGATCAATGTTTTCGTTTTGGGCGAAAACCAACGGCGCCGGTTTTGCGGTCACGCAAATAGCACCGTTTTTCAGCGCGGGTCCCAGTGCAGGGCCTGCATTCCCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGCCTGGGCCGCTTTCTAGGCCCGCGCGTGCTCGCCCCCGCCGTTACCAATCACTAATACTTTCAt                                                                                                                                                                                                      <  1:404645/244‑1 (MQ=255)
                                             ggtcgcgcccaccctttccccgaGCGGCGTTTCCGGGCCGACGATGGCCAGACCACTCTTTTCGTTTTTGGCGAAATCAAACAGCGCCGGGATATCGGTCACGCCAATAGCAGCGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGAAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTCGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCt                                                                                                                                                                                                    <  2:24755/229‑1 (MQ=255)
                                             ggTATCGACCACGCCTTTCACCAGCGGCGCTTCCGGGCCGACGATGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACAAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCACCCGCCGTTACCAATCACTAATACTTTCATCt                                                                                                                                                                                                    >  2:670814/1‑251 (MQ=255)
                                                aTCGACCACGCCTTTCACCAGCGGCGCTTCCGGGCCGACGATGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTAtt                                                                                                                                                                                                 >  2:594956/1‑251 (MQ=255)
                                                         gCCTTTCACCAGCGGCGCTTCCGGGCCGACGATGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTAGCGCCCGCCGTTACCAATCACTAATAc                                                                                                                                                                                                            >  2:642731/1‑231 (MQ=255)
                                                         cccTTTCCCCAGCGGCTCTCCCGGGCCGACTATGGCCAGATCAATCTTTTCGTTTTGTGAGAAATCCACCGGCGCGGGGATATCGGCCACGCCAATAGCAACGTTTTGCAGCGCGGGTCCCAGTGCAGTGCCTGCATTCCCCGGAGCAACAAAAACAGTCTCACCCAGCGGCGCCTGGGCCGCTTTCTAGGCCAGCGCTTGCCCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTa                                                                                                                                                                                        <  2:245840/250‑1 (MQ=255)
                                                             ttccccAGCGGCGCTTCCGGGCGGACGATGGTCAGATCAATCTTTTCGTTTGGGGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGGGCAGTGCCTGCATTACCCGGAGAAACAAAAACATTCTCAACCAGCGGCGACTGGGCGGCTTTGAAGCCGAGCGCGTTCTGGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAAt                                                                                                                                                                                      <  1:668068/245‑1 (MQ=255)
                                                                cccAGGGGCGCCTCCGGGCCTCCGAGGGTCATACAAATCTTTTTGTTTTGTGCAAAACCCACCAGCGCGGGGATGCCGGTCACGCCAATAGAACCGTTTGGCAGCGCGGTTCCCAGTGCGGGGCCTGCTTTCCCCGGAGAAACAAAAACAGTCCCACCCCGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCTGGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCgg                                                                                                                                                                                 <  2:1852/250‑1 (MQ=255)
                                                                  cAGCGGCGCTTCCGGGCAGACGATGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATCGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGa                                                                                                                                                                                >  2:372193/1‑250 (MQ=255)
                                                                   aGCGGCGCTTCCGGGCCGACGATGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGa                                                                                                                                                                                >  2:227110/1‑249 (MQ=255)
                                                                          cTTCCGGGCCGACGATGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGCCGGAAGTGgcgc                                                                                                                                                                        >  1:364446/1‑250 (MQ=255)
                                                                             gccGGGCCGACGAGGGGCAGATAAATCTTTCCGTTTTGTGCGAAACCCAACAGCGCCGGGATCTCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCACCCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCAt                                                                                                                                                                      <  1:364962/249‑1 (MQ=255)
                                                                                gggccgccgaggttcaaatcatttttttCGTTTGGTGCGAAACCCAACAGCGCGGTGATATCGGTCCCGCCAATAGCACCGTTTTTCAGCGCGGGTTCCAGTGCAGTGCCTGCATTCCCCGGAGCACCAAAAACAGTCTCACCCAGCGGCGCCTGGGCCGTTTTCTAGCCCAGCGCGTCCTCGCGCCCGCCGTTACCAATCACTAATACGTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCg                                                                                                                                                                  <  1:157779/227‑1 (MQ=255)
                                                                                 ggCCGACGATGGTCAGATAAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCACTAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGGGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCgg                                                                                                                                                                 >  1:37431/1‑251 (MQ=255)
                                                                                  gCCGACGAGGGTCAGATCAATCTTTTCTTTTGGTGCGAAATCCACCGGCGCCGGGATTTCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTCCCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGt                                                                                                                                                                <  1:541585/251‑1 (MQ=255)
                                                                                          gggtaggatcatttttctctttttGTGCAAAACCCACCAGCGCGGGTATTTGGGTCACGCCAATAGCAACTTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCTTTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGTCGGGGCCTTTTTCTAGCCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCa                                                                                                                                                        <  2:46305/231‑1 (MQ=255)
                                                                                                 aTCAATCTTTGCGTTTTGTGCGAAATCCAACAGCGCCGGGATCTCGGTCATGCCAATAGCAACGTTTTGACGCGCGGGCTCCAGTGCAGTGCCTGCATTCCCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAAt                                                                                                                                                 >  2:110015/1‑251 (MQ=255)
                                                                                                     attttttCTTTTTGTGGGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCGGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCg                                                                                                                                             <  1:483054/248‑1 (MQ=255)
                                                                                                          ttCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGCGCATCGCAATACCGTGCtc                                                                                                                                        >  1:638590/1‑251 (MQ=255)
                                                                                                                    gcaaaagccAACGGCGCCGGGATTTCGGCCACGCCAATAGCAACGTTTTGCAGCGGGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAACCAGTCACAACCAGCGGCGCCTGGGCCGCTTTCTAGGCCAGCGCGTGCCCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTcggcggcg                                                                                                                              <  1:16380/244‑1 (MQ=255)
                                                                                                                      gAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCggcggcggc                                                                                                                            >  1:127585/1‑251 (MQ=255)
                                                                                                                            caccAGCGCCGGGATATCGGTCGCGCCAATAGCAGCGTTTTGCAGCGCGGGTTCCAGTGCAGTCCCTGCATTCCCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGTCAATCAc                                                                                                                      <  1:183795/248‑1 (MQ=255)
                                                                                                                                  cgcCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAACAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCACTCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGCGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGAGGCGGCAATCACTTAGtc                                                                                                                >  1:319864/1‑251 (MQ=255)
                                                                                                                                   gcCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGtc                                                                                                                >  1:577332/1‑250 (MQ=255)
                                                                                                                                    ccGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGtcat                                                                                                              >  1:590380/1‑251 (MQ=255)
                                                                                                                                       ggATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGCAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCAc                                                                                                           >  1:368856/1‑251 (MQ=255)
                                                                                                                                         atatCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACgg                                                                                                         >  1:663800/1‑251 (MQ=255)
                                                                                                                                           agcGGTCACGCAAATAGCAACGTTTTGCAGCGCGGTTTCCAGTGCAGTGCCTGCATTACCCGGAGCAAAAAAAACAGTCTAAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGAt                                                                                                       <  1:610433/249‑1 (MQ=255)
                                                                                                                                                    gCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGACGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGcc                                                                                              >  1:387699/1‑251 (MQ=255)
                                                                                                                                                            cAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTgg                                                                                       >  1:496742/1‑250 (MQ=255)
                                                                                                                                                               cgttttgaggctcgGGTTCCAGTGCGGGCCCTGATTTCCCGAGGCAAAAAAAACCAGTCCCAACAGGCGGCGCTGGGCCTGGTTTCTAGGCCAGCGCTGGCCCCCCCCCGCCGTTACCAATCACTAAGACTTTCAGCTATGGCTCCATTAATGGCGGAAGGGGAGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGGGGCGGCAATCCCTGCGTCACCACGGATATAACCGCCAGGCTGGAtt                                                                                    <  2:233132/238‑1 (MQ=255)
                                                                                                                                                                 ttttGCAGCGCGGGTTCCAGGGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCAGGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGGCAgatggattaagc                                                                                 >  2:551475/1‑240 (MQ=255)
                                                                                                                                                                    tGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGCTAGAACCGCCAGGCTGGATTACGc                                                                                >  1:191913/1‑249 (MQ=255)
                                                                                                                                                                     gCAGCGCGGGTTCAGTGGCAGGGCCGGCATTCCCCGGGACAACAAAAACAGTCTCAACCAGCGGCGACGGGGCCGCTTTCAAGCCCACCGGGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAAGGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAg                                                                              <  1:655082/250‑1 (MQ=255)
                                                                                                                                                                         cgcggttcccAGTACAGTCCCTGCATCCCCCGGAGCAAAAAAAACAGTCTCAACCAGCGGCGCCGGGGCCGCTTTCTAGGCCCGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAGGTCAc                                                                         <  1:383228/243‑1 (MQ=255)
                                                                                                                                                                         cgcgGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCCGGCTGGATTACGCAGGTccc                                                                         >  2:141551/1‑249 (MQ=255)
                                                                                                                                                                            gggTTCCAGTGCAGTGCCACCATTACCCGGGGCAACAAAAACATTCTCAACCGGCGCCCGCGGGACCGTTTTCTAGCCCAGCGCGTGCCCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGGCGGTGAAGATCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCCCTGCGTCATCCCGGAGAGAACCGCCAGGCTGGATTACGCAGGTCACGcc                                                                      <  1:580909/251‑1 (MQ=255)
                                                                                                                                                                            gggTTCCAGTGAAGTCCCTGAATTCCCCGGAGCAAAAAAAACAGTCCCAACCAGGGGCGACGGGGCCTCTTTCTAGGCCGGCGCGTGCCCGCCCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAGGTCACGcc                                                                      <  1:206481/251‑1 (MQ=255)
                                                                                                                                                                               ttCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACAAGAGGCGACTCGGCCGCTTTCTAGGACAGCGAGTACTCGAGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGACGCATGTCGGTAAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCCATCACTTCGTCATCACGGATAGAACCGACAGGCTGGATTACGCAGGACACGCCCGc                                                                   >  1:92179/1‑251 (MQ=255)
                                                                                                                                                                                  cagtgaagtgCCTGCATTACCCGGAGCAACAAAAAAAGTCTCAACCAGCGGCGACTGGGCAGCTTTCTAGGCCAGCGCGAGCTCGCGCCCGCCGTTACCAATCCCTACTACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCCTGTCGGTGAAGAGCACCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCCCGGCTAGAACCGCCAGGCTGGATTACGCCGGTCAAGCCCCCAgc                                                                >  2:425415/1‑251 (MQ=255)
                                                                                                                                                                                            gcCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCCCTTCGTCATCACGGATAGAACCGCCCGGCTGGATTACGCAGGTCACGCCCGCAGCggcggcggc                                                       >  1:24802/1‑251 (MQ=255)
                                                                                                                                                                                             ccTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCggcggcggcc                                                      >  2:261602/1‑250 (MQ=255)
                                                                                                                                                                                                gCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGACAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGATGGCTTACGCAGGGCACGCCCGCAGCGGCGGCGGCATc                                                    >  2:446640/1‑250 (MQ=255)
                                                                                                                                                                                                   ttACCCGGAGCAACAAACACAGGCTCCACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATAACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGCGGCGCATGTCGGTGACGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCCTCACGGATAGAACCGCCCAGCTGGGTTAAGCAAGTCACGCCCCCCGCGGAGGCGGCAtgcatc                                                >  2:325538/1‑246 (MQ=255)
                                                                                                                                                                                                    tACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCCGCGGCGGCGGCCTCAATAc                                               >  1:628643/1‑251 (MQ=255)
                                                                                                                                                                                                      cccGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCg                                             >  1:394085/1‑251 (MQ=255)
                                                                                                                                                                                                                  aaaaCAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACggg                                 >  1:597209/1‑251 (MQ=255)
                                                                                                                                                                                                                     aacggTCTAAACCAGGGGCGCCTGGGCCGCTTTTTAGGCCGGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGgaa                               <  1:545769/247‑1 (MQ=255)
                                                                                                                                                                                                                      cAGTCTCAACCAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCACCCGCCGTTACCAATCCCTAATACTTTCATCTATTGCTCCATTACTGGCGGAAGTGGCGCATGCCGGTGAAGAGCATCGCCATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACAGGTAGAACCGCCAGGATGGATTACGCCGGTCACGCCCGCAGAGGCGGCGGCCTCAATACCGACGCGGAACGGga                                >  1:707095/1‑248 (MQ=255)
                                                                                                                                                                                                                                ccAGCGGCGACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAAc                           >  1:44939/1‑244 (MQ=255)
                                                                                                                                                                                                                                     ggcgACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAgg                                                                             >  1:221539/1‑189 (MQ=255)
                                                                                                                                                                                                                                        gACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCa                                                                                             >  1:568/1‑170 (MQ=255)
                                                                                                                                                                                                                                        gACTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATc            >  1:254781/1‑251 (MQ=255)
                                                                                                                                                                                                                                         aCTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCCGGCTGGATTACGCCGGTCACGCCCGCCGCGGCGGCGGCCTCCCTACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCg           >  1:245945/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cTGGGCCGCTTTCTAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGATGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGACATCg           >  1:614808/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   cTTTCTAGGCCAGCGCGTACTCGCGCGCGCCGTTACCTATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAATCGCGCATGTCGGAGAAAAGCCTCGCAATAACGTGATCGTGGGCGGCAGCAATCCCTTCGTCATCACGGATAGAACCGCCCCGGTGGATTAAGCCAGTCACGCCCGAAGCGGACGTGGCATCAAGACCGTCAAGGAACGGGGAGAACGAGTCAAAAGcactcgaaggacctt  >  2:672314/1‑237 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
CGCACCTGCGGTTGGGCCGAAGATTTTCAGCCCGGCGGCGCGGAAGGTATCGACCACGCCTTTCACCAGCGGCGCTTCCGGGCCGACGATGGTCAGATCAATCTTTTCGTTTTGTGCGAAATCCAACAGCGCCGGGATATCGGTCACGCCAATAGCAACGTTTTGCAGCGCGGGTTCCAGTGCAGTGCCTGCATTACCCGGAGCAACAAAAACAGTCTCAACCAGCGGCGACTGGGCCGCTTTCCAGGCCAGCGCGTGCTCGCGCCCGCCGTTACCAATCACTAATACTTTCATCTATTGCTCCATTAATGGCGGAAGTGGCGCATGTCGGTGAAGAGCATCGCAATACCGTGCTCGTCGGCGGCGGCAATCACTTCGTCATCACGGATAGAACCGCCAGGCTGGATTACGCAGGTCACGCCCGCAGCGGCGGCGGCATCAATACCGTCGCGGAACGGGAAGAACGCGTCAGAAGCCATCGAGGAACCTT  >  NC_000913/4205638‑4206127

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