Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,997,022 C→T 6.5% intergenic (‑190/+40) sdiA ← / ← yecC quorum‑sensing transcriptional activator/putative ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,997,0220CT 6.5% 188.8 / 2.4 61intergenic (‑190/+40)sdiA/yecCquorum‑sensing transcriptional activator/putative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (28/29);  new base T (2/2);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAAACTATCACATAAATAAGATTTATATATAATTTATATTATTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTGGTCCATAAAG  >  NC_000913/1996880‑1997170
                                                                                                                                              |                                                                                                                                                    
tCAAACTATCACATAAATAAGATTTATATATAATTTATATTATTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGcccc                                                                                                                                                >  1:12599/1‑149 (MQ=255)
 cAAACAATAACATAAATAAGATTTATATATAATTTATATTATTAAGGCAATTAAATACTTTTTAAATACATCGCATTATCTTATGTTAGTAATGATTTATGCTGATATGTTAAGCAGTTTTATCAGGCTTGAAATTGAGTGCAGaccc                                                                                                                                                <  5:202698/148‑1 (MQ=255)
  aaaCTATCACATAAATAAGATTTATATATAATTTATATTATTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCTAGCCCCGa                                                                                                                                              >  7:172153/1‑149 (MQ=255)
  aaaCTATCACATAAATAAGATTTATATATAATTTATATTATTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCTAGCCCCGa                                                                                                                                              >  7:129047/1‑149 (MQ=255)
   aaCTATCACATAAATAAGATTTATATATAATTTATATTATTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGAc                                                                                                                                             <  6:164607/149‑1 (MQ=255)
      tATCACATAAATAAGATTTATATATAATTTATATTATTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAgg                                                                                                                                          >  6:169098/1‑149 (MQ=255)
        tCACATAAATAAGATTTATATATAATTTATATTATTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGt                                                                                                                                         <  3:188659/148‑1 (MQ=255)
            ataaataaGATTTATATATAATTTATATTATTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAAtc                                                                                                                                    <  6:11553/149‑1 (MQ=255)
            ataaataaGATTTATATATAATTTATATTATTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAAtc                                                                                                                                    >  2:72443/1‑149 (MQ=255)
             taaataaGATTTATATATAATTTATATTATTCAGGCAATGAATTACTTTCGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAgg                                                                                                                                          <  3:108549/142‑1 (MQ=255)
             taaataaGATTTATATATAATTTATATTATTCAGGCAATGAATTACTTTCGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAgg                                                                                                                                          >  4:108549/1‑142 (MQ=255)
               aataaGATTTATATATAATTTATATTATTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAAtcgtc                                                                                                                                 >  8:18673/1‑149 (MQ=255)
                     aTTTATATATAATTTATATTATTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCt                                                                                                                           >  4:53654/1‑149 (MQ=255)
                            tataATTTATATTATTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGAttttt                                                                                                                    <  7:224239/149‑1 (MQ=255)
                             ataATTTATATTATTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGAttttt                                                                                                                    >  6:218997/1‑148 (MQ=255)
                                aTTTATATTATTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCTAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCt                                                                                                                 <  8:172153/148‑1 (MQ=255)
                                aTTTATATTATTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCTAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCt                                                                                                                 <  8:129047/148‑1 (MQ=255)
                                aTTTATATTATTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCt                                                                                                                 >  8:236/1‑148 (MQ=255)
                                  ttATATTATTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCttat                                                                                                              >  8:107900/1‑149 (MQ=255)
                                     tattatTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTtattat                                                                                                           <  2:12599/149‑1 (MQ=255)
                                     tattatTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTGATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTtattat                                                                                                           >  7:228874/1‑149 (MQ=255)
                                                cAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGaaa                                                                                                <  5:169098/149‑1 (MQ=255)
                                                      aTTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTtctc                                                                                          >  7:145859/1‑149 (MQ=255)
                                                       ttACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTtctc                                                                                          >  4:160035/1‑148 (MQ=255)
                                                          cTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAg                                                                                       >  7:108996/1‑148 (MQ=255)
                                                          cTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAg                                                                                       <  5:218997/148‑1 (MQ=255)
                                                           ttttGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAgg                                                                                      >  1:117098/1‑148 (MQ=255)
                                                             ttGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAAc                                                                                   >  6:40139/1‑149 (MQ=255)
                                                              tGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACATGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAAc                                                                                   <  2:221463/148‑1 (MQ=255)
                                                               gCAAGACATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACt                                                                                  <  5:186042/148‑1 (MQ=255)
                                                                 aaGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGc                                                                               <  3:84975/149‑1 (MQ=255)
                                                                 aaGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGc                                                                               <  5:62449/149‑1 (MQ=255)
                                                                 aaGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGc                                                                               <  1:72443/149‑1 (MQ=255)
                                                                       tCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGc                                                                         >  2:22013/1‑149 (MQ=255)
                                                                       tCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGc                                                                         >  5:186702/1‑149 (MQ=255)
                                                                                 ttATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGgc                                                                    >  6:170942/1‑144 (MQ=255)
                                                                                 ttATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGgc                                                                    <  5:170942/144‑1 (MQ=255)
                                                                                   aTGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCg                                                              <  1:22013/148‑1 (MQ=255)
                                                                                   aTGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTAgg                                                             >  8:116817/1‑149 (MQ=255)
                                                                                   aTGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATAGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCgg                                                             >  1:34664/1‑149 (MQ=255)
                                                                                    tGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCggg                                                            >  8:66395/1‑149 (MQ=255)
                                                                                      ttattaATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGt                                                          >  8:134199/1‑149 (MQ=255)
                                                                                         ttaATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCg                                                        <  3:53654/148‑1 (MQ=255)
                                                                                               aGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATa                                                 <  7:107900/149‑1 (MQ=255)
                                                                                                 ttATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAaacaac                                               <  7:18673/149‑1 (MQ=255)
                                                                                                 ttATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAAc                                               >  3:57623/1‑149 (MQ=255)
                                                                                                 ttATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAAc                                               <  2:218404/149‑1 (MQ=255)
                                                                                                      cTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCttt                                          >  8:208993/1‑149 (MQ=255)
                                                                                                      cTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCtt                                           <  8:108996/148‑1 (MQ=255)
                                                                                                       tGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCtttt                                         >  4:107424/1‑149 (MQ=255)
                                                                                                        gATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTg                                        >  3:78318/1‑149 (MQ=255)
                                                                                                         aTTTTTTAAGCAGTTTTATCAGGCGTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCtttt                                         <  5:206184/147‑1 (MQ=255)
                                                                                                         aTTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTg                                        <  7:236/148‑1 (MQ=255)
                                                                                                                        ttATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGa                         <  8:228874/148‑1 (MQ=255)
                                                                                                                          ctCAGGCTTGAAATGGCGGCCAGCCCCGACAGGTGAAGCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTa                      <  7:71727/148‑1 (MQ=255)
                                                                                                                          aTCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACt                       <  7:134199/148‑1 (MQ=255)
                                                                                                                             aGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGAcga                      >  8:71727/1‑144 (MQ=255)
                                                                                                                             aGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATc                    <  7:208993/148‑1 (MQ=255)
                                                                                                                                  tGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGcccc              >  8:35263/1‑149 (MQ=255)
                                                                                                                                         gCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTGGTCCa       <  4:78318/149‑1 (MQ=255)
                                                                                                                                           gTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTGGTCCATa     <  8:145859/149‑1 (MQ=255)
                                                                                                                                             ccAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTGGTCCATaaa   <  6:186702/149‑1 (MQ=255)
                                                                                                                                              cAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTGGTCCATAAAg  <  7:116817/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TCAAACTATCACATAAATAAGATTTATATATAATTTATATTATTCAGGCAATGAATTACTTTTGCAAGCCATCGCATTCTCTTATGTTATTAATGAGTTATGCTGATTTGTTAAGCAGTTTTATCAGGCTTGAAATGGCGTCCAGCCCCGACAGGTGAATCGTCGGGGCTGATTTTTTCTTATTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTGGTCCATAAAG  >  NC_000913/1996880‑1997170

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