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.C95.8% 214.0 / 0.5 72intergenic (+586/+56)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base C (22/46);  total (25/47)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01
*NC_0009134,296,3802.G95.8% 203.6 / ‑0.0 72intergenic (+586/+56)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base G (23/45);  total (25/47)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

CGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTG  >  NC_000913/4296192‑4296619
                                                                                                                                                                                               ||                                                                                                                                                                                                                                                   
cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGCCGCTTAACGCGTCTTATCAGGCCTACGCCAGCCGGCGCTATAGTTTGCTTTAGCTTGATTTTGTATGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGAAGCTTGCCG‑CGTCTTATCAGGCCTATCTTAAc                                                                                                                                                                                     >  2:508107/1‑251 (MQ=11)
    ttgtaccagtcccaccccAGTCAGCGAAATGGCCTGATGTAGCTTGATTTTGTGGCTTGTACAAGGCGTGTATTCCTCATCCGACATCACTGCCTGATGCGACGCTTAACGCGTCTTATCAGGCTTACGCCAGACAGCGCAGTAGCCTGATTTAGCGTGGTTTTGGAGGGCGGTTTAGGCGTTT‑‑ACGC‑C‑‑TCCTCCGACATGAATGCCTGATGCTACGCTTGCCG‑CGTCTTATCAGTCCTATCTTAACCGt                                                                                                                                                                                  <  2:830024/234‑1 (MQ=2)
        aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGtt                                                                                                                                                                             >  1:404317/1‑251 (MQ=18)
           aggcttccGCCAGTCAGCGCAATAGTCTGAATTATCGTGATTTTGTAGGTCGGATAAGGCTTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                                                                                                                                                                          <  2:615005/244‑1 (MQ=11)
            ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                                                                                                                                                         <  1:829595/251‑1 (MQ=25)
                 aagacagccAGCGCAATAGCCGAATTTAGCGTTATTTTGTAGGACGGGTAAGGGGTGTTTGCCGCATCCGACATCAACGCCGGATGCGACGTTGAACTCGTCTTATCAGGCCTACGCCAGACATCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGAGAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGGCTTATCAGGCCTATCTTAACCTTTGGTTAATTTTTc                                                                                                                                                                     <  2:413786/242‑1 (MQ=11)
                  agtagtgaggcgctatagtcagatgtAGCGTGATTTTGTAGGTCTGATAAGGGGTTTATGCCGCATCAGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                                                                                                                                                                   <  1:417475/226‑1 (MQ=11)
                    ccAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                                                                                                                                                                 <  1:627854/251‑1 (MQ=255)
                          aGCGTAATAGCCGGATTTAGCCTGATTTTGTAGGTTGGATAAGGCGTTTATGCCGCATCCCACATCAATGCCTGATGCGACGCTTAACGCGTCTTATGATGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGTCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                                                                                                                           <  2:534970/251‑1 (MQ=14)
                             gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCacac                                                                                                                                                        <  1:29274/251‑1 (MQ=255)
                               aaTAGCCTGATTTAGCTGGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc                                                                                                                                                      <  1:874423/251‑1 (MQ=255)
                                aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                                                                                                                                     >  1:186330/1‑251 (MQ=255)
                                  aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGt                                                                                                                                                   <  1:162562/251‑1 (MQ=255)
                                      tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑CGCA‑CCATCCGACATCAATGACTGATGCGACGCTTTCCG‑TGTCTTATCAGGCCTATCTTCACCTTTAGTTAATTTTTccacctcagcacacccattaac                                                                                                                                                >  2:776053/1‑229 (MQ=11)
                                       gatgtAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCa                                                                                                                                               <  1:160741/246‑1 (MQ=255)
                                       gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCa                                                                                                                                               >  1:719402/1‑250 (MQ=255)
                                          ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                                                                                                                            <  1:508107/250‑1 (MQ=255)
                                                 gagtttGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                                                                                                                                     <  1:246634/247‑1 (MQ=255)
                                                 gATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACGTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCTTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCAt                                                                                                                                    <  1:531953/251‑1 (MQ=255)
                                                 gATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCAt                                                                                                                                    >  1:720936/1‑251 (MQ=255)
                                                  aTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGAAGCTTGCCG‑CGTCTTATCAGGCCTCTCTTAACCGTTGGTTAATTTTTCAATCTCTTCAAACCCGGTATCAAACCCTTCCAt                                                                                                                                    >  2:323522/1‑250 (MQ=17)
                                                         aGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAATGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCCGGCCTCTCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCAAACAGCTcc                                                                                                                            >  2:682727/1‑250 (MQ=17)
                                                             cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGtt                                                                                                                        <  1:625761/251‑1 (MQ=255)
                                                              ggATAAGACGTTTATGCCGCATCCGACATCAACGCCTGATGCGACACTTAACGCGTCTTATCAGGACAACGCCAGACAGCGCTATAGCCTGTTTTAGCGTGATTTTGTAGGCCGGAAAAGGCGTTT‑‑ATGC‑T‑‑CCATCCGACATCCATGCCTGATCCGACGATTGCCC‑CCCCTTTTCACGACACCCTTAACCGATGGTTACTTTTTCAACCTCGTTAAACCCGGTACCAAACCCATCCATACAGCTCAGGTTc                                                                                                                       >  2:52953/1‑251 (MQ=2)
                                                               gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGCCATCAATGCCTGATGCGACGCTTGCCGACG‑CTTACCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGCTATCAAACCCTTCAATCCAACTAAGGTTc                                                                                                                       >  2:470333/1‑250 (MQ=17)
                                                                  aaGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc                                                                                                                    <  1:631597/250‑1 (MQ=255)
                                                                         gtATGCCCCACCCGACATCAACGCCGGATGCGACTCTTAACCCGTCTTATCGGGCCTACGCCGGACAGCGCAATAGCCTGATTTAGCTTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑CCCGCATCCGACATCAATGCATGCTGCGACGCTTGCCG‑CTTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAg                                                                                                            <  2:832313/250‑1 (MQ=255)
                                                                             gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgcgc                                                                                                         <  1:804101/250‑1 (MQ=255)
                                                                                gCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCAACGCTTGCCG‑CGTCTTATCAGGCCTCTCTTCACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATAAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCttt                                                                                                     >  2:541527/1‑251 (MQ=18)
                                                                                 cATCCGACATCAACTCATGATGCTACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAACCTGAATTAGCGTGATGTTGTAGGTCGGATAAGGCTGTT‑‑ACTC‑C‑‑GCTTCCTACATCAATGCCTGATGCGACGCTTGCCT‑CGTCTTATCAGGTCTATCGTAACCGTTGGCGAATTTTTCAAACTCTTCACACCCGGTATCAAACCTTTCCATACAGCTCAGGTTCAACCAGTGCAGCGACtttt                                                                                                    <  1:837312/251‑1 (MQ=18)
                                                                                 cATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCtttt                                                                                                    >  2:314560/1‑251 (MQ=255)
                                                                                   cccGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGc                                                                                                  <  1:123783/250‑1 (MQ=255)
                                                                                     cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCCGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCACCCAGTGCAGCGCCTTTTGCtt                                                                                                >  2:313785/1‑251 (MQ=255)
                                                                                     cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCtt                                                                                                >  1:748013/1‑251 (MQ=255)
                                                                                         aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCCG‑‑C‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGAAGCGCCTTTTGCTTGttc                                                                                            >  2:121940/1‑251 (MQ=255)
                                                                                           cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGttctt                                                                                          <  1:639545/251‑1 (MQ=255)
                                                                                             aCGCCTGATGCGACGCTTAACGCTGCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTtcttct                                                                                        <  2:375930/251‑1 (MQ=255)
                                                                                              cGCCTGATGCGACGTTTAACGCGTCTTGTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GAATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTAAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTtctc                                                                                       <  1:134193/251‑1 (MQ=21)
                                                                                              cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTtctc                                                                                       >  1:6379/1‑251 (MQ=255)
                                                                                              cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGACTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTtctc                                                                                       <  1:817117/251‑1 (MQ=255)
                                                                                               gCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACATCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCg                                                                                      <  1:158934/251‑1 (MQ=255)
                                                                                                    aTGCGACGGTTAACGCTTCTTATCAGGGCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATGTTGTAGGTCGGATAAGGCGGTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACTCCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATaa                                                                                  <  2:475162/250‑1 (MQ=21)
                                                                                                      gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATaaaa                                                                                <  1:121940/250‑1 (MQ=255)
                                                                                                       cgctgCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATATCCTGATTTAGCGTTATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAAcc                                                                              <  2:77004/247‑1 (MQ=255)
                                                                                                        gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATCGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAACCct                                                                             >  2:424669/1‑251 (MQ=255)
                                                                                                          cGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACcttc                                                                           <  1:640790/251‑1 (MQ=255)
                                                                                                          cGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACcttc                                                                           <  1:464471/251‑1 (MQ=255)
                                                                                                                   cgTCTTATCAGGCCTACGCCAGACAGTCCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTACCCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGtt                                                                  >  1:319572/1‑251 (MQ=255)
                                                                                                                    gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACATTCTTCACCGttt                                                                 >  1:550820/1‑251 (MQ=255)
                                                                                                                        tttcatgCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTaa                                                               <  2:790316/243‑1 (MQ=255)
                                                                                                                           cgggCCTACGCCAGAAATCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAca                                                          <  2:186330/249‑1 (MQ=255)
                                                                                                                           ccggCCTACGCCAGAAAGCGCAATTGTCTGATTTAGCGTTAGTTTGTAGGTCGGATAAGGCGCTT‑‑ACGC‑C‑‑GCATCTGACATTAATGCCTGATGCGACTCTTGCCG‑CGGCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACCCTGCACCCCCGGTATCAATCCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCTCCGTTTAAGAAca                                                          <  2:35773/249‑1 (MQ=255)
                                                                                                                            atgCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCTCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAcat                                                         <  2:308012/249‑1 (MQ=255)
                                                                                                                            aGGCCTACGCCATACAGCGCTATAGCCTGATTTAGCGTGATTTTGGAGGGCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCGTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAcat                                                         <  2:17483/251‑1 (MQ=37)
                                                                                                                             ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAcatc                                                        >  1:22850/1‑251 (MQ=255)
                                                                                                                                  aCGCCAGACAGCGCAATAGCCTGATTTTGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAcatcat                                                      <  2:109087/248‑1 (MQ=255)
                                                                                                                                   cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGTCGTTT‑‑ACGC‑C‑‑GCATCCGCCATCAATGCCTGATGCGCCGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGCATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTccc                                                  >  1:551231/1‑251 (MQ=37)
                                                                                                                                   cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                                                                                                                                                                 >  2:287218/1‑140 (MQ=255)
                                                                                                                                   cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                                                                                                                                                                 <  1:287218/140‑1 (MQ=255)
                                                                                                                                         aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGTAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTGCACCGTTTAAGAACAGCATTCCCGCCCAg                                            <  1:710978/251‑1 (MQ=255)
                                                                                                                                               caagagcatGATTTCGCGTGATTTTGTAGGTAGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTGTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGGCCGGTACTCg                                      <  2:116951/243‑1 (MQ=37)
                                                                                                                                                 aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTCACCGTTGGTTACTTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGCTAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTAc                                         >  2:778637/1‑246 (MQ=255)
                                                                                                                                                 aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACAGCATTCCCGCCCAGTAc                                         <  1:778637/246‑1 (MQ=255)
                                                                                                                                                     ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATaa                                <  1:234430/251‑1 (MQ=255)
                                                                                                                                                        gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAAc                               <  1:224486/249‑1 (MQ=255)
                                                                                                                                                         aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATGCCCGCCCAGTACTCGGAATAACCgg                            <  1:314560/251‑1 (MQ=255)
                                                                                                                                                            tAGCGTGATTTTGGAGGTCGGATAAGTCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGtt                          <  2:404973/250‑1 (MQ=255)
                                                                                                                                                             aGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCg                        <  1:464557/251‑1 (MQ=255)
                                                                                                                                                                     tttGTAGGTCGGCTAAGGCGATT‑‑ACGC‑C‑‑GCATCCGACATCATTGACTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAAACGTTGGTTAATTTTTCAAACTCTTCACACAAGGTATCAAACCCTTCCATAAAGCTCAGGTTCAACCAGTGCAGCGCATTTTGATTGTTCTTCTCGATAAAACATTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACAGGTTCGGGACATTg                >  2:680007/1‑251 (MQ=37)
                                                                                                                                                                      gttgaGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑TCATCCGACTTCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGc               <  2:591282/247‑1 (MQ=255)
                                                                                                                                                                          aGGTCGGATAAGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCaga            <  1:313785/250‑1 (MQ=255)
                                                                                                                                                                                    aGGCGTTT‑‑ACGC‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTg  >  2:628611/1‑250 (MQ=255)
                                                                                                                                                                                               ||                                                                                                                                                                                                                                                   
CGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑C‑‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCG‑CGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTG  >  NC_000913/4296192‑4296619

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