Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,600,288 C→T S128L (TCG→TTG)  gcvR → transcriptional repressor, regulatory protein accessory to GcvA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,600,2880CT97.8% 307.6 / ‑2.9 90S128L (TCG→TTG) gcvRtranscriptional repressor, regulatory protein accessory to GcvA
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (48/40);  total (49/41)
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.82e-01

GTCTGGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTCGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGAGC  >  NC_000913/2600191‑2600383
                                                                                                 |                                                                                               
gtctgGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACCTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTtgc                                                                                               >  2:194397/1‑100 (MQ=255)
  ctgGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTgcgc                                                                                             >  1:1283216/1‑100 (MQ=255)
  ctgGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTATTCGACGAGCATCATATGAACATTGCGAAGCTGGTGTTgcgc                                                                                             >  1:84634/1‑100 (MQ=255)
   tgGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgca                                                                                            >  2:2020677/1‑100 (MQ=255)
    gggTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcac                                                                                           >  2:1021246/1‑100 (MQ=255)
     ggTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcacg                                                                                          >  1:1229922/1‑100 (MQ=255)
      nttCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcacgc                                                                                         >  1:1237551/2‑100 (MQ=255)
      gTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcacgc                                                                                         >  2:1162723/1‑100 (MQ=255)
      gTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcacgc                                                                                         >  2:835197/1‑100 (MQ=255)
       ttCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGTTTCACAGAACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcacgca                                                                                        <  2:441902/100‑1 (MQ=255)
        tCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGAACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCaa                                                                                       >  2:194383/1‑100 (MQ=255)
         cAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAAc                                                                                      >  2:147493/1‑100 (MQ=255)
          agacagATGTGACAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAAcc                                                                                     <  2:148743/95‑1 (MQ=255)
          aGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAAcc                                                                                     <  2:1699275/100‑1 (MQ=255)
           ggTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACct                                                                                    >  1:1602424/1‑100 (MQ=255)
           ggTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACct                                                                                    >  1:1195622/1‑100 (MQ=255)
           ggTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACct                                                                                    >  1:1727741/1‑100 (MQ=255)
              cGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACctgct                                                                                 <  1:322031/100‑1 (MQ=255)
               gATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCTCAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACctgctg                                                                                <  2:1130506/100‑1 (MQ=255)
               gATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACctgctg                                                                                >  2:1602881/1‑100 (MQ=255)
                  gtgGCAGACTCCCCGCATTTAATTGAACGCTTCACGGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGaaa                                                                             <  1:1878548/100‑1 (MQ=255)
                  gtgGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGaaa                                                                             >  1:220810/1‑100 (MQ=255)
                   tgGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGaaaa                                                                            >  1:650832/1‑100 (MQ=255)
                    ggCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAAt                                                                           >  2:229184/1‑100 (MQ=255)
                    ggCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAAt                                                                           >  1:42189/1‑100 (MQ=255)
                    ggCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAAt                                                                           >  2:260715/1‑100 (MQ=255)
                     gCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATg                                                                          <  1:1773969/100‑1 (MQ=255)
                        gACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGagc                                                                       >  2:1972442/1‑98 (MQ=255)
                        gACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGagc                                                                       >  2:1972434/1‑98 (MQ=255)
                         aCTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACg                                                                      <  1:1170155/100‑1 (MQ=255)
                            ccccGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGc                                                                   >  1:1519137/1‑100 (MQ=255)
                            ccccGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGc                                                                   >  1:818333/1‑100 (MQ=255)
                              ccGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTg                                                                 >  1:1661617/1‑100 (MQ=255)
                              ccGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTg                                                                 >  2:640296/1‑100 (MQ=255)
                              ccGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTg                                                                 >  2:2026283/1‑100 (MQ=255)
                                  aTTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCa                                                             >  1:1028643/1‑100 (MQ=255)
                                  aTTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCa                                                             <  1:1602893/100‑1 (MQ=255)
                                  aTTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCa                                                             >  2:799802/1‑100 (MQ=255)
                                   tttAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAg                                                            <  2:1354974/100‑1 (MQ=255)
                                    ttAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGt                                                           <  1:617935/100‑1 (MQ=255)
                                     tAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGtt                                                          >  2:1330834/1‑100 (MQ=255)
                                     tAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGtt                                                          >  1:213228/1‑100 (MQ=255)
                                     tAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGtt                                                          <  1:1021251/100‑1 (MQ=255)
                                     tAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGCGCTTCGCAGtt                                                          <  2:1075956/100‑1 (MQ=255)
                                     tAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATGGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATCAACGGGCGGCGCa                                                             <  2:1259138/97‑1 (MQ=255)
                                      aaTTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTg                                                         >  1:1274120/1‑100 (MQ=255)
                                      aaTTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTg                                                         >  1:887519/1‑100 (MQ=255)
                                      aaTTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTg                                                         <  1:1344464/100‑1 (MQ=255)
                                      aaTTGAACGCTTCACAGCACTCTTCGACGCGCATCATATGAACATTGCGCAGCTGGTGTCGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTg                                                         >  2:1296689/1‑100 (MQ=255)
                                         tGAACGCTTCACAGCACTTTACGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCat                                                      >  2:1249915/1‑100 (MQ=255)
                                          gAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCata                                                     >  1:1026403/1‑100 (MQ=255)
                                           aaCGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCatat                                                    <  2:887516/100‑1 (MQ=255)
                                            aCGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATAtt                                                   <  1:1530509/100‑1 (MQ=255)
                                             cGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAATGAACGGGCTGCGCAGTTGCATAtt                                                   <  1:701393/98‑1 (MQ=255)
                                              gCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCa                                                 <  1:2020695/100‑1 (MQ=255)
                                              gCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCa                                                 <  1:1162729/100‑1 (MQ=255)
                                              gCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCa                                                 <  1:194398/100‑1 (MQ=255)
                                               cTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAg                                                <  1:1200712/100‑1 (MQ=255)
                                                ttCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGa                                               <  2:609464/100‑1 (MQ=255)
                                                ttCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGa                                               <  1:1788339/100‑1 (MQ=255)
                                                      gCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCg                                         >  2:728722/1‑100 (MQ=255)
                                                       cACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGc                                        >  1:1329776/1‑100 (MQ=255)
                                                         cTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGccc                                      >  1:1118754/1‑100 (MQ=255)
                                                           tttCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCcac                                    <  1:835200/100‑1 (MQ=255)
                                                             tCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAg                                  <  2:42189/100‑1 (MQ=255)
                                                               gACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGcc                                <  1:1172928/100‑1 (MQ=255)
                                                                 cgcgCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGcccc                              <  2:1229916/100‑1 (MQ=255)
                                                                  gcgcATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCg                             <  2:912970/100‑1 (MQ=255)
                                                                  gcgcATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCg                             <  1:1277679/100‑1 (MQ=255)
                                                                   cgcATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGc                            <  2:1727728/100‑1 (MQ=255)
                                                                   cgcATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGc                            <  1:55709/100‑1 (MQ=255)
                                                                     catcatATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCAt                          <  2:818330/100‑1 (MQ=255)
                                                                      atcatATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATc                         <  1:1970915/100‑1 (MQ=255)
                                                                       tcatATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCt                        <  1:728724/100‑1 (MQ=255)
                                                                       ncatATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCt                        >  1:364924/2‑100 (MQ=255)
                                                                          tatGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCg                     >  1:699412/1‑100 (MQ=255)
                                                                          tatGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCg                     <  2:838707/100‑1 (MQ=255)
                                                                          tatGAACATTGCGGAGCTGGTGTCGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCg                     <  1:749593/100‑1 (MQ=255)
                                                                           atGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCgg                    >  1:1325232/1‑100 (MQ=255)
                                                                             gAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGAc                  >  1:1360625/1‑100 (MQ=255)
                                                                             gAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGAc                  >  2:444911/1‑100 (MQ=255)
                                                                               aCATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACgc                >  2:300485/1‑100 (MQ=255)
                                                                                  ttGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACgcagn             <  2:1519128/100‑2 (MQ=255)
                                                                                  ttGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACgcagc             >  1:1579409/1‑100 (MQ=255)
                                                                                    gCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAGCGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCaa           <  1:1972452/100‑1 (MQ=255)
                                                                                    gCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAGCGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCaa           <  1:1972460/100‑1 (MQ=255)
                                                                                    gCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCaa           >  2:1324351/1‑100 (MQ=255)
                                                                                     cGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCaaa          >  1:1620694/1‑100 (MQ=255)
                                                                                           tGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGa    <  2:213227/100‑1 (MQ=255)
                                                                                           tGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGa    >  1:843848/1‑100 (MQ=255)
                                                                                             gtgtTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGAGc  >  1:1849039/1‑100 (MQ=255)
                                                                                                 |                                                                                               
GTCTGGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTCGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGAGC  >  NC_000913/2600191‑2600383

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