Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,380:2 +G 90.8% intergenic (+586/+56) gltP → / ← yjcO glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,3802.G90.8% 137.6 / 6.1 65intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (2/4);  new base G (29/27);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTT  >  NC_000913/4296172‑4296595
                                                                                                                                                                                                                     |                                                                                                                                                                                                                         
aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTCTTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCAACCGACATCAAt                                                                                                                                                                                                        >  1:222606/1‑228 (MQ=2)
                                           gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGGCAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGGGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcg                                                                                                                                                            >  1:390500/1‑228 (MQ=0)
                                           gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                                                                                                                                            >  1:518909/1‑229 (MQ=11)
                                                cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATCAGGCGTTTCCC‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGCCGCTGGCCCGGTCTTATCAGGCCACTCTTACCCGTTgg                                                                                                                                                       >  2:358848/1‑229 (MQ=11)
                                                cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCATGATTTTGTATGCCGGATAAGCCGTTT‑‑ACGC‑‑CGCATCCGACATCAATACCTGATGCGCCGCTTGCCGCTTCTTATCAGGCCTATCTTACCCGTTgg                                                                                                                                                       >  1:674098/1‑229 (MQ=11)
                                                   aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                                                                                                                                                    >  1:84299/1‑229 (MQ=19)
                                                          tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                                                                                                                                              >  1:609124/1‑228 (MQ=255)
                                                             tttAGCGTTATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGGTTAACGAGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                                                                                                                                           <  1:304443/228‑1 (MQ=255)
                                                              ttgactgttttttGTAGGGGGATTAAGGCGTTTGTGCCGACGCCGACATAAAGGCCGGGTGCGCCGCTGAGGGTGTTTTATCAGTCCGACGCCAGACAGCCCAATAGCCGGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                                                                                                                                         <  2:532236/220‑1 (MQ=14)
                                                                    tGATTGTGTAGGTCGGATAAGGCGTTTATGCCGCATCCCACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                                                                                                   <  1:252799/229‑1 (MQ=255)
                                                                        tttGGAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                                                                                               <  1:194917/229‑1 (MQ=255)
                                                                          tGTAGGGCGGATAAGGGGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                                                                                                             <  1:248102/229‑1 (MQ=255)
                                                                                  ggagtaGGCGTTTAAGGCCGATGCGACATCAACGCCTGGTGCGAAGGTTAACGGGTCTTATAAGGCCTACGCCAGACAGCGCAATAGCCTTATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                                                                                                     <  1:299704/224‑1 (MQ=17)
                                                                                  ggATAAGTGGTTTAGGCCGCATCCGACATGAACGCCTGATGCGACGCTTAACTCGTTTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                                                                                                     <  1:701840/229‑1 (MQ=17)
                                                                                  ggATAAGGGGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGGTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                                                                                                     <  1:573335/229‑1 (MQ=255)
                                                                                  ggATAAGGCTTTTAAGCCGCATCCGACATCAACGCCGGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                                                                                                     <  1:620986/229‑1 (MQ=255)
                                                                                  ggATAAGGCGTTTATGCCGCATCCGACATCAACGACTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                                                                                                     <  1:644654/229‑1 (MQ=255)
                                                                                  ggATAAGGCGTTTAGGCCGCATCCGACATCAACGCCTGCTGCGACGCTTAACGCGTCTTATCAGGGCTACGCAAGACAGCGCAATAGCCTGATTTAGCTGTATTTTGTAGGTCGGAGAAGGTGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                                                                                                     <  1:358848/229‑1 (MQ=17)
                                                                                  ggATAAGGCGTTTAAGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa                                                                                                                     <  1:41490/229‑1 (MQ=255)
                                                                                     tAAGGCTTTTTTGCCGCATCCGAACTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCACGTAAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAccc                                                                                                                  <  1:499751/229‑1 (MQ=17)
                                                                                        ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTGGGGTAATTTTTCAAACTTTTCACCCACGGTATCAAACCCTTc                                                                                                               >  2:410712/1‑229 (MQ=14)
                                                                                              tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCCATGCCTGATGCGACGCTTGCAGCGTCTTATCAGGCCTATCTTACCCGTTGGTTAATTTTTCAAACTCTTCCAACCCGGTATCaaa                                                                                                                     >  2:120035/1‑217 (MQ=14)
                                                                                              tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAACAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaag                                                                                                                     <  1:120035/217‑2 (MQ=255)
                                                                                                 tccgctacccgcATCAACGCCTGATGCGTCGCTTAACGCGTTTTAGCAGGACTAAGCCAGATAGCGCAATCGCCTGGTTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTT‑‑ACGC‑‑CGCAGCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCt                                                                                                      <  2:205486/218‑1 (MQ=18)
                                                                                                  ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGTTCGGATAAGGCGTTT‑‑ACGC‑‑CGAATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTCATCAGGCCCATCTTAACCGTTGGTTAATTTTTCAACATCTTCACACCCGGTATCAAAACCTTCCATTCAACTc                                                                                                     >  2:307487/1‑229 (MQ=17)
                                                                                                   cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCTTACAGCTCa                                                                                                    >  1:157095/1‑229 (MQ=255)
                                                                                                    gCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGCTAAGGCGTTT‑‑ACGC‑‑CGCATCCGACACAAATTCCTGATCCGACGATTGCCTCGATCGATCAGCCCTTTCTTAACCTTTGGTTACTTTTTCAAACACTTCACACCCGTAATCAAACCCTTCCATACacccccg                                                                                                   >  2:251989/1‑223 (MQ=14)
                                                                                                    gCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCCAACCCGGTATCAAACCCTTCAATACAGTTCAg                                                                                                   >  1:22226/1‑229 (MQ=255)
                                                                                                    gCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTCCCCCCCCGGTATAAAAACCCTCcacacaaccccg                                                                                                   >  2:748365/1‑223 (MQ=14)
                                                                                                      aTCCTACAGCAACGCCTGATGAGACGCTGAACGGGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt                                                                                                 <  1:748881/229‑1 (MQ=18)
                                                                                                                   gccttacgcagcGCTTGATGTGTTTTTTCAGGTCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAg                                                                                    <  1:733733/218‑1 (MQ=21)
                                                                                                                   gCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAg                                                                                    <  1:78321/229‑1 (MQ=255)
                                                                                                                      tGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCCCACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgcgc                                                                                 >  1:594989/1‑229 (MQ=255)
                                                                                                                        aTGCGACGCTTAGCGCGTTGTATCAGTCCTACGCCAGACAGCGCAATAGCCGGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACTCCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCt                                                                               <  1:348928/229‑1 (MQ=21)
                                                                                                                             aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTg                                                                           >  1:426796/1‑228 (MQ=255)
                                                                                                                                        gTCTTATGAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTTGTTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTTATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTtcttct                                                                <  2:609123/228‑1 (MQ=255)
                                                                                                                                          cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTtctc                                                               >  1:592274/1‑227 (MQ=255)
                                                                                                                                           ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTCTCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCgtt                                                            >  1:399749/1‑227 (MQ=255)
                                                                                                                                            tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATa                                                           <  1:135650/229‑1 (MQ=255)
                                                                                                                                              ttcggCTTACGCCAGACAGCTAAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACCCCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATaaa                                                         <  2:518909/226‑1 (MQ=25)
                                                                                                                                                    cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGCTATCAAACCCTTCCATAAAGGTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACcttc                                                   >  1:269687/1‑229 (MQ=255)
                                                                                                                                                       cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTTTCAAACCCTTCCCTCACGCCTAGGTTCAACCAGTGCACCGCCTTTTGCTTGTTTTTCTCGCTAAAACCttcttc                                                >  2:294451/1‑229 (MQ=21)
                                                                                                                                                            gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGCCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCACACCCTTCCCTGCCGCTCTGGTTCAACCCGTGCACACCCTTTTGCCTGTTCTTCTCGATAAAACCTACTTCAcct                                            >  2:358949/1‑227 (MQ=21)
                                                                                                                                                               aGCGCCATAGCCTGATTTGGCGGTATTTTTTAGGTCGGATAAGGCGTTT‑‑A‑‑CGTCGCATCCTACATCAATGCCTGATGCTACTCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCGTCTTCACCGTTTa                                        <  2:157095/229‑1 (MQ=255)
                                                                                                                                                                 cgcCATAGCCGGATTTGGCGGTGTTTTTGTGGTCGGATAAGGCGTTT‑‑ACGC‑‑CCCATCCGACACCAACTCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACCCCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAg                                      <  2:355016/229‑1 (MQ=37)
                                                                                                                                                                 cgcAATAGTCTGATTTGGCGTGATTTTTGAGGTCGGATAAGGCGTTTTCG‑‑C‑‑CGCAGCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACGCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAg                                      <  2:426796/229‑1 (MQ=37)
                                                                                                                                                                 cgcAATAGCCTGATTTAGCGTGATTTTGTATGTCGGATAATGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGTCTGATGAGACTCTTGCCGAGTCTTAGCAGCACTATCTTACACGCCTGTTAATTTTTCAACCCCGTCACCCCCGGTATCACACCATTCCATACCGCTACGGTTCACACAACGCCACGCCTTTCTCTAGTTCTTCTCCCTCAGACATTCTTCCCCGTTTaa                                       >  2:512632/1‑228 (MQ=9)
                                                                                                                                                                 cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGTTCAGGt                                                                                                 >  1:229319/1‑170 (MQ=255)
                                                                                                                                                                 cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt                                                                                                 <  2:229319/170‑1 (MQ=255)
                                                                                                                                                                 cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAg                                      <  1:431064/229‑1 (MQ=255)
                                                                                                                                                                 cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGCTTaa                                       >  1:280310/1‑228 (MQ=255)
                                                                                                                                                                 cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTAAGGTTCAACCAGGGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTaa                                       >  1:687773/1‑228 (MQ=255)
                                                                                                                                                                 cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAATGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTaa                                       >  1:39973/1‑228 (MQ=255)
                                                                                                                                                                 cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCGTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAACTTTTCAAACTCCTCACACCCGGTATCAAACCCTTCCATACAGCTTAGGTTCAACCAGACCAGCGCCCTTTCCCTGTTCTTCCACGTAAAACCTTCTTCACAGTTTaa                                       >  2:383582/1‑228 (MQ=21)
                                                                                                                                                                      ttgCCTGATTTAGCGTGATTTTGGAGGTCGGATAAGGCGTTT‑‑ACCC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCGTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAcat                                 <  1:533392/227‑1 (MQ=255)
                                                                                                                                                                           tGATTTAGCGTGATTTTGTAGGGCGGGTAAGGCGTTTCCG‑‑C‑‑CGCATCCGACATCAATGCCTGATCCGACGCTTCCCGAGTCTTATCAGGCCTATCTTAACCGTTGGTTACTTTTTCAAACTCTTCACCCACGGTATCAAACCCTTCCATACAGCTCAGTTTCAGCAATCACAGGCCCATCTCCTTGGTCTTTTCGATAAAACCTTCCACACCGATTAAGAACACCAtt                             >  2:712622/1‑228 (MQ=255)
                                                                                                                                                                                tAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACCC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGcc                       <  1:410712/229‑1 (MQ=255)
                                                                                                                                                                                   cGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCCGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAg                    >  2:95873/1‑229 (MQ=255)
                                                                                                                                                                                   cGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATATTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAACGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTCCTTAACCGTTTAAGAACATACTTCcccccccg                    >  2:3296/1‑227 (MQ=37)
                                                                                                                                                                                       ttttttttAGGTGGGGTATGGGGTTT‑‑T‑‑CCGCGCCTCCGACATCAATGCCTGATGCGAAGCTTGCCGCTTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACTCCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTAc                 <  2:84299/227‑1 (MQ=255)
                                                                                                                                                                                            tgtgGGTCGGTTAAGGGGTTTCCG‑‑C‑‑CGCATCCGACCTCAATGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGa            <  2:39973/229‑1 (MQ=255)
                                                                                                                                                                                                 ggcGGATAAGGCTTTT‑‑ACGC‑‑CGCATCCGACATCAATTCCTGATGCCACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAAc       <  1:272117/227‑1 (MQ=255)
                                                                                                                                                                                                 gTCGGATTAGGCGTTT‑‑ACCC‑‑CGCATCCGACAATAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAAc       <  1:336498/229‑1 (MQ=255)
                                                                                                                                                                                                   cGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCg     >  1:463050/1‑229 (MQ=255)
                                                                                                                                                                                                      aTAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGtt  <  1:466047/229‑1 (MQ=255)
                                                                                                                                                                                                                     |                                                                                                                                                                                                                         
ACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑‑CGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTT  >  NC_000913/4296172‑4296595

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