Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,296,428 G→A intergenic (+634/+8) gltP → / ← yjcO glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,4280GA89.7% 85.9 / 3.5 39intergenic (+634/+8)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base A (17/18);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.17e-01

CGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCT  >  NC_000913/4296192‑4296650
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                              
cgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAAc                                                                                                                                                                                                                    >  1:234084/1‑251 (MQ=2)
              ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAAtttt                                                                                                                                                                                                      >  1:47512/1‑251 (MQ=255)
                 aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCTTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAAttttt                                                                                                                                                                                                     <  1:262884/249‑1 (MQ=255)
                                atatccTGTTTTAGCGTTATTTTGTAGTGCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCg                                                                                                                                                                                    <  2:258482/247‑1 (MQ=14)
                                agtcctagATTGAGCGTGAGTATGTGGGGCGGTTAAGGCGTTGATGCCGACTCCCTCATCCACTCCTTGTGCGACGTCTAACGCGTCTGATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATGTTGTAGGTCGGATAAGGCGTTTACTC‑‑C‑GCAACCGACATCAATGCCTGATGCGACTCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCacac                                                                                                                                                                                       <  2:236842/241‑1 (MQ=9)
                                aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTCATTTTTCAAACTCTTCACACCCg                                                                                                                                                                                    >  1:284782/1‑251 (MQ=255)
                                      tGATTTAGCGTGATTTTGTAGGTCGGATAAGGAGTTTATGCCGCATCCGACATCAACGACTGATGCGCCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTTATTTAGCGTAATTTTGTCGATCGGATAAGGCGGTAA‑‑CG‑CAGCATCTGACATCACTGCGTGCTGCGACGCTGGCAGCGTCTTATCAAGCCTCTCTTCTCCGTTGGTAAATTATTCAAACTATTCACACCCAGTATc                                                                                                                                                                               >  2:361770/1‑250 (MQ=2)
                                              cGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTc                                                                                                                                                                      <  1:1085/251‑1 (MQ=255)
                                               gTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                                                                                                                                                                     >  1:266299/1‑251 (MQ=255)
                                                  aTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                                                                                                                                                                  <  1:322833/251‑1 (MQ=255)
                                                   ttttGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTACTTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                                                                                                                 >  1:183509/1‑251 (MQ=255)
                                                             cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGtt                                                                                                                                                       <  1:88962/251‑1 (MQ=255)
                                                                 tAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc                                                                                                                                                   >  1:303231/1‑251 (MQ=255)
                                                                    gacGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATAAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGATATCAAACCCTTCCCTACAGCTCAGGTTCAACCa                                                                                                                                                 >  1:358332/3‑250 (MQ=18)
                                                                         ttATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGACTCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGc                                                                                                                                             >  2:382051/1‑249 (MQ=18)
                                                                                     cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTCTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCtt                                                                                                                               <  1:364687/251‑1 (MQ=255)
                                                                                            aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCACTGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGATTCAACCAGTGCAGCGCCTTTTGCCTGttcttc                                                                                                                        >  2:134783/1‑251 (MQ=21)
                                                                                               gCCTGATGCGACGATTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGACTGATTTAGCGTGATTTTGTAGGACGGATAAGGCGTTTA‑‑CG‑CAGCATCCGACATCAATGCCTGATGCGACGCGTGCCGAGTCTTATCAAGCATCTCTTAACCGTTGGTTAATTTTTCAACCAATTAACACCCGGTATCCTACCCTTCCCAAAAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTGGttcttcacg                                                                                                                     >  2:306455/1‑248 (MQ=17)
                                                                                                       cgacgaTTAACGAGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAAcc                                                                                                             <  1:275442/251‑1 (MQ=255)
                                                                                                        gACGCTTAACGCGTCTTATCCGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACct                                                                                                            >  1:330229/1‑251 (MQ=255)
                                                                                                         aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACctt                                                                                                           <  1:382051/251‑1 (MQ=255)
                                                                                                                   cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGtt                                                                                                 >  2:13551/1‑251 (MQ=255)
                                                                                                                     tCTTATCAGGCCTACGCCAGACAGCGCAATAGACTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCATATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTAAATAAAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTATCGATAAAACCTTCTTAACCGttt                                                                                                >  1:328334/1‑250 (MQ=21)
                                                                                                                        tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTaa                                                                                              <  1:13551/249‑1 (MQ=255)
                                                                                                                                  gcGCCAGACAGCGCTATAGCCTGATTTAGCGTGATTTTGTAGGTCGGTTAAGTCGTTTACGC‑‑C‑GCATCCTACATCAATGCCTTATTCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAATCCTTCTTCACCGTTTAAGAACATCATTc                                                                                   <  2:303231/249‑1 (MQ=37)
                                                                                                                                      cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGcc                                                                              <  1:306455/251‑1 (MQ=255)
                                                                                                                                           agggcAATAGCCTGATTTAGGTTGATTTTGTAGGTCGGAAAATGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCGTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTACCGCCCAGTa                                                                         <  2:234084/248‑1 (MQ=37)
                                                                                                                                                   aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAAt                                                                 <  1:272894/251‑1 (MQ=255)
                                                                                                                                                     ccTGATTTAGTGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCGCACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATa                                                                >  1:114736/1‑250 (MQ=255)
                                                                                                                                                      cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCATTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGACCAGTACTCGGAATAAc                                                              >  2:324308/1‑251 (MQ=255)
                                                                                                                                                             agtgTGATTTTGTAGGTCGGATAAGGCGTGTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACTCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCGCACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTc                                                        <  2:114736/247‑1 (MQ=255)
                                                                                                                                                              gcttGATTTTGTAGGTCGGATAAGGCGCTTACGC‑‑C‑GCATCCGACGTCAATGCCGGATGCGACTCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTTGTTAACTTTTCCAACTCTTCACACCCTGTATCCAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTGTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGATCATCATTCCCGCCCAGTACGCGGAATAACCGGTTCgg                                                      <  2:293963/248‑1 (MQ=255)
                                                                                                                                                                  gattgtggaGGTCGGATAAGGCTGTGA‑‑CCAC‑GCATCCGACAGCAATGCCTGATGCTACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCCCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGaa                                                   <  2:38962/242‑1 (MQ=255)
                                                                                                                                                                  gATTTTGTAGGTCTGATAAGTCGTTTACGC‑‑C‑GCAGCCGACATCAATGCTTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCATCGCCTTTTGCGTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCTTTCCCGCCCAGTACTCGGAATAACCGGTTCGGGaa                                                   <  2:124590/250‑1 (MQ=255)
                                                                                                                                                                                gagaaGGCGTTTACGC‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCTTATCTTAACCGTTGGTTAATTTTTCAACCTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGTAAATTGCAGAGCTGCGt                                    <  1:84032/248‑1 (MQ=255)
                                                                                                                                                                                                         cGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGt              >  1:376258/1‑251 (MQ=255)
                                                                                                                                                                                                              tCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCAtcg         >  1:310390/1‑251 (MQ=255)
                                                                                                                                                                                                                    ccTGATGCGACGCTTGCCGCGTCTTATCAAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGcc   <  1:324308/251‑1 (MQ=255)
                                                                                                                                                                                                                       gATGCGACGCTTGCCGCGTCTTATCTAGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCt  <  2:266299/249‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                              
CGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCT  >  NC_000913/4296192‑4296650

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