Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,380:1 +CG 100% 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,3801.C92.5% 155.9 / 4.7 53intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base C (31/18);  total (32/21)
Fisher's exact test for biased strand distribution p-value = 2.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*NC_0009134,296,3802.G92.5% 157.7 / 4.3 53intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base G (31/18);  total (32/21)
Fisher's exact test for biased strand distribution p-value = 2.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAA‑CCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATT  >  NC_000913/4296209‑4296604
                                                                                                                                                                              ||                                                                                                                                                                                                                                      
aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACACTGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑CTT‑‑A‑‑C‑‑CGCCGACTCCGACCTCACTCCCTTATGCGACACTTGCCCCGCCTTTCCAGGCCCAACTTAAcgccgg                                                                                                                                                                  >  2:33893/1‑234 (MQ=11)
      gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCGGATGCGACGCTTTCCGCGACTTATCAGGCCTATCTTAACCATTGGTTAAt                                                                                                                                                            >  1:647411/1‑239 (MQ=255)
       aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCATTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CCC‑‑A‑TCCGACAGCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACCTTAACCGTTGGTTAAtt                                                                                                                                                           >  1:891551/1‑239 (MQ=255)
           cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGAACTCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                                                                                                                                                       >  1:843885/1‑239 (MQ=255)
              aaTAGCCTGATTTAGGTGGATTTGGTGAGTCGGATAATGGGTTTATGCCGGATCCGACATAAACGCCTGATGCGACGCTTACCGCGTCTTACCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                                                                                                                                                    <  2:176831/239‑1 (MQ=14)
                     tGATTTAGCGTTATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                                                                                                             <  1:194370/239‑1 (MQ=255)
                         ttAGCGTTATTTTGTAGGTCGGATAAGGCTTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc                                                                                                                                         <  1:72691/239‑1 (MQ=255)
                             cGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCCAACTCTTCACACCCGGt                                                                                                                                     >  1:494871/1‑239 (MQ=255)
                                 aTTTTTTAGTGCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGGGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCa                                                                                                                                 <  1:187427/239‑1 (MQ=255)
                                            cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACCCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGTTATCAAACCCTTCCAt                                                                                                                      >  1:738841/1‑239 (MQ=255)
                                                  aGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCt                                                                                                                <  1:748201/239‑1 (MQ=255)
                                                     cGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTc                                                                                                               <  1:44319/237‑1 (MQ=255)
                                                       tttATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTACCCGCTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCCGGt                                                                                                           >  1:8231/1‑239 (MQ=255)
                                                       tttATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACCGCTCAGGt                                                                                                           >  1:186940/1‑239 (MQ=255)
                                                        ttATGCCTCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGt                                                                                                           <  1:33893/238‑1 (MQ=255)
                                                         tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCTCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTACCAAACCCTTCCATACAGCTAAGGTTc                                                                                                         >  1:864231/1‑239 (MQ=255)
                                                         tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATCCAGCCCAGGTTc                                                                                                         >  1:112887/1‑239 (MQ=255)
                                                         tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCAATCTTACCCTTTGCTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACCGCCCAGGTTc                                                                                                         >  1:371597/1‑239 (MQ=17)
                                                         tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAACGCCTGATGCGACGCTTGCCGCGGCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTc                                                                                                         >  1:123330/1‑239 (MQ=255)
                                                              cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGTGATCAAACCCTTCCATCCAGCCCAGGTTCAACCa                                                                                                    >  1:90984/1‑239 (MQ=255)
                                                                    cgtactCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGACTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑CGACGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAg                                                                                              <  1:659595/234‑1 (MQ=255)
                                                                         tCATCGCATGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGGCAGCCAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCt                                                                                         <  1:272303/239‑1 (MQ=21)
                                                                             cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGAAGCGCCtttt                                                                                      >  1:327473/1‑238 (MQ=255)
                                                                              gCCTGATAAGACGCGTTAAGCGTCGCATCAGGCCTACGCCAGACAGGGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATACGACGCTTCCCGCGTCTTATAAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACCCTTCACACCCGGTCTCAAACCCTTCCATCAAGCTCAGGTCCAACCAGTGCCACCCCtttttc                                                                                    >  2:774311/1‑237 (MQ=18)
                                                                                  gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTATCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCTTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAACCTCTTCACACCCGGTACCAAACCCTTCCATACAGCTCAGGTTCTACACGTGCACCACCCTTTGCTTGt                                                                                >  1:533469/1‑239 (MQ=18)
                                                                                         cGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCGGAGTCTTATCAGGCCTATCTTAACCGTTGGTAAATTTTTACAACTCTTCACACCCGGATTCAAACCCTTCCAGACAACCCAGTTTCAACCAGTGCAGCGCCTTTTTCTTGttttcttc                                                                         >  2:762095/1‑233 (MQ=18)
                                                                                         cGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCTTACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTtcttct                                                                          >  1:286042/1‑238 (MQ=255)
                                                                                               aagcgTCGTATCAGGCCTACGCCAGACAGCGGAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATaaa                                                                   <  1:774311/237‑1 (MQ=255)
                                                                                                    tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCCCCTTTTGCTTGTTCTTCTCGATAAAA‑Cctt                                                             >  2:535259/1‑239 (MQ=255)
                                                                                                             gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTTTCGATAAAA‑CCTTCTTAACCGt                                                    >  1:91786/1‑239 (MQ=255)
                                                                                                                tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATCAAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAA‑CCTTCTTCACCGTTTa                                                 >  1:104113/1‑239 (MQ=255)
                                                                                                                 aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAA‑CCTTCTTCACCGTTTaa                                                <  1:570473/239‑1 (MQ=255)
                                                                                                                  ggCCAGACAGGGCAATTGCTTGATTTGGCTTTATTTTGTAGGGCGGATAAGGCGTTTCCG‑‑C‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAA‑CCTTCTTCACCGTTTAAg                                               <  2:588447/238‑1 (MQ=37)
                                                                                                                  cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAACCCCTTACATACAGCTCAGGTTCAACCAGTTCAGCGCCTTTTGCTTGTTTTTCTC‑ATAAAACCCTCCTTCACGGCTTAAg                                               >  2:131960/1‑239 (MQ=21)
                                                                                                                            cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAA‑CCTTCTTCACCGTTTCAGAACATCAtt                                      >  1:190657/1‑238 (MQ=255)
                                                                                                                            cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAA‑CCTTCTTCACCGTTTAAGAACATCAtt                                      >  1:404995/1‑238 (MQ=255)
                                                                                                                            cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAA‑CCTTCTTCACCGTTTAAGAACATCAtt                                      >  2:687891/1‑238 (MQ=255)
                                                                                                                             gcAATTGCCTGATTTCGCGTTATTTGGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGATATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACGCCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAA‑CCTTCTTCACCGTTTAAGAACATCATTc                                     <  2:619755/238‑1 (MQ=37)
                                                                                                                                 tggCCTGATTTGGCTTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAA‑CCTTCTTCACCGTTTAAGAACATCATTCCCGcc                                <  1:805207/237‑1 (MQ=255)
                                                                                                                                 tAGCCTGATTTGGCGTGATTTTGGAAGTCGGATAAGGCGGTT‑‑ACGC‑‑CAC‑‑A‑TCCGACATCAATGCCTGATTCCACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAA‑CCTTCTTCACCGTTTAAGAACATCATTCCCGcc                                <  1:488210/239‑1 (MQ=255)
                                                                                                                                       gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAA‑CCTTCTTCCCCGTTTAAAAACATCATTCCCGCCCAgaa                           >  2:158525/1‑236 (MQ=255)
                                                                                                                                       gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAA‑CCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTAc                          <  1:202458/239‑1 (MQ=255)
                                                                                                                                        aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGATTCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAC‑ACTTCTTAACCGTTTAACAACATCATTCCCGCCCAGTAcc                         >  2:657825/1‑238 (MQ=37)
                                                                                                                                              cGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAA‑CCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAAt                   >  1:704454/1‑239 (MQ=255)
                                                                                                                                                  aTTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGTTAAAA‑CCTTCTTCACCGTTTATGAACATCATTCCCGCCCAGAACTCGGAATAAcc               >  1:364572/1‑239 (MQ=255)
                                                                                                                                                      tGTAGGTCGGATAAGGCGTTTTCC‑‑C‑‑CCC‑‑A‑TCCGATATTAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTCGTTCTTCTCGATAAAA‑CCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGtt           <  2:336013/239‑1 (MQ=39)
                                                                                                                                                        tAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAA‑CCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCg         <  1:302613/239‑1 (MQ=255)
                                                                                                                                                            tCGGATAAGGCGTTT‑‑ACGC‑‑CGT‑‑C‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAA‑CCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGa      <  1:426025/238‑1 (MQ=255)
                                                                                                                                                             cGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAA‑CCTTCTTCa                                                        >  2:370846/1‑187 (MQ=255)
                                                                                                                                                             cGGATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAA‑CCTTCTTCa                                                        <  1:370846/187‑1 (MQ=255)
                                                                                                                                                               gATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCt                                                                                                                >  2:546632/1‑130 (MQ=255)
                                                                                                                                                               gATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCt                                                                                                                <  1:546632/130‑1 (MQ=255)
                                                                                                                                                               gATAAGGCGTTT‑‑ACGC‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGTTAAAA‑CCTTCTTCACCGTTTAAGAACATCATTCCCGCCCATTACTCGGAATAACAGGTTCGGGAAAtt  >  2:424213/1‑239 (MQ=255)
                                                                                                                                                                              ||                                                                                                                                                                                                                                      
ACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑‑CGC‑‑A‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAA‑CCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATT  >  NC_000913/4296209‑4296604

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