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.C96.7% 196.3 / ‑0.7 63intergenic (+586/+56)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (25/33);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-01
*NC_0009134,296,3802.G95.2% 192.1 / 1.5 63intergenic (+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 (26/33);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 5.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

TCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑A‑‑C‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTG  >  NC_000913/4296196‑4296619
                                                                                                                                                                                          ||                                                                                                                                                                                                                                                     
tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGtt                                                                                                                                                                                 >  1:19607/1‑251 (MQ=14)
   tATCAGGTCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGtt                                                                                                                                                                                 <  1:494598/248‑1 (MQ=12)
    aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                                                                                                                                                                             >  1:441210/1‑251 (MQ=18)
    aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                                                                                                                                                                              <  1:463351/250‑1 (MQ=18)
       aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                                                                                                                                                                          >  1:698982/1‑251 (MQ=21)
              cGCCAGACAGCGCAATTGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                                                                                                                                                                   <  1:109735/251‑1 (MQ=255)
                 cAGACAGCGCAATAGCCTGATTTAGCGTACTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTACCCCAGCCAGCGCTATAGCCTGCTTTAGCTTCACTTTCTGGATCGGATCCGGCGTTTCAG‑C‑AG‑AA‑‑‑TCCGACACCATTGCCCGATGGTTCG‑TTGACCGTGTCTTATCAGACCTCTCTAAACCGTAGGTTAATTTTTCCAAcc                                                                                                                                                                >  2:174669/1‑250 (MQ=255)
                      aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                                                                                                                           <  1:304446/251‑1 (MQ=255)
                       gttcaataaatggATTTAGCATCATTTTTGAGATCGGATAAGGTGTTTATCCCGAATCGTTCATAAATGCCTGATTAGGCGTTTATCGTGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTTGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ATGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAATCTCTTcac                                                                                                                                                          <  2:623308/239‑1 (MQ=9)
                               gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                                                                                                                                  <  2:576646/251‑1 (MQ=255)
                                     tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                                                                                                                            >  2:163464/1‑251 (MQ=255)
                                            tGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAACCTCTTCACACCCGGTATCAAACCCTTcc                                                                                                                                      >  2:217545/1‑250 (MQ=255)
                                              aTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCAt                                                                                                                                    >  1:388343/1‑250 (MQ=255)
                                                          ggATAAGGCGTTTCTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGCTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCATACCCTTCCATACAGCTCAGGTTc                                                                                                                       >  1:110810/1‑251 (MQ=255)
                                                          ggATAAGGCGTTTATGCCGCATCCGTCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTc                                                                                                                       <  1:105104/251‑1 (MQ=255)
                                                            aTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCa                                                                                                                      >  1:257372/1‑250 (MQ=255)
                                                              acggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAcc                                                                                                                   <  2:388866/249‑1 (MQ=255)
                                                              aaGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCaccc                                                                                                                   >  2:308270/1‑248 (MQ=255)
                                                               aGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCa                                                                                                                  <  1:101783/251‑1 (MQ=255)
                                                                ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAg                                                                                                                 <  1:4180/251‑1 (MQ=255)
                                                                 gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGt                                                                                                                <  2:612667/251‑1 (MQ=255)
                                                                              aTCCGCCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTACACCAGTGCAGCGCCTTTTg                                                                                                   >  2:656281/1‑251 (MQ=255)
                                                                                 cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCtt                                                                                                <  1:390968/251‑1 (MQ=255)
                                                                                    cATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGtt                                                                                             <  2:115985/251‑1 (MQ=255)
                                                                                    cATCAACGCCTGATGCGACGCTTAACGCGTCNTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGtt                                                                                             <  2:369036/251‑1 (MQ=255)
                                                                                     aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGttc                                                                                            >  2:165693/1‑251 (MQ=255)
                                                                                        aaCGCCTGATGCGACGCTTATCGTGTCTTATCAGTCCTACGCCAGACAGCGCAATAGCTTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑GCCGACAACAATGCCTGATGTGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTTGTTAATTTTTCAAACTCTTCACACCCGGTATCAATCCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGttcttc                                                                                         <  1:715768/251‑1 (MQ=21)
                                                                                         aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTtcttct                                                                                        >  2:268061/1‑251 (MQ=255)
                                                                                          cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTtctc                                                                                       <  2:441210/251‑1 (MQ=255)
                                                                                               gATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTCATTTAGCGTGATTTTGTAGGTCGGATAAGGAGTTT‑A‑‑C‑‑G‑CAGACTCCGACATCAATGCCTGATGCGACGCCTG‑CCGCGTCTTATCCGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTAAGGTTCAACCAATGCAGCGCCTTTTGCTTGTTCTTCTCGATa                                                                                   >  1:650132/1‑250 (MQ=21)
                                                                                               gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATaa                                                                                  <  2:154837/251‑1 (MQ=255)
                                                                                                 gtgtaaccttacagcggcTTATCCGGCCTACGTCAGTCAGCGCAATATCCTGATTTGGCGTGGTTTTGTAGGTCGGATAAGGCGTTT‑ACTC‑CG‑CA‑‑‑CCCGGCATCAATGCCTGATGCGACGCTTT‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACCCCCGGTATCAAACACTTCCATACAGCTCAGCTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTGCGCGATaaaa                                                                                <  2:572090/234‑1 (MQ=255)
                                                                                                     gcgcgtaactcggcTTATCAGGCGCACGCCAGACATCGCAATAGCCGGATTTAGCGTGAGTTTGTAGGACGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACctt                                                                            <  2:110810/238‑1 (MQ=25)
                                                                                                       gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGAGATTTTGTAGGTCGGATAAGGCGTGT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGATTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACcttct                                                                          >  1:327414/1‑251 (MQ=25)
                                                                                                        cTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACcttctt                                                                         >  1:1514/1‑251 (MQ=255)
                                                                                                           aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCa                                                                       <  1:148279/250‑1 (MQ=255)
                                                                                                            aCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCAcc                                                                     <  2:121547/251‑1 (MQ=255)
                                                                                                              gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGt                                                                   <  2:388343/251‑1 (MQ=255)
                                                                                                                gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGttt                                                                 >  1:571252/1‑251 (MQ=255)
                                                                                                                    tattaGGCCTTCGCCAGGTAGCTCTATAGCCTGAGTGAGCGTGTTTTTGTATGTCGGATAAGGCGGTT‑A‑‑CGCGGCA‑‑‑TCCGATATCAATGCCTGATGCTACGCTTT‑CCGCGTCTTATCCGGCTTATCTTAACCGTTGGTTCATGTTTCAAGCTCGTCACACCCGGTAGCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCGTTTGCTTGTTCTTCGCGATAAAGCCTTCTTCACCGTTTa                                                                <  2:508985/244‑1 (MQ=255)
                                                                                                                     aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAgg                                                             <  2:151762/250‑2 (MQ=255)
                                                                                                                     aTCAGGCCTACGACAGACAGCGCAATAGCCTGATTTAGCGTGCTTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCCGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTAACACCCGGTATCAAACCCTTCCATACAGATCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACAGTTTAAGa                                                             >  2:591187/1‑250 (MQ=25)
                                                                                                                        aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAcat                                                         <  2:445089/251‑1 (MQ=255)
                                                                                                                            ccaggCCAGACAGCGCAATAGCCTGATTTAGCGTGATTGTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCatcagcttttcgcaaaacg                                                     <  2:525841/247‑19 (MQ=25)
                                                                                                                              aCGCTAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCAtt                                                     <  1:585767/249‑1 (MQ=255)
                                                                                                                                   aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGccc                                              >  2:294139/1‑251 (MQ=255)
                                                                                                                                   aGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGccc                                              <  1:431201/251‑1 (MQ=255)
                                                                                                                                       cgcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGCTAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGt                                           >  2:140903/2‑250 (MQ=255)
                                                                                                                                       aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGt                                           >  2:534473/1‑250 (MQ=255)
                                                                                                                                          gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTc                                       <  2:698982/251‑1 (MQ=255)
                                                                                                                                           cAATAGCGTGCTTTAGCGTAATTTTGTAGGGNGGATAGGGTGTTT‑ACGC‑CG‑CA‑‑‑CCCGCCATCAATGCCTGATGCGACGCTTG‑CCGAGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCGTGTGCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCg                                      <  2:353876/251‑1 (MQ=255)
                                                                                                                                               aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAAt                                  <  1:308270/251‑1 (MQ=255)
                                                                                                                                                    gatgtAGCGTGAGTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCTACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAAcc                              <  2:257372/246‑1 (MQ=255)
                                                                                                                                                    gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCg                             <  1:406200/251‑1 (MQ=255)
                                                                                                                                                              gATTTTGTAGGTCGGCTAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGaa                    >  2:51144/1‑250 (MQ=255)
                                                                                                                                                              cgtaccgTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTACTTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAAAGCTCAGGTTCAACCCGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAACCCTTCTTCACCGTTTCAGAACATCCTTCCCGCCCAGTACTCGGAATAACCGGTTCGGGaac                   >  2:90380/7‑250 (MQ=38)
                                                                                                                                                               aTTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAAt                  >  1:668392/1‑251 (MQ=255)
                                                                                                                                                                 tttGTAGGTCGGATAAGGCGTTT‑AAGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAAACTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTg                >  2:504875/1‑251 (MQ=255)
                                                                                                                                                                  ttGTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTg                <  1:607008/250‑1 (MQ=255)
                                                                                                                                                                    gTAGGTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCTGTTCGGGAAATTGCag             <  2:408890/251‑1 (MQ=255)
                                                                                                                                                                        gTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATCAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAACTTGCAGAGCt         >  2:411690/1‑251 (MQ=255)
                                                                                                                                                                        gTCGGATAAGGCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCt         >  1:287385/1‑251 (MQ=255)
                                                                                                                                                                               agggCGTTT‑ACGC‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTAATCAGGCCTATCTTAACCGTTGGTTACTTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACCGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTCCTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTg  >  2:482127/3‑251 (MQ=255)
                                                                                                                                                                                          ||                                                                                                                                                                                                                                                     
TCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑A‑‑C‑CG‑CA‑‑‑TCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTG  >  NC_000913/4296196‑4296619

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