Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,140,312 G→A M68I (ATG→ATA frwC → PTS system fructose‑like enzyme IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,140,3120GA99.7% 1276.1 / ‑6.0 352M68I (ATG→ATAfrwCPTS system fructose‑like enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (187/164);  total (187/165)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 365 total aligned reads displayed.

GTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAG  >  NC_000913/4140222‑4140401
                                                                                          |                                                                                         
gttgtATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGt                                                                                         <  2:4457356/93‑1 (MQ=255)
 ttgtATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTg                                                                                        >  1:7334477/1‑93 (MQ=255)
 ttgtATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTg                                                                                        >  2:706032/1‑93 (MQ=255)
 ttgtATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTCACGCTGATAGTg                                                                                        >  2:6972752/1‑93 (MQ=255)
   gtATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGcc                                                                                      <  2:6748729/93‑1 (MQ=255)
   gtATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGcc                                                                                      <  1:5692764/93‑1 (MQ=255)
    tATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCt                                                                                     >  1:114843/1‑93 (MQ=255)
      tGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCttt                                                                                   >  2:6428650/1‑93 (MQ=255)
      tGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCttt                                                                                   >  1:7756152/1‑93 (MQ=255)
         cAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCt                                                                                <  2:7376946/93‑1 (MQ=255)
         cAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCt                                                                                >  1:460053/1‑93 (MQ=255)
          aaaGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTc                                                                               >  2:3903837/1‑93 (MQ=255)
           aaGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCg                                                                              >  2:6024310/1‑93 (MQ=255)
             ggCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCa                                                                            <  1:321993/93‑1 (MQ=255)
             ggCGCAGTGCCGGATGCCGTAGCCGATCCAAATCAGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCa                                                                            <  1:4833991/93‑1 (MQ=255)
                 cAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGt                                                                        <  1:835897/93‑1 (MQ=255)
                  aGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTa                                                                       <  2:1378582/93‑1 (MQ=255)
                   gTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTAc                                                                      >  2:2205138/1‑93 (MQ=255)
                    tGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACa                                                                     >  2:1138606/1‑93 (MQ=255)
                      ccGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATc                                                                   <  2:5247141/93‑1 (MQ=255)
                      ccGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATc                                                                   <  2:3706672/93‑1 (MQ=255)
                       cGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCg                                                                  <  1:6980150/93‑1 (MQ=255)
                        ggATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCgg                                                                 >  1:2089939/1‑93 (MQ=255)
                        ggATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCgg                                                                 <  1:1912214/93‑1 (MQ=255)
                         gATGCCGTAGCCGATCCAAATCTGAAAAATCTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGt                                                                <  2:1731540/93‑1 (MQ=255)
                         gATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCGGCGTACATCGGt                                                                >  2:3116/1‑93 (MQ=255)
                         gATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGt                                                                >  1:4795455/1‑93 (MQ=255)
                         gATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGt                                                                >  1:1194323/1‑93 (MQ=255)
                         gATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGt                                                                >  2:2627280/1‑93 (MQ=255)
                           tGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTa                                                              <  2:183722/93‑1 (MQ=255)
                               gTAGCCGATCCAAATCTGAAAAAACTGTTTGATATGGGCGTTGCGGGCTTGGCGCTGATAGTGCCTTTCCTCGCACCCCACATCGGTTATTcg                                                          >  1:4968852/1‑92 (MQ=255)
                               gTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTcc                                                          >  2:665285/1‑93 (MQ=255)
                               gTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTcc                                                          >  1:5908634/1‑93 (MQ=255)
                                tAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCa                                                         >  2:5507890/1‑93 (MQ=255)
                                tAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCa                                                         <  1:6796779/93‑1 (MQ=255)
                                 aGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCAt                                                        <  1:6640224/93‑1 (MQ=255)
                                  gCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCAtt                                                       <  2:7550296/93‑1 (MQ=255)
                                  gCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTCTCCTCGCAGCGTACATCGGTTATTCCAtt                                                       >  2:6279807/1‑93 (MQ=255)
                                    cGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGc                                                     >  1:3605327/1‑93 (MQ=255)
                                     gATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCa                                                    <  2:2886458/93‑1 (MQ=255)
                                     gATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCa                                                    >  2:3300005/1‑93 (MQ=255)
                                      aTCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCag                                                   >  2:2731019/1‑93 (MQ=255)
                                      aTCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCag                                                   >  1:5037896/1‑93 (MQ=255)
                                      aTCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCag                                                   >  1:2526647/1‑93 (MQ=255)
                                      aTCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCCGTTATTCCATTGCag                                                   >  1:3512236/1‑93 (MQ=255)
                                        ccAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCagag                                                 <  1:6326547/93‑1 (MQ=255)
                                        ccAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCagag                                                 >  2:2546407/1‑93 (MQ=255)
                                         cAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGc                                                >  1:4193192/1‑93 (MQ=255)
                                         cAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGc                                                <  2:2987135/93‑1 (MQ=255)
                                           aaTCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGt                                              <  2:4818616/93‑1 (MQ=255)
                                           aaTCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGt                                              <  1:5136088/93‑1 (MQ=255)
                                            aTCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGtt                                             <  2:6319366/93‑1 (MQ=255)
                                            aTCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGtt                                             >  2:3207855/1‑93 (MQ=255)
                                             tCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTc                                            >  1:4625136/1‑93 (MQ=255)
                                              cTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCt                                           >  2:7710909/1‑93 (MQ=255)
                                               tGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTg                                          >  2:657941/1‑93 (MQ=255)
                                                gAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTgc                                         <  2:6255432/93‑1 (MQ=255)
                                                 aaaaaaCTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTgcg                                        <  2:3798530/93‑1 (MQ=255)
                                                    aaaCTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTg                                     <  1:2930506/93‑1 (MQ=255)
                                                     aaCTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTgg                                    >  1:4020952/1‑93 (MQ=255)
                                                     aaCTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGGACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTgg                                    >  1:4397835/1‑93 (MQ=255)
                                                      aCTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGc                                   >  1:6921213/1‑93 (MQ=255)
                                                       cTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCt                                  <  2:7932973/93‑1 (MQ=255)
                                                        tGTTTGATATCAGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGAAGAGCGTTCTGCGCTGGCTc                                 >  2:2312210/1‑93 (MQ=255)
                                                          tttGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCg                               >  1:5620024/1‑93 (MQ=255)
                                                           ttGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGt                              >  1:3794509/1‑93 (MQ=255)
                                                            tGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTg                             >  1:3348595/1‑93 (MQ=255)
                                                               tatCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTgcgc                          >  1:1452419/1‑93 (MQ=255)
                                                               tatCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTgcgc                          <  2:7254867/93‑1 (MQ=255)
                                                                 tCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTa                        <  1:5362790/93‑1 (MQ=255)
                                                                  cGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTAt                       <  2:1512664/93‑1 (MQ=255)
                                                                  cGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTAt                       <  2:3561688/93‑1 (MQ=255)
                                                                  cGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTAt                       <  2:5795540/93‑1 (MQ=255)
                                                                  cGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTAt                       >  1:4091111/1‑93 (MQ=255)
                                                                    gCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCg                     <  2:5668650/93‑1 (MQ=255)
                                                                    gCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCg                     >  2:4999432/1‑93 (MQ=255)
                                                                    gCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCGCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCg                     <  1:7641569/93‑1 (MQ=255)
                                                                      gTTGCGGGCTTGGCGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGt                   <  1:5512475/93‑1 (MQ=255)
                                                                       ttGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTg                  >  1:2723652/1‑93 (MQ=255)
                                                                         gCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGcc                >  1:7903095/1‑93 (MQ=255)
                                                                          cGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCt               <  1:3038200/93‑1 (MQ=255)
                                                                          cGGGCTTGACGCTGATAGTGCCATTCCTGGCAGCGTAAATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCTGGGCCt               <  2:422612/93‑1 (MQ=255)
                                                                            ggCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTgg             <  2:1202763/93‑1 (MQ=255)
                                                                              cTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGt           >  2:4630800/1‑93 (MQ=255)
                                                                              cTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGt           >  1:26490/1‑93 (MQ=255)
                                                                               ttGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGtt          >  2:2116019/1‑93 (MQ=255)
                                                                               ttGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGtt          <  1:101889/93‑1 (MQ=255)
                                                                               ttGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGtt          >  1:6507796/1‑93 (MQ=255)
                                                                               ttGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGtt          >  2:4073235/1‑93 (MQ=255)
                                                                                tGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTg         >  1:7486654/1‑93 (MQ=255)
                                                                                 gACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTgg        <  2:5387700/93‑1 (MQ=255)
                                                                                    gCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTaa     >  1:6049360/1‑93 (MQ=255)
                                                                                    gCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTaa     <  2:6832934/93‑1 (MQ=255)
                                                                                    gCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTaa     >  2:1240762/1‑93 (MQ=255)
                                                                                    gCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTaa     >  2:3452578/1‑93 (MQ=255)
                                                                                     cTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGaaaa    >  2:4411201/1‑92 (MQ=255)
                                                                                     cTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAAc    >  1:6242914/1‑93 (MQ=255)
                                                                                      tGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACa   >  1:4497048/1‑93 (MQ=255)
                                                                                      tGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACa   <  1:159573/93‑1 (MQ=255)
                                                                                       gATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAg  >  2:1584828/1‑93 (MQ=255)
                                                                                          |                                                                                         
GTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAG  >  NC_000913/4140222‑4140401

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