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,3120GA98.2% 577.6 / ‑2.0 166M68I (ATG→ATAfrwCPTS system fructose‑like enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (81/82);  total (83/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01

Only 100 of 172 total aligned reads displayed.

TGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAG  >  NC_000913/4140221‑4140401
                                                                                           |                                                                                         
      aTGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCtt                                                                                    >  1:3528557/1‑93 (MQ=255)
        ggCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTc                                                                                  >  2:4515281/1‑93 (MQ=255)
           aaaGGCGCAGTGCGGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTc                                                                               <  1:816757/93‑1 (MQ=255)
           aaaGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTc                                                                               <  2:332540/93‑1 (MQ=255)
             aGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGc                                                                             <  2:4150725/93‑1 (MQ=255)
             aGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGc                                                                             >  1:3709715/1‑93 (MQ=255)
              ggCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCa                                                                            <  1:2881276/93‑1 (MQ=255)
                 gcAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCATCg                                                                         >  2:2025527/1‑93 (MQ=255)
                 gcAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCg                                                                         >  1:3147540/1‑93 (MQ=255)
                 gcAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCg                                                                         >  2:1527968/1‑93 (MQ=255)
                   aGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTa                                                                       <  2:847656/93‑1 (MQ=255)
                   aGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTa                                                                       >  1:2094280/1‑93 (MQ=255)
                     tGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACa                                                                     <  2:113622/93‑1 (MQ=255)
                     tGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACa                                                                     >  1:4528391/1‑93 (MQ=255)
                      gCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACAt                                                                    <  2:3674694/93‑1 (MQ=255)
                       ccGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATc                                                                   <  1:3071446/93‑1 (MQ=255)
                        cGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCg                                                                  >  1:1746303/1‑93 (MQ=255)
                           aggCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGtt                                                               <  1:641028/91‑1 (MQ=255)
                           aTGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGtt                                                               <  1:3891335/93‑1 (MQ=255)
                           aTGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGtt                                                               >  2:2649895/1‑93 (MQ=255)
                           aTGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGtt                                                               <  1:3036949/93‑1 (MQ=255)
                            tGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTa                                                              <  1:4161235/93‑1 (MQ=255)
                            tGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTa                                                              >  1:3097919/1‑93 (MQ=255)
                             gCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTAt                                                             <  1:4482510/93‑1 (MQ=255)
                              ccGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTAtt                                                            <  2:42912/93‑1 (MQ=255)
                              ccGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTAtt                                                            <  1:2199351/93‑1 (MQ=255)
                               cGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTc                                                           <  2:1232973/93‑1 (MQ=255)
                                gTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTcc                                                          <  1:2539218/93‑1 (MQ=255)
                                 tAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCa                                                         <  1:3630529/93‑1 (MQ=255)
                                   gCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTTATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCAtt                                                       >  2:2488409/1‑93 (MQ=255)
                                   gCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCAtt                                                       <  2:3555084/93‑1 (MQ=255)
                                   gCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCAtt                                                       <  2:3450375/93‑1 (MQ=255)
                                   gCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCAtt                                                       <  2:992964/93‑1 (MQ=255)
                                   gCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCAtt                                                       <  1:3477753/93‑1 (MQ=255)
                                   gCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCAtt                                                       >  1:2609140/1‑93 (MQ=255)
                                     cGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGc                                                     <  2:3726742/93‑1 (MQ=255)
                                      gATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCa                                                    <  2:2744928/93‑1 (MQ=255)
                                      gATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCa                                                    >  1:3967651/1‑93 (MQ=255)
                                      gATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCa                                                    >  2:2542741/1‑93 (MQ=255)
                                      gATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCa                                                    <  2:3405095/93‑1 (MQ=255)
                                       aTCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCag                                                   >  2:1161976/1‑93 (MQ=255)
                                       aTCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCag                                                   >  1:1669744/1‑93 (MQ=255)
                                       aTCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCag                                                   >  1:1455381/1‑93 (MQ=255)
                                        tCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCaga                                                  <  1:1210865/93‑1 (MQ=255)
                                        tCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCaga                                                  <  2:1833636/93‑1 (MQ=255)
                                        tCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCaga                                                  >  1:2751867/1‑93 (MQ=255)
                                        tCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCaga                                                  >  2:243291/1‑93 (MQ=255)
                                        tCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCaga                                                  >  2:2509061/1‑93 (MQ=255)
                                        tCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTCGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCaga                                                  >  2:1947353/1‑93 (MQ=255)
                                          cAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGc                                                >  1:4397985/1‑93 (MQ=255)
                                           aaaTCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGAAGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCg                                               >  2:4175727/1‑93 (MQ=255)
                                             aTCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGtt                                             >  1:4155576/1‑93 (MQ=255)
                                              tCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTc                                            <  2:28047/93‑1 (MQ=255)
                                              tCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTc                                            >  1:1613041/1‑93 (MQ=255)
                                               cTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCt                                           <  1:2334726/93‑1 (MQ=255)
                                                 gAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTgc                                         >  1:4150782/1‑93 (MQ=255)
                                                 gAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTgc                                         <  1:2462216/93‑1 (MQ=255)
                                                 gAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTgc                                         <  2:3763985/93‑1 (MQ=255)
                                                   aaaaaCTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTgcgc                                       >  2:812720/1‑93 (MQ=255)
                                                    aaaaCTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCt                                      <  1:1464223/93‑1 (MQ=255)
                                                    aaaaCTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCt                                      <  2:1874230/93‑1 (MQ=255)
                                                    aaaaCTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCt                                      >  2:3847777/1‑93 (MQ=255)
                                                     aaaCTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTg                                     >  2:2586136/1‑93 (MQ=255)
                                                      aaCTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTgg                                    >  1:916893/1‑93 (MQ=255)
                                                       aCTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGc                                   >  1:4747566/1‑93 (MQ=255)
                                                        cTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCt                                  >  2:1918100/1‑93 (MQ=255)
                                                        cTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCt                                  >  1:793787/1‑93 (MQ=255)
                                                         tGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTc                                 >  1:1081523/1‑93 (MQ=255)
                                                         tGTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTc                                 <  1:2052552/93‑1 (MQ=255)
                                                          gTTTGATATCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTcc                                >  2:1114138/1‑93 (MQ=255)
                                                               atatCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTgcg                           >  1:3405593/1‑93 (MQ=255)
                                                                tatCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTgcgc                          >  1:1774173/1‑93 (MQ=255)
                                                                tatCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTgcgc                          <  2:3097292/93‑1 (MQ=255)
                                                                tatCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTgcgc                          >  1:488959/1‑93 (MQ=255)
                                                                  tCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTa                        <  2:2243384/93‑1 (MQ=255)
                                                                  tCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTa                        >  1:959234/1‑93 (MQ=255)
                                                                  tCGGCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTa                        >  2:3446111/1‑93 (MQ=255)
                                                                     gCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCg                     >  2:2116844/1‑93 (MQ=255)
                                                                     gCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCg                     >  2:1098330/1‑93 (MQ=255)
                                                                     gCGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCg                     <  2:1501617/93‑1 (MQ=255)
                                                                      cGTTGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCgg                    <  1:877651/93‑1 (MQ=255)
                                                                        ttGCGGGCTTGCCGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTg                  >  1:1734418/1‑93 (MQ=255)
                                                                        ttGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTg                  >  1:1064832/1‑93 (MQ=255)
                                                                        ttGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTg                  <  2:659511/93‑1 (MQ=255)
                                                                        ttGCGGGCTTGACGCTGATAGTGCCATTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTg                  >  2:4682969/1‑93 (MQ=255)
                                                                         tGCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGc                 <  2:2915523/93‑1 (MQ=255)
                                                                          gCGGGCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGcc                >  2:436707/1‑93 (MQ=255)
                                                                            gggCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTg              <  2:2697028/93‑1 (MQ=255)
                                                                            gggCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTg              <  2:1005197/93‑1 (MQ=255)
                                                                             ggCTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTgg             <  2:4253796/93‑1 (MQ=255)
                                                                              tcTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATTGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGCGCGCTATCGGTGCCTggg            >  1:688268/2‑93 (MQ=255)
                                                                               cTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGt           <  2:4528644/93‑1 (MQ=255)
                                                                               cTTGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGt           <  1:3096744/93‑1 (MQ=255)
                                                                                 tGACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTg         >  2:4118807/1‑93 (MQ=255)
                                                                                  gACGCTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTgg        >  2:800255/1‑93 (MQ=255)
                                                                                      cTGATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAAc    >  2:3328357/1‑93 (MQ=255)
                                                                                       tGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACa   >  2:4039031/1‑93 (MQ=255)
                                                                                       tGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACa   >  1:54740/1‑93 (MQ=255)
                                                                                       tGATAGTGCCTTTCCTCGCGGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACa   <  1:59076/93‑1 (MQ=255)
                                                                                        gATAGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAg  >  1:4303733/1‑93 (MQ=255)
                                                                                           |                                                                                         
TGTTGTATGGCAAAGGCGCAGTGCCGGATGCCGTAGCCGATCCAAATCTGAAAAAACTGTTTGATATCGGCGTTGCGGGCTTGACGCTGATGGTGCCTTTCCTCGCAGCGTACATCGGTTATTCCATTGCAGAGCGTTCTGCGCTGGCTCCGTGCGCTATCGGTGCCTGGGTTGGTAACAG  >  NC_000913/4140221‑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: