Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,296,381 +GC intergenic (+587/+55) gltP → / ← yjcO glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,3801.C92.9% 160.3 / 4.4 56intergenic (+586/+56)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base C (23/29);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01
*NC_0009134,296,3802.G92.7% 161.7 / 4.5 56intergenic (+586/+56)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (3/1);  new base G (23/28);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

GATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTG  >  NC_000913/4296178‑4296619
                                                                                                                                                                                                           ||                                                                                                                                                                                                                                                   
gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                                                                                                    <  1:333383/251‑1 (MQ=11)
                         aGGCCTCCGCCAGACAGCGCACTGACCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                                                                                                                                                                           <  2:540800/251‑1 (MQ=11)
                         aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                                                                                                                                                                           >  1:68166/1‑251 (MQ=21)
                               ccGCCAGACAGCACAATCTGCTGATTTAGGGTGATTGTGTAGGTCGTATAATGGTTTTATGCCTGATCAGACATCACCGCCTGATGCAACGCGTAACGCGTCGTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCTTGATTTTGTATGTCGGATAAGGCGTTTT‑‑CG‑CAGCATTCGACATCAATGCCTTATGCGACGCTTGTCGCGTCTGATCAGGCCTATCTTAACCTTTTGTTAAttttt                                                                                                                                                                       <  2:319424/248‑1 (MQ=2)
                                     gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct                                                                                                                                                               <  1:575267/251‑1 (MQ=255)
                                      aCAGCGCAATAGCCTGATTTAGCGTGATTTTTTAGGTCGGATAAGGCGTTTATGCCGCATCTGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                                                                                                                              <  2:115270/251‑1 (MQ=255)
                                       cAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                                                                                                                                             >  1:234178/1‑251 (MQ=255)
                                         gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTGTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                                                                                                                                                           <  2:693181/251‑1 (MQ=255)
                                              aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                                                                                                                                      >  1:346471/1‑251 (MQ=255)
                                               tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCTAGCCAGCGCAATAGCCTGATTTCGCGTGACTTTGTACATCGGACACTGCGTTTACGC‑‑C‑TCATCCGACATAGATGCCTGAATCTACGTTGGCCGCGTCTTAACAGCCCTATCTTAACCTTTGGTTAATTTTTCACACTCTTCACCCCCg                                                                                                                                                      >  2:7425/1‑250 (MQ=9)
                                                   cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATc                                                                                                                                                 >  1:598437/1‑251 (MQ=255)
                                                    tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATc                                                                                                                                                 >  1:516335/1‑250 (MQ=255)
                                                     gagttAGCGTGATTTTGTAGGGCGGATAAGGCGTTTATGCCGCATCCGCCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTGTTCAAACTCTTCACACCCGGTATc                                                                                                                                                 <  2:527425/246‑1 (MQ=14)
                                                      aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGGCATGAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCGACGCCAGACAGGGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGGGGTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                                                                                                                              <  1:694645/251‑1 (MQ=17)
                                                      aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaa                                                                                                                                               >  1:673655/1‑250 (MQ=255)
                                                          agactgagtttGGAGGTCGGAGAAGGCGTTTATGCCGCATACGACACCAACGCCGGATGCGACGCTTAACGCGTCTTATCAGGTCTGCGCCAGCCAGCGCAATAGCCTGATTTAGCGTTATTTTGTAGGTCGGATAAGGGGTTTACGC‑‑C‑GCATCGGCCATCAATGCCTGATGCGGCGCTGGCGTCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACGCTTCACACCCGGTATCAAAccc                                                                                                                                           <  1:447599/242‑1 (MQ=255)
                                                            cGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTc                                                                                                                                       <  1:153493/251‑1 (MQ=255)
                                                                      tAGGTCGGATAATGCGTTTATGCCGCAGCCGACATCAACTCTTGATGCGACGCTTACCTCGTCTTATCAGGGCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCt                                                                                                                              <  2:460394/250‑1 (MQ=17)
                                                                           cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTAAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGtt                                                                                                                        >  1:120250/1‑251 (MQ=255)
                                                                            ggATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTc                                                                                                                       <  1:580240/251‑1 (MQ=255)
                                                                                        tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCACGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCCGGTTCAACCAGTGCAgc                                                                                                           >  2:558140/1‑251 (MQ=255)
                                                                                             cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCtt                                                                                                      <  1:78514/251‑1 (MQ=255)
                                                                                              gCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTGATCAGGCCTACGACAGACAGCGCAATAGCCTGATTTAGCGTGATTTAGTAGGTCGGATAAGGCGTTTA‑‑CGAC‑GCATCCAACATCAATGCCGGATGAGACGCTGTACTAGTCTTATCAGGCCTATCTTAACCGTTGGTTTCTTTTTCAAACTCTTCAAACCCGTTA‑CAAACCCACTCCAAACAGCACCGAGTCAACAAGTGCAGCGCAttt                                                                                                     >  2:432935/1‑251 (MQ=9)
                                                                                                  ccGACATCAACTCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTcc                                                                                                                                      <  1:285371/214‑1 (MQ=255)
                                                                                                  ccGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCTCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTcc                                                                                                                                      >  2:285371/1‑214 (MQ=255)
                                                                                                  ccGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACATCTCAGGTTCAACCAGTGCAGCGCCTTTTGCt                                                                                                 >  2:228062/1‑251 (MQ=255)
                                                                                                   cgtcATCAACGCCTGATGCGACGCTTAACGCGTCGTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGGTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCGCAGGTTCAACCAGTGCAGCGCCTTTTGCtt                                                                                                <  2:68166/248‑1 (MQ=255)
                                                                                                   cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCtt                                                                                                <  1:157683/251‑1 (MQ=255)
                                                                                                         cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGttctt                                                                                          >  1:419929/1‑251 (MQ=255)
                                                                                                           aCGCCTGATGCGACGCTTAACGTGTCTTATCCGGCCTACGATAGACAGCGCAATAGCCTGATGTAGCGTGATTTTGTAGGTCGGATAGGGCGTCTACGC‑‑C‑GCACCCGATATCAATGCCTGATGCGACGCGTTCCTCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCA‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCGTGttctt                                                                                          <  2:235733/249‑1 (MQ=255)
                                                                                                             gCCTGATGCGACGCTTCACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTACTTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGcc                                                                                                        >  1:508615/1‑233 (MQ=255)
                                                                                                             gCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTGGGTTAATTTTTCAAACGCGTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGcc                                                                                                        <  2:508615/233‑1 (MQ=255)
                                                                                                               cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGACTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGAGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGAt                                                                                    >  1:454152/1‑251 (MQ=255)
                                                                                                                 gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGAt                                                                                    <  1:73452/249‑1 (MQ=255)
                                                                                                                  aTGCGACTCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGTCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATGTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTGGCTTGTTCTTCTCGATaaa                                                                                 <  1:585498/251‑1 (MQ=21)
                                                                                                                         gcttgccgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGATAATTTTTCAATCTCTTGACACCGGATTTCAGACCA‑TTCCATACTGCTCAGGTTCACCCCAAGCACCGCCTTTTGCTTGTTCTTCTCGATTAAGCcttac                                                                          >  2:247532/7‑249 (MQ=17)
                                                                                                                             aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCAc                                                                      >  1:678580/1‑251 (MQ=255)
                                                                                                                               cgcttCTTATCAGTCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACcttctgccccg                                                                    <  2:346471/247‑7 (MQ=255)
                                                                                                                                gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGt                                                                   <  1:452414/251‑1 (MQ=255)
                                                                                                                                      tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTa                                                                <  1:55169/248‑1 (MQ=255)
                                                                                                                                       aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGa                                                             <  1:24670/250‑1 (MQ=255)
                                                                                                                                         cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCTTTTAAGAAc                                                           <  2:673655/250‑1 (MQ=255)
                                                                                                                                         cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAc                                                           >  2:42560/1‑250 (MQ=255)
                                                                                                                                                aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTcc                                                   <  1:379156/251‑1 (MQ=255)
                                                                                                                                                            gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CGTC‑ACATCCGACATCAAAGTCTGATGCGACCCTTCCCGCGTCTTATCAGGCCTATCTTACCCGTTGGTTAATTTTTCACACTCTTCACACACGGTATCACACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTc                                       >  1:97476/1‑251 (MQ=37)
                                                                                                                                                               aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTACTTTTTCAAACTCTTCACACCCGGTATCAAACCA‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCgg                                     >  1:585052/1‑250 (MQ=255)
                                                                                                                                                                 aGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGACCATCATTCCCGCCCAGTACTCggtat                                  >  2:522527/1‑248 (MQ=255)
                                                                                                                                                                 aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGaa                                   <  1:42560/250‑1 (MQ=255)
                                                                                                                                                                          tAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCAAACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGtt                          <  2:520512/250‑1 (MQ=255)
                                                                                                                                                                                    ttGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGc               >  2:68932/1‑251 (MQ=255)
                                                                                                                                                                                    ttGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGc               >  1:321008/1‑251 (MQ=255)
                                                                                                                                                                                        aggacGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCCACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCag             >  1:64561/5‑249 (MQ=255)
                                                                                                                                                                                         ggTCGGATAAGGCGTTTACGC‑‑C‑GAATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGc          >  1:175229/1‑251 (MQ=255)
                                                                                                                                                                                             ggATAAGGCGTTTACTC‑‑C‑GCATCCGATATCAATGCCTGATGCGACGCTTTCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATGTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCACAGGTTCAACCAGTGCAGCGCCTGTTGCTTGTGCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCg      <  2:414427/251‑1 (MQ=255)
                                                                                                                                                                                             ggATAAGGCGTTTA‑‑CG‑CAGCATCCGACATCAATGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGAGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCCCCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCg      >  2:483722/1‑251 (MQ=255)
                                                                                                                                                                                                 aaGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATTCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCAGTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTg  <  2:532064/251‑1 (MQ=255)
                                                                                                                                                                                                           ||                                                                                                                                                                                                                                                   
GATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCC‑TTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTG  >  NC_000913/4296178‑4296619

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: ‑