Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,965,001 T→C M486V (ATG→GTG)  ptsP ← PtsP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,965,0010TC99.4% 574.8 / ‑4.6 177M486V (ATG→GTG) ptsPPtsP
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (87/89);  total (88/89)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

Only 100 of 185 total aligned reads displayed.

CTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCAT  >  NC_000913/2964911‑2965090
                                                                                          |                                                                                         
ctTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCAccc                                                                                         <  1:970108/93‑1 (MQ=255)
ctTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCAccc                                                                                         >  2:1293666/1‑93 (MQ=255)
ctTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCAccc                                                                                         >  1:2512008/1‑93 (MQ=255)
 ttCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCAcccc                                                                                        >  2:5142/1‑93 (MQ=255)
 ttCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCAcccc                                                                                        >  2:5141/1‑93 (MQ=255)
  tCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCGTTGAACATTTGCAGCACcccc                                                                                       <  2:61490/93‑1 (MQ=255)
   cGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCt                                                                                      <  2:1745302/93‑1 (MQ=255)
   cGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCt                                                                                      <  1:2222188/93‑1 (MQ=255)
   cGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCt                                                                                      <  2:4052226/93‑1 (MQ=255)
     cTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGa                                                                                    <  2:2508024/93‑1 (MQ=255)
       gATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGACCATTTGCAGCATCCCCTGATa                                                                                  >  1:4331551/1‑93 (MQ=255)
       gATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATa                                                                                  >  1:193114/1‑93 (MQ=255)
       gATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATa                                                                                  >  1:305282/1‑93 (MQ=255)
        aTCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATAc                                                                                 <  1:4226344/93‑1 (MQ=255)
            gCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTgcg                                                                             >  2:2877055/1‑93 (MQ=255)
              aTGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGcc                                                                           >  1:1900449/1‑93 (MQ=255)
              aTGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGcc                                                                           <  1:91637/93‑1 (MQ=255)
               tGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGccc                                                                          >  1:2533020/1‑92 (MQ=255)
                gTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCAc                                                                         >  2:2754667/1‑93 (MQ=255)
                gTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCAc                                                                         >  1:3675616/1‑93 (MQ=255)
                gTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCAc                                                                         <  2:2964969/93‑1 (MQ=255)
                 tAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCAcc                                                                        >  2:455905/1‑93 (MQ=255)
                   aGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTg                                                                      <  2:346887/93‑1 (MQ=255)
                   aGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTg                                                                      >  2:1327457/1‑93 (MQ=255)
                   aGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTg                                                                      >  2:245430/1‑93 (MQ=255)
                    ggCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGt                                                                     >  1:4179270/1‑93 (MQ=255)
                    ggCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGt                                                                     >  2:2143041/1‑93 (MQ=255)
                      cAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttc                                                                   <  1:1088280/93‑1 (MQ=255)
                      aaGCTGCTTATCTGTTGCGACATCCAGCGTCCGCAAGGTGACGGGTTTATCATTGACCCTTTGCAGCACCCCCTGATACTGCTGCGCCTGttc                                                                   <  2:4744704/92‑1 (MQ=255)
                       aGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttct                                                                  >  1:2871933/1‑93 (MQ=255)
                       aGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttct                                                                  >  2:3516548/1‑93 (MQ=255)
                        gctgctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttctt                                                                 >  1:4520268/1‑93 (MQ=255)
                        gctgctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttctt                                                                 >  2:2096664/1‑93 (MQ=255)
                         ctgctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttc                                                                <  2:4170867/93‑1 (MQ=255)
                         ctgcgTAGCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTGTATCATTGAACATTTGCAGCTCCCACTGATACTGCGCCCCCTGttcttc                                                                <  2:2884445/93‑1 (MQ=255)
                          tgctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttct                                                               >  1:2978177/1‑93 (MQ=255)
                          tgctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttct                                                               >  1:2734368/1‑93 (MQ=255)
                           gctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttctt                                                              <  1:3707332/93‑1 (MQ=255)
                           gctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttctt                                                              <  1:2050822/93‑1 (MQ=255)
                            ctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttcttc                                                             >  1:3480693/1‑93 (MQ=255)
                            ctTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGttcttcttc                                                             <  1:4661980/93‑1 (MQ=255)
                             ttATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTcc                                                            <  2:1368222/93‑1 (MQ=255)
                                tCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGAc                                                         <  1:2301729/93‑1 (MQ=255)
                                 cTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACg                                                        >  2:2963492/1‑93 (MQ=255)
                                 cTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACg                                                        <  2:855590/93‑1 (MQ=255)
                                  tGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACgg                                                       <  1:3690770/93‑1 (MQ=255)
                                  tGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGAACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACgg                                                       <  2:4082646/93‑1 (MQ=255)
                                   gCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACggg                                                      <  1:2395937/93‑1 (MQ=255)
                                   gCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATATGCAGCACCCCCTGATACTGCGCCACCTGGTCTTCTTCCGACggg                                                      >  1:4216074/1‑93 (MQ=255)
                                    cTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGa                                                     <  1:1366347/93‑1 (MQ=255)
                                       cGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAAc                                                  <  2:3238507/93‑1 (MQ=255)
                                       cGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAAc                                                  <  2:3426085/93‑1 (MQ=255)
                                         aCATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCa                                                >  2:631802/1‑93 (MQ=255)
                                         aCATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCa                                                <  1:326240/93‑1 (MQ=255)
                                         aCATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCa                                                >  1:3592653/1‑93 (MQ=255)
                                          cATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCAc                                               >  1:3982623/1‑93 (MQ=255)
                                          cATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCAc                                               >  1:732639/1‑93 (MQ=255)
                                           aTCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACt                                              >  1:3174332/1‑93 (MQ=255)
                                            tCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACtt                                             >  1:449430/1‑93 (MQ=255)
                                            tCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACtt                                             >  1:3568027/1‑93 (MQ=255)
                                            tCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACtt                                             >  2:445416/1‑93 (MQ=255)
                                               aGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTgc                                          <  1:2251876/93‑1 (MQ=255)
                                               aGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTgc                                          >  1:1808839/1‑93 (MQ=255)
                                                   tACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTgcagca                                      <  1:3622955/93‑1 (MQ=255)
                                                   tACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTgcagca                                      >  2:1218033/1‑93 (MQ=255)
                                                    aCGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCAt                                     >  1:1226688/1‑93 (MQ=255)
                                                       cAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGaa                                  >  2:3985474/1‑93 (MQ=255)
                                                       cAAGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGaa                                  <  2:3763053/93‑1 (MQ=255)
                                                        aaGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAAt                                 >  2:1589935/1‑93 (MQ=255)
                                                        aaGGTGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAAt                                 <  2:3841728/93‑1 (MQ=255)
                                                           gttgcGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATggg                              <  1:1603515/89‑1 (MQ=255)
                                                            tGACGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATgggt                             >  2:3899942/1‑92 (MQ=255)
                                                              aCGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGAtt                           >  1:1645952/1‑93 (MQ=255)
                                                              aCGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGAtt                           <  1:4732214/93‑1 (MQ=255)
                                                              aCGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGAtt                           <  1:4732213/93‑1 (MQ=255)
                                                               cGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGAttt                          >  1:1155062/1‑93 (MQ=255)
                                                               cGGGTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGAttt                          <  1:821664/93‑1 (MQ=255)
                                                                gggTTTATCATTGAACATTTGCATCTCCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTc                         <  2:2595228/93‑1 (MQ=255)
                                                                gggTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTc                         <  1:4007923/93‑1 (MQ=255)
                                                                gggTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTc                         <  1:3926096/93‑1 (MQ=255)
                                                                 ggTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCa                        <  2:1853843/93‑1 (MQ=255)
                                                                  gTTTATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAg                       <  2:1463260/93‑1 (MQ=255)
                                                                   tttATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGt                      <  2:689883/93‑1 (MQ=255)
                                                                   tttATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGt                      <  2:2492353/93‑1 (MQ=255)
                                                                   gttATCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGt                      <  2:2492352/92‑1 (MQ=255)
                                                                      aTCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCg                   >  1:3099145/1‑93 (MQ=255)
                                                                       tCATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGa                  <  2:820725/93‑1 (MQ=255)
                                                                        cATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCATTGCGAt                 <  2:4608561/93‑1 (MQ=255)
                                                                        cATTGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGAt                 <  1:1722000/93‑1 (MQ=255)
                                                                          ttGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGAtca               >  1:3826362/1‑91 (MQ=255)
                                                                          ttGAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATaa               >  2:2439251/1‑93 (MQ=255)
                                                                            gAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAg             <  2:4248128/93‑1 (MQ=255)
                                                                            gAACATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAg             <  1:3190486/93‑1 (MQ=255)
                                                                             gaCATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGa            <  1:2819285/92‑1 (MQ=255)
                                                                              aCATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGAc           <  2:2302974/93‑1 (MQ=255)
                                                                              aCATTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGAc           <  2:2358750/93‑1 (MQ=255)
                                                                                aTTTGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCg         >  1:2807176/1‑93 (MQ=255)
                                                                                 tttGCAGCACCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGa        >  1:1741297/1‑93 (MQ=255)
                                                                                    gcagcaCCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGc     >  2:3755701/1‑93 (MQ=255)
                                                                                      cgcaCCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCa   <  2:3808787/92‑1 (MQ=255)
                                                                                          |                                                                                         
CTTCGCTGATCGGCATGTAAGGCAGCTGCTTATCTGCTCCGACATCCAGCGTACGCAAGGTGACGGGTTTATCATTGAACATTTGCAGCATCCCCTGATACTGCGCCACCTGTTCTTCTTCCGACGGGAAACCACTTTGCAGCATGAATGGGATTTCAGTGCGATAAAGACCGATGCCAT  >  NC_000913/2964911‑2965090

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