Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,251,782 A→G V755A (GTG→GCG)  plsB ← glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,251,7820AG99.8% 1317.5 / NA 413V755A (GTG→GCG) plsBglycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (187/224);  minor base C (1/0);  total (189/224)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 435 total aligned reads displayed.

TGATAAACCTTCATCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCG  >  NC_000913/4251692‑4251871
                                                                                          |                                                                                         
tGATAAACCTTCATCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGcc                                                                                         <  2:5798635/93‑1 (MQ=255)
 gATAAACCTTCATCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCg                                                                                        <  1:478076/93‑1 (MQ=255)
 gATAAACCTTCATCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCg                                                                                        >  2:4612021/1‑93 (MQ=255)
 gATAAACCTTCATCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCg                                                                                        >  2:940737/1‑93 (MQ=255)
  aTAAACCTTCATCGTTTCTGCCGGTTCGACATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAAcgccgc                                                                                       <  1:9881893/93‑1 (MQ=255)
    aaaCCTTCATCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCt                                                                                     >  1:256756/1‑93 (MQ=255)
     acccTTCACCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCtt                                                                                    <  1:4408939/91‑1 (MQ=255)
        cTTCATCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTc                                                                                 <  1:6314998/93‑1 (MQ=255)
        cTTCATCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTc                                                                                 <  2:6514832/93‑1 (MQ=255)
          tCATCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCga                                                                               <  1:8529479/93‑1 (MQ=255)
           cATCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCgaa                                                                              >  2:8330458/1‑93 (MQ=255)
            aTCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCgaag                                                                             >  1:7829593/1‑93 (MQ=255)
             tCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCgaaga                                                                            >  2:3654818/1‑93 (MQ=255)
             tCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCgaaga                                                                            <  1:1852679/93‑1 (MQ=255)
             tCGTGTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCgaaga                                                                            <  1:6582750/93‑1 (MQ=255)
             tCGGTTCTGTCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCgaaga                                                                            <  2:3633995/93‑1 (MQ=255)
              cGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCgaagaa                                                                           <  1:3067749/93‑1 (MQ=255)
              cGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCgaagaa                                                                           <  1:968941/93‑1 (MQ=255)
              cGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCgaagaa                                                                           <  1:6046749/93‑1 (MQ=255)
                tttCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACt                                                                         >  1:2901173/1‑93 (MQ=255)
                 ttCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTc                                                                        <  2:8743960/93‑1 (MQ=255)
                  tCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTcc                                                                       <  1:5833837/93‑1 (MQ=255)
                    tGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCgg                                                                     >  2:2434209/1‑93 (MQ=255)
                    tGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCgg                                                                     <  2:5000266/93‑1 (MQ=255)
                    tGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCgg                                                                     <  1:7166401/93‑1 (MQ=255)
                    tGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCgg                                                                     >  2:9059947/1‑93 (MQ=255)
                      ccGGTGCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGgcg                                                                   <  2:1576224/93‑1 (MQ=255)
                         gTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGt                                                                <  1:8080138/93‑1 (MQ=255)
                            cGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGGAGAACTCCGGCGCGTTGa                                                             <  1:3342961/93‑1 (MQ=255)
                            cGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGa                                                             <  1:6919799/93‑1 (MQ=255)
                            cGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGa                                                             <  2:5665120/93‑1 (MQ=255)
                             ggCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGAt                                                            <  2:7874916/93‑1 (MQ=255)
                                aTCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGcc                                                         >  1:10006327/1‑93 (MQ=255)
                                 tCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCg                                                        >  2:6847904/1‑93 (MQ=255)
                                 tCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCg                                                        <  2:3716834/93‑1 (MQ=255)
                                 acgccgcTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCg                                                        <  2:7012209/92‑1 (MQ=255)
                                  cgccgcTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGt                                                       <  2:9297497/93‑1 (MQ=255)
                                    ccgcTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCGGCCTTGTCGAAGAACTCCGGCGCGGTGATGCCGTgc                                                     >  1:402145/1‑93 (MQ=255)
                                     cgcTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTgca                                                    >  1:1143936/1‑93 (MQ=255)
                                      gcTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTgcaa                                                   >  1:9305543/1‑92 (MQ=255)
                                       cTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTgcagc                                                  >  1:1570048/1‑93 (MQ=255)
                                       cTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTgcagc                                                  <  2:2820622/93‑1 (MQ=255)
                                       cTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTgcagc                                                  >  2:621503/1‑93 (MQ=255)
                                        tATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTgcagca                                                 <  1:7414106/93‑1 (MQ=255)
                                         aTCGCTGATATTCCCTTCATCACGCAGTGTCAGCAACAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCAc                                                <  2:4881123/93‑1 (MQ=255)
                                          tCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACg                                               <  2:7695209/93‑1 (MQ=255)
                                          tCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACg                                               >  2:7386383/1‑93 (MQ=255)
                                            gCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGga                                             <  2:2143385/93‑1 (MQ=255)
                                             cTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGgag                                            <  1:6993020/93‑1 (MQ=255)
                                             cTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGgag                                            >  2:192780/1‑93 (MQ=255)
                                             cTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGgag                                            <  1:8782479/93‑1 (MQ=255)
                                              tGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGgaga                                           <  1:9087873/93‑1 (MQ=255)
                                              tGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGgaga                                           >  1:3773098/1‑93 (MQ=255)
                                              tGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGgaga                                           <  1:2122702/93‑1 (MQ=255)
                                               gATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGgagag                                          <  2:5267932/93‑1 (MQ=255)
                                               gATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGgagag                                          >  1:5177195/1‑93 (MQ=255)
                                               gATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGgagag                                          <  2:3881400/93‑1 (MQ=255)
                                               gATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGgagag                                          >  1:2073343/1‑93 (MQ=255)
                                                atatACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGgagaga                                         <  1:2287636/93‑1 (MQ=255)
                                                   tgcccTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGt                                      <  2:10007102/91‑1 (MQ=255)
                                                   taCCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGt                                      <  2:2723726/93‑1 (MQ=255)
                                                    aCCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGtt                                     <  1:6406491/93‑1 (MQ=255)
                                                    aCCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGtt                                     <  1:7240693/93‑1 (MQ=255)
                                                    aCCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGtt                                     <  2:558132/93‑1 (MQ=255)
                                                    aCCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGtt                                     >  1:1313649/1‑93 (MQ=255)
                                                     cccTTCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTg                                    <  2:3233652/93‑1 (MQ=255)
                                                        ttCATCACGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTgcgc                                 >  1:7342743/1‑93 (MQ=255)
                                                         tcatcaCGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGcgcg                                >  2:6350934/1‑93 (MQ=255)
                                                         ccatcaCGCAGTGTCAGCATCAGAGAACTGAATGCCGCCTTGTCGAAGAACTCCGGCGCTTTGAAGCCTTTCATCTCGGAGAGACGTTGcgcg                                >  1:8974678/2‑93 (MQ=255)
                                                           atcaCGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGag                              >  1:7554767/1‑92 (MQ=255)
                                                           atcaCGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGAc                              >  2:2564128/1‑93 (MQ=255)
                                                            tcaCGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACg                             >  2:9702564/1‑93 (MQ=255)
                                                             caCGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACgg                            <  2:9761146/93‑1 (MQ=255)
                                                               cGCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTg                          <  2:4013446/93‑1 (MQ=255)
                                                                gCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGACGAACTCCGGCGCGTTCATGCCGTGCCGCAACGAGAGACGTTGCTCGACgggcg                         >  2:1415936/1‑90 (MQ=255)
                                                                gCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGc                         <  1:1351692/93‑1 (MQ=255)
                                                                gCAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGc                         <  2:6617334/93‑1 (MQ=255)
                                                                 cAGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCg                        <  1:1687875/93‑1 (MQ=255)
                                                                  aGTGTCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCgg                       >  2:1025732/1‑93 (MQ=255)
                                                                   gtgtCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGc                      <  2:8915226/93‑1 (MQ=255)
                                                                    tgtCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCTCGACGGAGCGGcc                     >  2:4261128/1‑92 (MQ=255)
                                                                    tgtCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGct                     >  1:8320751/1‑93 (MQ=255)
                                                                    tgtCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGcg                     >  2:5484421/1‑92 (MQ=255)
                                                                      tCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGctct                   >  2:4431796/1‑93 (MQ=255)
                                                                      tCAGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGctct                   >  2:2030247/1‑93 (MQ=255)
                                                                        aGCACCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGctctct                 <  1:5744816/93‑1 (MQ=255)
                                                                           aCCAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTtc              <  1:498811/93‑1 (MQ=255)
                                                                            ccAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTtct             >  1:3232331/1‑93 (MQ=255)
                                                                             cAGAGAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTtctc            >  2:7247847/1‑93 (MQ=255)
                                                                              agagAACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTcc           <  2:8202427/93‑1 (MQ=255)
                                                                               gagaACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCa          <  1:4125090/93‑1 (MQ=255)
                                                                               gagaACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCa          >  1:3641951/1‑93 (MQ=255)
                                                                                 gaACTGACCGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGc        <  2:51068/93‑1 (MQ=255)
                                                                                 gaACTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGCGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGc        >  1:3562501/1‑93 (MQ=255)
                                                                                    cTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTa     <  2:1162028/93‑1 (MQ=255)
                                                                                    cTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTa     >  1:5269798/1‑93 (MQ=255)
                                                                                    cTGAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTa     <  2:9643032/93‑1 (MQ=255)
                                                                                     tGAACGCCGCCTTGTCGAAGACCTCCGGCGCGTTGATGCCGTGCAGCACGGATAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTAc    <  2:6000597/93‑1 (MQ=255)
                                                                                     tGAACGCCGCCTTGTCGAAGAACTCCGGCGCGCTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGATCCCTTTCTCCACCGTAc    >  2:2279476/1‑93 (MQ=255)
                                                                                      gAACGCCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTAcc   <  1:4543206/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGATAAACCTTCATCGTTTCTGCCGGTTCGGCATCGCCGCTATCGCTGATATACCCTTCATCACGCAGTGTCAGCACCAGAGAACTGAACACCGCCTTGTCGAAGAACTCCGGCGCGTTGATGCCGTGCAGCACGGAGAGACGTTGCGCGACGGTGCGGCTCTCTTTCTCCAGCGTACCG  >  NC_000913/4251692‑4251871

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