Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,115,605 A→C 6.2% intergenic (‑115/+57) mdtG ← / ← lpxL putative drug efflux system protein/lauryl‑acyl carrier protein (ACP)‑dependent acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,115,6050AC 6.2% 227.6 / 6.1 76intergenic (‑115/+57)mdtG/lpxLputative drug efflux system protein/lauryl‑acyl carrier protein (ACP)‑dependent acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (34/37);  new base C (2/3);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

TATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACTGCTCTGGTGCCATCATGATGCATTTTT  >  NC_000913/1115461‑1115750
                                                                                                                                                |                                                                                                                                                 
tataGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAgg                                                                                                                                                 <  8:181819/147‑1 (MQ=255)
    ggggTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGa                                                                                                                                             <  4:249061/147‑1 (MQ=255)
      ggTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGAttt                                                                                                                                          <  6:196927/148‑1 (MQ=255)
        tgtCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATCAGCAGAAAAGATACGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGc                                                                                                                                        <  2:347190/148‑1 (MQ=255)
         gtCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCt                                                                                                                                       >  2:271097/1‑148 (MQ=255)
            aTTTTCACAGGGTGACATAGCAATCCGCTGTTGGGGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCtgt                                                                                                                                     <  7:302687/147‑1 (MQ=255)
             ttttCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTgg                                                                                                                                   >  5:195210/1‑148 (MQ=255)
               ttCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGTGCCCGACTTCCAGGGGGCGCAATCCAGATAGCTTTTATAGCTAAATCAGGTGGATTTGCTGTGGt                                                                                                                                  <  2:372650/147‑1 (MQ=255)
                tCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTa                                                                                                                                 <  5:285401/147‑1 (MQ=255)
                tCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTa                                                                                                                                 <  2:357622/147‑1 (MQ=255)
                      ggTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGTGGGATTTGCTGTGGTAATGCCgg                                                                                                                          <  1:297882/148‑1 (MQ=255)
                      ggTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCgg                                                                                                                          <  1:176307/148‑1 (MQ=255)
                      ggTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCgg                                                                                                                          <  1:94323/148‑1 (MQ=255)
                      ggTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCgg                                                                                                                          <  1:303004/148‑1 (MQ=255)
                             tAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCAt                                                                                                                   >  7:302401/1‑148 (MQ=255)
                              aGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCAtt                                                                                                                  >  3:45400/1‑148 (MQ=255)
                              aGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCAtt                                                                                                                  >  5:240606/1‑148 (MQ=255)
                              aGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCAtt                                                                                                                  >  8:352740/1‑148 (MQ=255)
                              aGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCGAAATCAGGTGGATTTGCTGTgg                                                                                                                                   >  4:14898/1‑131 (MQ=255)
                               gCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTc                                                                                                                 >  8:222340/1‑148 (MQ=255)
                               gCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTc                                                                                                                 >  1:301563/1‑148 (MQ=255)
                                       cTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAg                                                                                                                                                  <  3:230468/107‑1 (MQ=255)
                                       cTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAg                                                                                                                                                  >  4:230468/1‑107 (MQ=255)
                                          ttGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCg                                                                                                      <  4:45400/148‑1 (MQ=255)
                                               gcgcCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCCGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATg                                                                                                 <  8:220931/148‑1 (MQ=255)
                                                 gcCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATgg                                                                                                >  7:349071/1‑147 (MQ=255)
                                                         cgGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTaa                                                                                       <  2:301563/148‑1 (MQ=255)
                                                           gTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTaa                                                                                       >  8:299134/1‑146 (MQ=255)
                                                           gTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTaa                                                                                       >  3:182934/1‑146 (MQ=255)
                                                            tGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAg                                                                                    <  6:240606/148‑1 (MQ=255)
                                                                 aTAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGa                                                                               <  2:86459/148‑1 (MQ=255)
                                                                  tAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTg                                                                                >  7:357238/1‑146 (MQ=255)
                                                                   aagaagAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAg                                                                             >  8:45746/1‑148 (MQ=255)
                                                                    agaagaAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAg                                                                             >  3:82796/1‑147 (MQ=255)
                                                                      aagaaAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGaa                                                                          >  2:155913/1‑148 (MQ=255)
                                                                      aagaaAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGaa                                                                          >  1:223707/1‑148 (MQ=255)
                                                                          aaaGATAAGCACCCTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATCCCGGATGCCATTATGAAGCATACGGGATGGGAGATTTAATAGAGTGAAGGAAAGcc                                                                      >  6:99058/1‑148 (MQ=255)
                                                                            aGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCtt                                                                    >  6:230748/1‑148 (MQ=255)
                                                                                 aGCACACTAATTATGCGCCCGACTTCCATGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCgg                                                                <  1:334945/147‑1 (MQ=255)
                                                                                 aGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCggg                                                               >  4:339631/1‑148 (MQ=255)
                                                                                  gCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGTATGCCATTCTGAAGCATCCGGAATGGGAGATTTAATAGCGTGAAGTAACGCCTTCCGGGc                                                              >  4:207508/1‑148 (MQ=255)
                                                                                    acacTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCgt                                                            <  7:352740/148‑1 (MQ=255)
                                                                                      acTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCgtgt                                                          >  1:300708/1‑148 (MQ=255)
                                                                                       cTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGtt                                                         <  5:140583/148‑1 (MQ=255)
                                                                                       cTAATTATGCGCCCGACTTCCAGGGGGCACAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGtt                                                         >  3:245225/1‑148 (MQ=255)
                                                                                         aaTTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGttt                                                        >  1:21545/1‑147 (MQ=255)
                                                                                          aTTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGtttt                                                       >  5:84288/1‑147 (MQ=255)
                                                                                                cgcCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCg                                                 <  5:323012/147‑1 (MQ=255)
                                                                                                   ccGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACgg                                             <  2:300708/148‑1 (MQ=255)
                                                                                                     gACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGgtg                                           <  8:349071/148‑1 (MQ=255)
                                                                                                       cTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTa                                         <  1:60723/148‑1 (MQ=255)
                                                                                                       cTTCCAGGGGGCGCAATCCAGAGAGCTTTAATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTa                                         <  6:362987/148‑1 (MQ=255)
                                                                                                         tCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAAc                                       >  7:173016/1‑148 (MQ=255)
                                                                                                          ccAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAAc                                       <  1:326189/147‑1 (MQ=255)
                                                                                                           cAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAAcc                                      <  4:182934/147‑1 (MQ=255)
                                                                                                           cAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACca                                     <  7:222340/148‑1 (MQ=255)
                                                                                                            aGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACcac                                    <  4:82796/148‑1 (MQ=255)
                                                                                                              ggggCGCAATCCAGAGAGCTTTTATCGCTAAATCCGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGAc                                               <  3:163984/135‑1 (MQ=255)
                                                                                                              ggggCGCAATCCAGAGAGCTTTTATCGCTAAATCCGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGAc                                               <  3:158297/135‑1 (MQ=255)
                                                                                                              ggggCGCAATCCAGAGAGCTTTTATCGCTAAATCCGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGAc                                               >  4:163984/1‑135 (MQ=255)
                                                                                                              ggggCGCAATCCAGAGAGCTTTTATCGCTAAATCCGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGAc                                               >  4:158297/1‑135 (MQ=255)
                                                                                                              ggggCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACat                                  <  7:45746/148‑1 (MQ=255)
                                                                                                               gggCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACata                                 >  2:221151/1‑148 (MQ=255)
                                                                                                                  cgcAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATAc                              >  8:68066/1‑148 (MQ=255)
                                                                                                                      aTCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCGTTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACTGct                          >  8:334150/1‑148 (MQ=255)
                                                                                                                       tCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACTGctc                         >  7:280765/1‑148 (MQ=255)
                                                                                                                        ccAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACTGctct                        >  1:344990/1‑148 (MQ=255)
                                                                                                                         cAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACTGCTCTg                       <  2:337528/148‑1 (MQ=255)
                                                                                                                          agagagCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACTGCTCTgg                      <  6:195210/148‑1 (MQ=255)
                                                                                                                                   ttATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACTGCTCTGGTGCcatc               <  6:84288/146‑1 (MQ=255)
                                                                                                                                   ttATCGCTAAATCAGGGGGATTTGCTGGGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACTGCTCTGGTGCcatc               <  5:224922/146‑1 (MQ=255)
                                                                                                                                    tATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACTGCTCTGGTGCcatcat             <  1:271097/147‑1 (MQ=255)
                                                                                                                                      tCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACTGCTCTGGTGCCATCatgat          >  2:377344/1‑148 (MQ=255)
                                                                                                                                           aaaTCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACTGCTCTGGTGCCATCATGATGCAtt     >  4:75227/1‑148 (MQ=255)
                                                                                                                                            aaTCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGTCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCAGATATACTGCTCTGGTGCCATCATGATGCAtt     <  3:207508/147‑1 (MQ=255)
                                                                                                                                              tCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACTGCTCTGGTGCCATCATGATGCAttttt  <  8:302401/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACTGCTCTGGTGCCATCATGATGCATTTTT  >  NC_000913/1115461‑1115750

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: