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,2880CT96.3% 276.6 / ‑3.9 84S128L (TCG→TTG) gcvRtranscriptional repressor, regulatory protein accessory to GcvA
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (41/38);  total (43/41)
Fisher's exact test for biased strand distribution p-value = 6.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCTGGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTCGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGAG  >  NC_000913/2600191‑2600382
                                                                                                 |                                                                                              
gtctgGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTtgc                                                                                              <  2:1899110/100‑1 (MQ=255)
 tctgGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTtgcg                                                                                             <  1:1944622/100‑1 (MQ=255)
 cttccgatctggTCGGTGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTtgcg                                                                                             <  2:1843140/90‑1 (MQ=255)
   tgGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgca                                                                                           >  1:1153003/1‑100 (MQ=255)
   tgGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgca                                                                                           <  2:2096009/100‑1 (MQ=255)
   tgGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgca                                                                                           <  1:797303/100‑1 (MQ=255)
    gggTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcac                                                                                          >  2:932923/1‑100 (MQ=255)
    gggTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcac                                                                                          >  1:652772/1‑100 (MQ=255)
    gggTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcac                                                                                          >  1:1841104/1‑100 (MQ=255)
     ggTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcacg                                                                                         <  1:1515123/100‑1 (MQ=255)
     ggTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcacg                                                                                         <  2:1309884/100‑1 (MQ=255)
     ggTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcacg                                                                                         <  1:1711577/100‑1 (MQ=255)
     ggTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcacg                                                                                         <  1:944033/100‑1 (MQ=255)
     ggTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcacg                                                                                         <  1:1860637/100‑1 (MQ=255)
       ttCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcacgca                                                                                       <  1:545919/100‑1 (MQ=255)
          aGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAAcc                                                                                    >  2:1666338/1‑100 (MQ=255)
           ggTCGGTGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTtgcgagatcggaag                                                                                   >  1:1843145/1‑90 (MQ=255)
           ggTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACct                                                                                   >  2:1473698/1‑100 (MQ=255)
           ggTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACct                                                                                   >  1:1805087/1‑100 (MQ=255)
           ggTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACct                                                                                   >  2:1710968/1‑100 (MQ=255)
             tCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACctgc                                                                                 <  2:205511/100‑1 (MQ=255)
              cGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACctgct                                                                                >  1:1242968/1‑100 (MQ=255)
                 tgtgGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGaa                                                                             >  1:2004835/1‑100 (MQ=255)
                   tgGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGaaaa                                                                           <  1:1864196/100‑1 (MQ=255)
                   tgGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGaaaa                                                                           >  2:1177840/1‑100 (MQ=255)
                    ggCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAAt                                                                          >  1:46723/1‑100 (MQ=255)
                     gCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATg                                                                         >  1:2050193/1‑100 (MQ=255)
                     gCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATg                                                                         <  1:858976/100‑1 (MQ=255)
                      cAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGa                                                                        >  2:891240/1‑100 (MQ=255)
                       aGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGaa                                                                       <  2:328771/100‑1 (MQ=255)
                            ccccGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGc                                                                  >  1:1278597/1‑100 (MQ=255)
                             cccGCATTTAATTGAACGCTTCACAGCACTTTTCGCCGCGCATCATATGAACATTGCTGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCt                                                                 >  2:1644042/1‑100 (MQ=255)
                             cccGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCt                                                                 >  2:1077390/1‑100 (MQ=255)
                             cccGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTCGCGCACGCAACCTGCTGAAAATGAACGGGCt                                                                 >  2:2151275/1‑100 (MQ=255)
                                 cATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTgcgc                                                             <  1:1666343/100‑1 (MQ=255)
                                    ttAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTAGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGt                                                          >  2:1277930/1‑100 (MQ=255)
                                     tAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGtt                                                         <  1:1473702/100‑1 (MQ=255)
                                     tAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGtt                                                         >  2:1318891/1‑100 (MQ=255)
                                      aaTTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTg                                                        >  1:1328598/1‑100 (MQ=255)
                                       aTTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGc                                                       <  1:1318894/100‑1 (MQ=255)
                                            aCGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATAtt                                                  <  2:1046989/100‑1 (MQ=255)
                                             cGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTc                                                 <  1:932925/100‑1 (MQ=255)
                                             cGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTc                                                 <  2:1617346/100‑1 (MQ=255)
                                              gCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCa                                                <  2:187577/100‑1 (MQ=255)
                                              gCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCa                                                <  1:1592532/100‑1 (MQ=255)
                                              gCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCa                                                <  1:1411191/100‑1 (MQ=255)
                                              gCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCa                                                >  1:1323282/1‑100 (MQ=255)
                                              gCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCa                                                <  1:1960245/100‑1 (MQ=255)
                                              gCTTCACAGCACTTTTCGACGCGCATCATATGAACATTCCGGAGCTGGTGTTGCGCACTCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCa                                                <  2:1278594/100‑1 (MQ=255)
                                               cTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTccg                                               >  1:775825/1‑98 (MQ=255)
                                                ttCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGa                                              <  1:786648/100‑1 (MQ=255)
                                                 tCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGAt                                             >  1:859560/1‑100 (MQ=255)
                                                 tCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGAt                                             <  2:652771/100‑1 (MQ=255)
                                                 tCACAGCACTTGTCGACGCGCATCATAGGAACATTGCGGAGCTGGTGTCGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAg                                               <  1:2183772/98‑1 (MQ=255)
                                                  cacaGCACTTTTCAACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATa                                            >  1:1587360/1‑100 (MQ=255)
                                                   acaGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATaa                                           <  2:1323279/100‑1 (MQ=255)
                                                     aGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCGGATAAcc                                         >  2:2251744/1‑100 (MQ=255)
                                                        aCTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGcc                                      <  2:1678969/100‑1 (MQ=255)
                                                         cTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGccc                                     >  1:1491792/1‑100 (MQ=255)
                                                             tCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAg                                 >  2:2267699/1‑100 (MQ=255)
                                                              cGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGc                                <  1:2155895/100‑1 (MQ=255)
                                                               gACGCGCATCATATGAACATTGCGGAGCTGGTGTGGCGCACGCAACCGGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCca                                    <  2:1914920/95‑1 (MQ=255)
                                                                aCGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGccc                              <  1:1983370/100‑1 (MQ=255)
                                                                 cgcgCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGACCGGGATGCGCAGTTGCATATTCAGATAACCGCCCACAGcccc                             >  2:372236/1‑100 (MQ=255)
                                                                 cgcgCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGcccc                             >  1:2161708/1‑100 (MQ=255)
                                                                  tcttccgatctGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCg                            <  2:830338/90‑1 (MQ=255)
                                                                  gcgcATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCg                            <  2:1607073/100‑1 (MQ=255)
                                                                    gcATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCa                          >  2:210617/1‑100 (MQ=255)
                                                                     catcatATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCAt                         >  1:1319957/1‑100 (MQ=255)
                                                                         tcttGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGc                     <  1:1051609/97‑1 (MQ=255)
                                                                         atatGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGc                     <  1:2212484/100‑1 (MQ=255)
                                                                          tatGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAGCGGCCCACAGCGCCGCATTTGCg                    <  2:465331/100‑1 (MQ=255)
                                                                           atGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCgg                   >  1:1565240/1‑100 (MQ=255)
                                                                            tGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGcaga                  >  2:1051607/1‑97 (MQ=255)
                                                                            tGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGa                  <  1:1710973/100‑1 (MQ=255)
                                                                             gAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTTCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGa                  >  1:690567/1‑99 (MQ=255)
                                                                             gAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGAc                 <  1:375106/100‑1 (MQ=255)
                                                                              aaCATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACg                >  2:859855/1‑100 (MQ=255)
                                                                               aCATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACgc               >  1:1873836/1‑100 (MQ=255)
                                                                               aCATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACgc               >  2:880956/1‑100 (MQ=255)
                                                                                cATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACgca              >  1:832197/1‑100 (MQ=255)
                                                                                     cGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCaaa         >  2:243546/1‑100 (MQ=255)
                                                                                        aGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAatat      >  1:901536/1‑100 (MQ=255)
                                                                                         tcTGGTGTTGCGCACGCAACCTGCTGAAAATGAATGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATAtt     <  2:1238306/99‑1 (MQ=255)
                                                                                         gCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATAtt     <  2:1587355/100‑1 (MQ=255)
                                                                                         gCTGGTGTCGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACATCCCCGCATCTGCGGACGCAGCAAATAtt     <  2:1767764/100‑1 (MQ=255)
                                                                                           tGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGa   >  2:2009337/1‑100 (MQ=255)
                                                                                            ggTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATAttaga  >  1:1238309/1‑97 (MQ=255)
                                                                                                 |                                                                                              
GTCTGGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTCGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGAG  >  NC_000913/2600191‑2600382

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: