Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,600,288 C→T 100% 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,2880CT98.9% 291.2 / ‑4.6 87S128L (TCG→TTG) gcvRtranscriptional repressor, regulatory protein accessory to GcvA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (39/47);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TCTGGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTCGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGAGC  >  NC_000913/2600192‑2600383
                                                                                                |                                                                                               
tctgGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTtgcg                                                                                              >  1:2066991/1‑100 (MQ=255)
  tgGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgca                                                                                            >  1:2164813/1‑100 (MQ=255)
    ggTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcacg                                                                                          <  2:527062/100‑1 (MQ=255)
     gTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcacgc                                                                                         >  2:2080973/1‑100 (MQ=255)
      ttCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTtgcgtgcgca                                                                                        <  1:1250663/100‑7 (MQ=255)
      ttCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcacgca                                                                                        <  2:558614/100‑1 (MQ=255)
      ttCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGcgcacgca                                                                                        >  1:501605/1‑100 (MQ=255)
         aGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAAcc                                                                                     <  1:160047/100‑1 (MQ=255)
         aGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAAcc                                                                                     <  2:519166/100‑1 (MQ=255)
          ggTCGATTTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGTACATTGCGGAGTTGGTGTTGCGCACGCAAcc                                                                                     <  2:585572/99‑1 (MQ=255)
          ggTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGTGCGCAAcac                                                                                    <  1:1564356/100‑3 (MQ=255)
          ggTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACct                                                                                    >  2:1711427/1‑100 (MQ=255)
             cGATGTGGCAGACTCCCCGCATTTAATTGATCGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACctgct                                                                                 <  1:584630/100‑1 (MQ=255)
             cGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGTGCGCAAcaccag                                                                                 <  2:203831/100‑6 (MQ=255)
             cGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACctgct                                                                                 <  2:38906/100‑1 (MQ=255)
             cGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACctgct                                                                                 <  1:1468167/100‑1 (MQ=255)
                tgtgGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGaa                                                                              <  1:573516/100‑1 (MQ=255)
                 gtgGCAGACTCCCCGCATTTAATTTAACGCTTCACAGCACTTGTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGTGCGCAAcaccagc                                                                                <  2:1979911/97‑7 (MQ=255)
                 gtgGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGaaa                                                                             <  2:546457/100‑1 (MQ=255)
                  tgGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGaaaa                                                                            >  2:1815619/1‑100 (MQ=255)
                  tgGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGaaaa                                                                            >  1:1887225/1‑100 (MQ=255)
                   ggCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAAt                                                                           >  1:335036/1‑100 (MQ=255)
                   ggCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAAt                                                                           >  1:2194518/1‑100 (MQ=255)
                   ggCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAAt                                                                           >  1:2386932/1‑100 (MQ=255)
                    gCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAAt                                                                           <  1:156724/99‑1 (MQ=255)
                     cAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGa                                                                         >  1:1701341/1‑100 (MQ=255)
                     cAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGa                                                                         >  1:1347496/1‑100 (MQ=255)
                        aCTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACg                                                                      <  1:1895006/100‑1 (MQ=255)
                        aCTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACg                                                                      <  2:1947720/100‑1 (MQ=255)
                         cTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACg                                                                      >  1:1332733/1‑99 (MQ=255)
                           ccccGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGc                                                                   >  1:1931502/1‑100 (MQ=255)
                            cccGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCt                                                                  >  1:1882273/1‑100 (MQ=255)
                            cccGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGATGTTGCGCACGCAAACTGCTGAAAATGAAc                                                                       >  1:2004918/1‑95 (MQ=255)
                                   ttAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGt                                                           <  1:2304368/100‑1 (MQ=255)
                                    tAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGtt                                                          <  1:2291769/100‑1 (MQ=255)
                                    tAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGtt                                                          >  1:287361/1‑100 (MQ=255)
                                     aaTTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTg                                                         <  1:1451240/100‑1 (MQ=255)
                                     aaTTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTg                                                         >  1:1190352/1‑100 (MQ=255)
                                      aTTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCCACCTGCTGAAAATGAACGGGCTGCGCAGTTGc                                                        >  1:508486/1‑100 (MQ=255)
                                      aTTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTg                                                         <  1:446144/99‑1 (MQ=255)
                                      aTTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGc                                                        <  2:1600230/100‑1 (MQ=255)
                                      aTTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGc                                                        >  1:1024718/1‑100 (MQ=255)
                                       ttGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCa                                                       >  2:2359236/1‑100 (MQ=255)
                                       ttGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCa                                                       >  1:1911161/1‑100 (MQ=255)
                                        tGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCat                                                      >  2:2203165/1‑100 (MQ=255)
                                          aaCGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGCACGGGCTGCGCAGTTGCatat                                                    >  1:1906618/1‑100 (MQ=255)
                                           aCGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATAtt                                                   >  1:857579/1‑100 (MQ=255)
                                             gCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCa                                                 <  1:1644449/100‑1 (MQ=255)
                                             gCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCa                                                 <  1:101882/100‑1 (MQ=255)
                                              cTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAg                                                >  2:1203542/1‑100 (MQ=255)
                                                 cacaGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATa                                             >  2:1805249/1‑100 (MQ=255)
                                                  acaGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTCGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATaa                                            >  1:1402657/1‑100 (MQ=255)
                                                      cACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCg                                         >  2:1618306/1‑99 (MQ=255)
                                                       aCTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGcc                                       >  1:1043628/1‑100 (MQ=255)
                                                         ttttgaactcgcATCATATGAACATTGTGTAGCTGGTGTTGCGCACGCAACTTGCTGAAAATGAACGGTTTGCTTAGTTGCATAGTCAGATAAc                                           <  2:1807620/85‑1 (MQ=255)
                                                         ttttCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCca                                     <  2:261852/100‑1 (MQ=255)
                                                          tttCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCcac                                    <  1:386864/100‑1 (MQ=255)
                                                          tttCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCcac                                    <  2:1190840/100‑1 (MQ=255)
                                                           ttCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCcaca                                   <  2:1698123/100‑1 (MQ=255)
                                                           ttCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCcaca                                   <  2:129998/100‑1 (MQ=255)
                                                            tCGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAg                                  >  1:730783/1‑100 (MQ=255)
                                                             cGACGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGc                                 >  1:1943729/1‑100 (MQ=255)
                                                               aCGCGCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGccc                               >  1:894693/1‑100 (MQ=255)
                                                                cgcgCATCATATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGcccc                              >  1:1146014/1‑100 (MQ=255)
                                                                     atcatATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATc                         <  1:1105605/100‑1 (MQ=255)
                                                                      tcatATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCt                        >  1:1210552/1‑100 (MQ=255)
                                                                       catATGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCt                        >  2:913222/1‑99 (MQ=255)
                                                                        atatGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGc                      <  1:1771782/100‑1 (MQ=255)
                                                                        atatGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGc                      <  2:641301/100‑1 (MQ=255)
                                                                         tatGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCg                     <  1:1244120/100‑1 (MQ=255)
                                                                         tatGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCg                     >  1:241022/1‑100 (MQ=255)
                                                                          atGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCgg                    <  2:501602/100‑1 (MQ=255)
                                                                           tGAACATTGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGa                   >  1:163238/1‑100 (MQ=255)
                                                                             aaCATTGCGGAGCTGGTGTTGCACACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACg                 <  1:1270534/100‑1 (MQ=255)
                                                                             aaCATTGCAACGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGTCTACTGCACCGCATCTGCGGACg                 <  2:175745/100‑1 (MQ=255)
                                                                                 ttGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACgcagc             <  1:1820411/100‑1 (MQ=255)
                                                                                  tGCGGAGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACgcagca            >  1:378097/1‑100 (MQ=255)
                                                                                       aGCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAatat       <  1:959030/100‑1 (MQ=255)
                                                                                        gCTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATAtt      <  1:2402684/100‑1 (MQ=255)
                                                                                         cTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTg     <  2:93139/100‑1 (MQ=255)
                                                                                         cTGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTg     <  2:1089489/100‑1 (MQ=255)
                                                                                          tGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGa    >  1:2053916/1‑100 (MQ=255)
                                                                                          tGGTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGa    >  1:24576/1‑100 (MQ=255)
                                                                                           ggTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGAg   <  1:1234691/100‑1 (MQ=255)
                                                                                           ggTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGAg   <  1:1713519/100‑1 (MQ=255)
                                                                                           ggTGTTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGAg   <  2:58246/100‑1 (MQ=255)
                                                                                            gtgtTGCGCACGCAACCTGCTGAAAATGAACGGGTTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGAGc  <  1:777579/100‑1 (MQ=255)
                                                                                            gtgtTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGAGc  <  2:1256216/100‑1 (MQ=255)
                                                                                            gtgtTGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGAGc  <  1:1785296/100‑1 (MQ=255)
                                                                                                |                                                                                               
TCTGGGTTCAGGTCGATGTGGCAGACTCCCCGCATTTAATTGAACGCTTCACAGCACTTTTCGACGCGCATCATATGAACATTGCGGAGCTGGTGTCGCGCACGCAACCTGCTGAAAATGAACGGGCTGCGCAGTTGCATATTCAGATAACCGCCCACAGCCCCGCATCTGCGGACGCAGCAAATATTGAGC  >  NC_000913/2600192‑2600383

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