Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,536,171 G→T R113L (CGC→CTC)  crr → glucose‑specific enzyme IIA component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,536,1710GT96.2% 266.0 / ‑1.0 81R113L (CGC→CTC) crrglucose‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (39/38);  total (40/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 = 1.00e+00

CGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGG  >  NC_000913/2536073‑2536266
                                                                                                  |                                                                                               
cGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAgcgc                                                                                                >  1:1429703/1‑100 (MQ=255)
 gAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGctcg                                                                                               <  2:614645/100‑1 (MQ=255)
      tGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGaaa                                                                                          >  1:1049532/1‑100 (MQ=255)
       naTAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAg                                                                                         >  2:788428/2‑100 (MQ=255)
       gATAGCGGCGTTGAACTGTTCGTCCACTTTTGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTAAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAg                                                                                         <  2:563297/100‑1 (MQ=255)
       gATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAg                                                                                         <  1:110378/100‑1 (MQ=255)
           gcggcgTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTgg                                                                                     >  1:657831/1‑100 (MQ=255)
           gcggcgTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTgg                                                                                     >  1:801517/1‑100 (MQ=255)
             ggcgTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCg                                                                                   <  1:884445/100‑1 (MQ=255)
              gcgTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGa                                                                                  >  1:796953/1‑100 (MQ=255)
                   gAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTg                                                                             >  1:689113/1‑100 (MQ=255)
                   gAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAACCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTg                                                                             >  2:1009720/1‑100 (MQ=255)
                   gAAATGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTg                                                                             <  2:317358/100‑1 (MQ=255)
                    aaCTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGt                                                                            <  2:600643/100‑1 (MQ=255)
                    aaCTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGt                                                                            <  2:1089227/100‑1 (MQ=255)
                       tGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCAt                                                                         >  1:174532/1‑100 (MQ=255)
                        gTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCAtt                                                                        >  1:1171317/1‑100 (MQ=255)
                         ttCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTg                                                                       >  1:1264821/1‑100 (MQ=255)
                         ttCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTg                                                                       <  1:1146554/100‑1 (MQ=255)
                         ttCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTCAAAGTTGGCGATACTGTCAtt                                                                        >  1:1265080/1‑99 (MQ=255)
                               cACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTg                                                                 <  1:1270053/100‑1 (MQ=255)
                                 cTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGa                                                                <  1:1050135/99‑1 (MQ=255)
                                  ttCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATc                                                              <  2:226275/100‑1 (MQ=255)
                                    cGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTg                                                            >  1:119330/1‑100 (MQ=255)
                                     ggTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGc                                                           >  2:199673/1‑100 (MQ=255)
                                      gTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGcc                                                          <  2:1085198/100‑1 (MQ=255)
                                       tATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCg                                                         <  2:101862/100‑1 (MQ=255)
                                       tATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCg                                                         <  2:102607/100‑1 (MQ=255)
                                       tATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCg                                                         <  1:274083/100‑1 (MQ=255)
                                       tATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCg                                                         <  1:74721/100‑1 (MQ=255)
                                          cGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTTAAGAAGGTAAGCTCGTGTAAGTTGGCGATACTGTCATTGAATTTGATCTGCCgctg                                                      >  2:723257/1‑100 (MQ=255)
                                          cGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCgctg                                                      >  2:1518697/1‑100 (MQ=255)
                                            acacCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAACTTGGCAATACTGTCGTTGAATTTGATCTGCCgctgct                                                    >  2:1064209/1‑100 (MQ=255)
                                             cacCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGctgctg                                                   <  2:119330/100‑1 (MQ=255)
                                               ccGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGa                                                 >  1:607918/1‑100 (MQ=255)
                                               ccGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGa                                                 >  1:731342/1‑100 (MQ=255)
                                                 gTTGAAGTGAATGGCGAGGGTTTCAAGCGTATGGCTGAGGAAGGTCAGCGCGGGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGa                                                 <  2:679628/98‑1 (MQ=255)
                                                  ttGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAaga                                              <  1:1131219/100‑1 (MQ=255)
                                                   tGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAagag                                             <  1:404015/100‑1 (MQ=255)
                                                   tGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAagag                                             >  2:285591/1‑100 (MQ=255)
                                                   tGAACTAGAAGGCGAAGGCTGCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAagag                                             <  2:245254/100‑1 (MQ=255)
                                                    gAACTGAAAGGGGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACGGTCATTGAATTTGATCTGCCGCTGCTGGAAgaga                                            >  2:830343/1‑100 (MQ=255)
                                                    gAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAgaga                                            <  2:562481/100‑1 (MQ=255)
                                                       cTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAg                                         >  1:1536261/1‑100 (MQ=255)
                                                        tGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGc                                        <  2:801512/100‑1 (MQ=255)
                                                         gAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGcc                                       >  2:470817/1‑100 (MQ=255)
                                                         gAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGcc                                       >  1:676224/1‑100 (MQ=255)
                                                         gAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTAAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGcc                                       >  2:846259/1‑100 (MQ=255)
                                                            aGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAg                                    >  1:634260/1‑100 (MQ=255)
                                                             ggCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGt                                   >  1:1477697/1‑100 (MQ=255)
                                                             ggCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGt                                   >  1:1266851/1‑100 (MQ=255)
                                                              gCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTc                                  <  1:1076447/100‑1 (MQ=255)
                                                                  aGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTAcc                              <  1:313424/100‑1 (MQ=255)
                                                                   ggCTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTAccc                             <  2:760833/100‑1 (MQ=255)
                                                                     cTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTg                           >  1:473173/1‑100 (MQ=255)
                                                                     cTTCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTg                           <  2:174532/100‑1 (MQ=255)
                                                                     cTTCAAGCGTATTGATGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTg                           <  2:1580730/100‑1 (MQ=255)
                                                                       tCAAGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGAc                         <  1:752617/100‑1 (MQ=255)
                                                                         aaGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTc                       >  1:894209/1‑100 (MQ=255)
                                                                          aGCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTcc                      >  2:1138196/1‑100 (MQ=255)
                                                                           gCGTATTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCg                     >  1:308988/1‑100 (MQ=255)
                                                                               aTTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGgttg                 <  2:1447362/100‑1 (MQ=255)
                                                                               aTTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGgttg                 <  2:208246/100‑1 (MQ=255)
                                                                               aTTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGgttg                 <  1:231322/100‑1 (MQ=255)
                                                                               aTTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGgttg                 <  2:208164/100‑1 (MQ=255)
                                                                               aTTGCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGgttg                 <  1:169883/100‑1 (MQ=255)
                                                                                 tGCTGAAGAAGGTCAGCCCGTGAAGGTTGGCGACACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAACGCCAAGTCTACCCTGACTCCGgttgt                <  1:305256/99‑1 (MQ=255)
                                                                                  gCTGAAGAAGGTCAGCTCGTGAATGTTGGCGATACTGTCATTGAATTTGATATGCCGCTGCTGGAAGAGAAAGCCAAGTATACCCTGACTCCGGTTGTTa              >  2:624002/1‑100 (MQ=255)
                                                                                  gCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTa              <  2:689111/100‑1 (MQ=255)
                                                                                  gCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTa              >  1:95353/1‑100 (MQ=255)
                                                                                  gCTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTa              <  1:243530/100‑1 (MQ=255)
                                                                                   cTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTAt             >  2:1323873/1‑100 (MQ=255)
                                                                                   cTGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTAt             >  2:941789/1‑100 (MQ=255)
                                                                                    tGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTAtc            >  2:276540/1‑100 (MQ=255)
                                                                                    tGAAGAAGGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTAtc            >  1:362847/1‑100 (MQ=255)
                                                                                        gaagGTAAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCa        <  2:74886/100‑1 (MQ=255)
                                                                                          agGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAAc      <  1:1367296/100‑1 (MQ=255)
                                                                                          agGTCAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAAc      <  1:1366283/100‑1 (MQ=255)
                                                                                              cAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATgg  <  2:556577/100‑1 (MQ=255)
                                                                                              cAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATgg  <  2:1395404/100‑1 (MQ=255)
                                                                                              cAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATgg  >  2:1255395/1‑100 (MQ=255)
                                                                                              cAGCTCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATgg  >  2:497921/1‑100 (MQ=255)
                                                                                              cAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATgg  >  2:1458293/1‑100 (MQ=255)
                                                                                                  |                                                                                               
CGAATCTGATAGCGGCGTTGAACTGTTCGTCCACTTCGGTATCGACACCGTTGAACTGAAAGGCGAAGGCTTCAAGCGTATTGCTGAAGAAGGTCAGCGCGTGAAAGTTGGCGATACTGTCATTGAATTTGATCTGCCGCTGCTGGAAGAGAAAGCCAAGTCTACCCTGACTCCGGTTGTTATCTCCAACATGG  >  NC_000913/2536073‑2536266

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: ‑