Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,296,381 +GC intergenic (+587/+55) 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.C93.0% 161.2 / 3.4 59intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base C (29/24);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
*NC_0009134,296,3802.G91.5% 160.1 / 6.0 59intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base G (29/25);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTA‑‑C‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAATGCCCACGCCGTTGGCGT  >  NC_000913/4296148‑4296671
                                                                                                                                                                                                                                          ||                                                                                                                                                                                                                                                                                                        
gCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCTCTTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTg                                                                                                                                                                                                                                    <  1:431763/300‑1 (MQ=9)
 cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGTTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCAACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                                                                                                                                                                                                                   <  1:103126/300‑1 (MQ=255)
 cGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGCTTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTCGGt                                                                                                                                                                                                                                  >  2:57447/1‑301 (MQ=9)
   ttgtagccgcgacgaaCATAAACTCCTGATGCGACGGTTGCCGCGTCTTACAAGGCCTACGCCAGACAGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                                                                                                                                                                                                                                <  2:99659/286‑1 (MQ=11)
     aacgcagCATCACACATCAAGGCCCGAAGCGACGCTGGCGTCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGTTTTGGCGTGATTTTGTAGGTCGAATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATTCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                                                                                                                                                                                                                              <  1:347658/295‑1 (MQ=11)
      ccgccgccACCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                                                                                                                                                                                                             <  1:57448/299‑1 (MQ=255)
                 gACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGCCGCCTAACGCGTTTTATCAGGCCTACGCCAGACCTCGCAATAGCCCGATTTAGCGTGATTTTGTCGGTCGGATAAGCCG‑CTTACGC‑‑C‑GC‑ATCCGACATCAATTCCTGATGCGATGCTTG‑CCGCGTCTCTACAGGCCTTTTTTAACCGTGTGTTAATTTCTAAAAct                                                                                                                                                                                                                    >  1:243504/1‑299 (MQ=11)
                    tctacaccccgaagagttgctagtcgagtgttATCAGGCATAGACAAGAAGGAGCAATAGCCTGATTTGGCGTGATTTTGTATGTCGGATATGGAGTTTATGCCGCATACGACATCAACGCCTGATGCGGCGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                                                                                                                                                                                                               <  2:373351/271‑1 (MQ=14)
                           ccTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACCGCGCAATACCCTGATTTCGCTTGATTTTGGAGGTCGGAGAAGGCG‑TTTA‑‑CG‑C‑CCAATCTGCCATCAATGCCTGCTGCGAC‑CTTGACCGCTTTCTATCCGGCCTATGTTACACGTTGGTTTATTTTTCatatctttcttcatctt                                                                                                                                                                                                        >  2:267357/1‑285 (MQ=1)
                              gaggCGACGCTTGCCGGGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCATATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                                                                                                                                                                                      <  1:191091/297‑1 (MQ=255)
                                     cGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATACGGCGTTTATGCCGCATCCGAAATCAACGCCTGATGCGACTATTACGGCGGCTTATCAAGACTATGCCCGACAGCCCAATAGCCTGATTCAGGGTGACTTTGTTGGACGGATCACACG‑TTTA‑‑CGAC‑GC‑CTCCTACATGAATACCCGATGCCACCCTTT‑CCGCATCATAACATGCATATCCTCCACCTTGGATAATTGTTCaccctcatacccgcaggtttcacaca                                                                                                                                                                                               >  2:490657/1‑275 (MQ=2)
                                     cGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTC‑ACGCGTCTTCTCAGGCCTATATTAACCGTTGGTTAATTTTTAAAACTTTTAACCCCCAGTTTCAAAcc                                                                                                                                                                                               >  2:490661/1‑300 (MQ=11)
                                          tcagacgtgtgctcttccgatctcAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCACACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCa                                                                                                                                                                                         <  1:394301/278‑1 (MQ=17)
                                            cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGTCAGACAGCGCAATAGCCTGCTTTATCGTGATTTTGTAGGTCGGATAAGTCG‑CTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACCCTGT‑CCGCGTCTTATCAGGCCTATTTTAACCGTTGGTTAATTNTTCAAACCCTTCCCACCCGGtcacaacacccacacca                                                                                                                                                                                       >  2:70553/1‑285 (MQ=14)
                                                tcttagacggCCTTCGCCACACGGCGCAATAGCCTGATTTGGCCGTATTTTGTTGGTCGGAATAGGCGTTTATGCCGCAGCCCAACTCAACGCCTTAAGCGACGCTTTACGCGTCTTATCAGGCCTAGGCCAGACAGAGGAATAGCTCGATTTAGCGTGAGTTTGTAGGTCGGTTAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGc                                                                                                                                                                                   <  1:488256/293‑1 (MQ=17)
                                                        ggCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACACTTCCATACAGCTCAGGtt                                                                                                                                                                            >  1:489064/1‑300 (MQ=255)
                                                           cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCaat                                                                                                                                                                        <  1:98527/301‑2 (MQ=255)
                                                                 cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCACACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACCCCCGGTATCAAACCCTTcaaagatcggaagaactgctgttagg                                                                                                                                                                  >  2:394292/1‑276 (MQ=14)
                                                                     aaggcgaaTAGCCTGATTAGGCGTGGATTTGAAGGGCGGAAAAGGGTAGTAGTCCGAATCAGACAGGAACGCCGGATGCGACGCTTAGCGCTTCTGATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgcg                                                                                                                                                              <  2:21113/295‑1 (MQ=18)
                                                                         gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCtt                                                                                                                                                          >  1:181074/1‑301 (MQ=255)
                                                                          ctttccctccacgccggtcttccgatctggTCGGAGAAGGCTTTCATGCCGCTTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCttt                                                                                                                                                         <  2:145560/273‑1 (MQ=18)
                                                                                acTGATTTAGCGTGATTTTGTAGGTCGTATAAGGCGTTTATGCCGCATCCGACATCAACGCCTTATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCt                                                                                                                                                     <  1:267361/298‑1 (MQ=255)
                                                                                  tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACATCTCCGGTTCAACCAGTGCAGCGCCTTTTGCTTGtt                                                                                                                                                 >  2:24533/1‑301 (MQ=255)
                                                                                          cgggATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCg                                                                                                                                          <  1:267953/297‑1 (MQ=255)
                                                                                                      ggTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTAGCtcggaagagcacacgtctgaactc                                                                                                                              >  1:145562/1‑277 (MQ=18)
                                                                                                                 gCGTTTATGCAGCAGCCGACATAAACGCCGGATGCGACGCTTAACGCGTCTTATCAAGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTaa                                                                                                                   <  2:181072/300‑1 (MQ=255)
                                                                                                                  cGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGa                                                                                                                 <  1:24533/301‑1 (MQ=255)
                                                                                                                   ttttATGCCGTATCGGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCATTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGaa                                                                                                                <  1:2813/300‑1 (MQ=255)
                                                                                                                         gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTCCAGCGCCTTTTGCTTTTTCTTTCCGATAAACCCTTCTTCACCGTTTAAGAActtcat                                                                                                          >  2:368748/1‑301 (MQ=255)
                                                                                                                                 cgatctCAACGCCTGATGAGGCGCTTAAAGCGTCTTATCAGGCCTACGGCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑GTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCGTATCTTAGCCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTGGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGccc                                                                                                  <  1:173303/296‑1 (MQ=25)
                                                                                                                                 cgaaagcAACGCATGATGCGACGGTGAACGCGGCTTGTGAGGCATACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTGTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GA‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGccc                                                                                                  <  2:159455/295‑1 (MQ=25)
                                                                                                                                      tCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCATGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCTAACCCTTCCATACAGCTCAGGTTCAACCCATGCAGCGCCTTTTTCTTGTTCTTCTCGATAAAACCTTCTTCCCCTGTTAAGAACATTCTTCCCGCCCAGTa                                                                                              >  2:299365/1‑300 (MQ=21)
                                                                                                                                       cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTACTTTTTCAAACTCTTCACACACGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTGAAGAACATCATTCCCACCCCGTa                                                                                              >  1:322461/1‑299 (MQ=255)
                                                                                                                                       cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTCAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCTTTCCATACAGCTCAGGTTCAACCCGTCCAGCGCCTTTTGCTTGTTCTTCTCGATAAACCCTTTTTACCCGTTTAATAACATCATTCCTGCCCAgatgg                                                                                            >  2:173301/1‑297 (MQ=21)
                                                                                                                                        aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCTGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACACTTCCATACAGCTCAGGTTCAACAAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCTCGCCCAGTACt                                                                                            >  1:329701/1‑300 (MQ=255)
                                                                                                                                             cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCTTTTTGCTTGTTCTTCTTGATAAAACCTTCTTCACCGTCTAAGAACATCATTCCCGCCCATTACCCGGaaa                                                                                      >  1:97703/1‑300 (MQ=255)
                                                                                                                                                                 tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCa                                                                  >  1:48063/1‑301 (MQ=255)
                                                                                                                                                                   ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCaga                                                                >  1:61724/1‑301 (MQ=255)
                                                                                                                                                                        aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTg                                                            >  1:48276/1‑300 (MQ=255)
                                                                                                                                                                         cgcccgcGCCACACACTGCAATAGCTTGATTTAGCTTGATTTTGTAGGGTGGATATGGCGATTCAG‑T‑‑C‑GC‑ATCCGACAACAACGCATGATGCGACGCTGG‑CCGCGTCTTATCAGGCCTATCATAACCGTTAGTGAATTTTTCAAACTCGCCACACCCGGTATCAAACCCCTCCATACAGCTCAGGTTCAACCAGAGCAGCGCCTTTCGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACAGCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCg                                                          <  1:490668/295‑1 (MQ=255)
                                                                                                                                                                          gCCTACGCCAGACCGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCg                                                          <  1:490672/300‑1 (MQ=255)
                                                                                                                                                                              acaccaagcAGCGCAGTCGTCTGATTTGGCGGGGTTTTTGCGGTGGGGTACGGCG‑TTTA‑‑CG‑CGGC‑ATCCGACAGCATTGCCTGGTGCGACGCTTG‑CCGCGTCTTATCAGGCCGATCTTACCCGTTGGGTAATTTTTCAACCTCTTCACGCCCGGTATCAAGCCCTTCCATACAGGTCAGGTTCAACCGGTGCAGCGGCTTTTGCTTGTTCTTCTCGGTAAGACCTTCTTCGCCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGa                                                     <  2:243501/293‑1 (MQ=255)
                                                                                                                                                                                gcaaGACAGCGAAATCGGCTGACTTAGGGTAATTTTGTAGGGCAGATAAGGCG‑CTTAAGC‑‑C‑GC‑GTCGAACATCAATGCCTGATGCGGCGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTAGTTAATTGTTCTAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCGGCGGCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACGTCATTCCCGCCCAATACTCGGAATAACCGGTTCGGGAAATTGCAGGGCTGCGTTTGaa                                                    <  2:97702/297‑1 (MQ=255)
                                                                                                                                                                                gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGaaa                                                   >  1:411644/1‑301 (MQ=255)
                                                                                                                                                                                          gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCTGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTg                                         >  1:51407/1‑301 (MQ=255)
                                                                                                                                                                                                     aTTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCCTACAGCTCATGTTCAACCAGTGCAGCGCATTTTTCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAATCACATCATTCCATCTCAGTAATGGGCATACCCGGTTCGGGAAATTGCAAAACTGGGTTTCAAATTCCAGGTTCCCttttagtat                              >  2:448904/1‑294 (MQ=25)
                                                                                                                                                                                                           ctgtaTTTTGTAGAGCGGATAATGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGc                         <  2:94136/296‑1 (MQ=255)
                                                                                                                                                                                                             tGATTTTGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCt                       <  1:70554/300‑1 (MQ=255)
                                                                                                                                                                                                                   tggaGGGCGTATAAGTCG‑TTTACGC‑‑C‑TC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAATGcc                <  2:51406/298‑1 (MQ=255)
                                                                                                                                                                                                                   tGTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCCTCGTCGGCCTTAATGc                 >  1:419787/1‑300 (MQ=255)
                                                                                                                                                                                                                    gTAGGTCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGTTTGCCTTGTCATCGTCGGCCTTACTGcc                >  2:283509/1‑300 (MQ=255)
                                                                                                                                                                                                                      atctCGGATAAGGCT‑TGTA‑‑CGAC‑TC‑ATCCAACATCAATGCCTGATTCGACGCTTA‑CCGCGTATTATCAGGCCTATCTTAACCTTTGGTTAATTGTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAATGcc                <  1:132169/295‑1 (MQ=255)
                                                                                                                                                                                                                         tCGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CTGCGTCTTATCAGGTCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGCATCAAACCCTCCCATACAGCTCAGGTTCAACCAGTGTAGTGCATTTTGCTTGTTCTCTTCGATAAAATCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACCTGGAATAAACGGTTCGGGAAATTGCCGAGCTGCGTCTGAAATACGAGGTCGCCCAGTAGTCGTCGGCTTTacttctctcgc           >  2:91973/1‑290 (MQ=255)
                                                                                                                                                                                                                          cGGATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTCAAGAACATCATTCCCGCCCCGTACTCGGAATAACCGGTTCGGGAAATTGAAGAGCTGCGCTTGAATCACCTGGTTGCCTATACATCCTCGGCCTTAAttacataacgg         >  2:132167/1‑291 (MQ=38)
                                                                                                                                                                                                                            gATAAGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAATGCCCACGCCGtt       >  1:482774/1‑301 (MQ=255)
                                                                                                                                                                                                                              tctggCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAATGCCCACGCCGTTgg     <  1:4946/298‑1 (MQ=255)
                                                                                                                                                                                                                                aGGCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCCCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTTTCATCGTGGGCCTTAATGCCCAAGCCGTTgggg   >  2:1526/1‑299 (MQ=255)
                                                                                                                                                                                                                                 ggCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATGTTCGGCCTTAATGCCAACGCCGTTggaga  >  2:4946/1‑298 (MQ=255)
                                                                                                                                                                                                                                 ggCG‑TTTACGC‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAATGCCCACGCCGTTGGCGt  >  1:140791/1‑301 (MQ=255)
                                                                                                                                                                                                                                          ||                                                                                                                                                                                                                                                                                                        
GCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCG‑TTTA‑‑C‑‑C‑GC‑ATCCGACATCAATGCCTGATGCGACGCTTG‑CCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAATGCCCACGCCGTTGGCGT  >  NC_000913/4296148‑4296671

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