Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,011,430 A→G H263R (CAT→CGT)  metE → tetrahydropteroyltriglutamate methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,011,4300AG99.6% 779.7 / ‑5.4 242H263R (CAT→CGT) metEtetrahydropteroyltriglutamate methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (122/119);  minor base C (0/1);  total (122/120)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 254 total aligned reads displayed.

AGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATG  >  NC_000913/4011340‑4011519
                                                                                          |                                                                                         
 gggACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCgtgt                                                                                        >  2:2854158/1‑93 (MQ=255)
 gggACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCgtgt                                                                                        >  1:6185819/1‑93 (MQ=255)
 gggACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCgtgt                                                                                        >  2:4247962/1‑93 (MQ=255)
 gggACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCgtgt                                                                                        >  2:6041328/1‑93 (MQ=255)
  ggACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGtt                                                                                       >  1:502707/1‑93 (MQ=255)
  ggACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGtt                                                                                       >  1:413384/1‑93 (MQ=255)
   gACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTg                                                                                      >  2:2651232/1‑93 (MQ=255)
    aCAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGa                                                                                     >  2:402541/1‑93 (MQ=255)
    aCAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGa                                                                                     <  1:1588900/93‑1 (MQ=255)
    aCAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGa                                                                                     >  2:3836082/1‑93 (MQ=255)
    aCAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGa                                                                                     <  2:1493645/93‑1 (MQ=255)
      aGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGAcc                                                                                   >  1:5517188/1‑93 (MQ=255)
         tGAAGCTGCTGCTGACCACCTATTGTGAAGGCGTAACGCCAAATCGCGACGCGATGACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCg                                                                                <  1:408519/93‑1 (MQ=255)
           aacctgctgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTCCTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTa                                                                              <  2:5437897/90‑1 (MQ=255)
           aaaCTGCTTCTGACTACCTATTTTGAAGTCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGTCCTCGTa                                                                              <  1:2248455/93‑1 (MQ=255)
           aaaCTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTa                                                                              >  2:2228603/1‑93 (MQ=255)
           aaaCTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTa                                                                              <  1:2768931/93‑1 (MQ=255)
              ctgctgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACAt                                                                           <  2:3059315/93‑1 (MQ=255)
              ctgctgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACAt                                                                           >  2:1227130/1‑93 (MQ=255)
              ctgctgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACAt                                                                           >  1:1508228/1‑93 (MQ=255)
              ctgctgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACAt                                                                           >  2:4418690/1‑93 (MQ=255)
                gctgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATgg                                                                         <  2:2689327/93‑1 (MQ=255)
                   gctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCTTACATGGTaa                                                                      <  2:4174089/93‑1 (MQ=255)
                   gctgACACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTaaa                                                                     <  2:6236274/93‑1 (MQ=255)
                    ctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTaaa                                                                     <  1:3413027/93‑1 (MQ=255)
                     tgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAg                                                                    >  2:3984966/1‑93 (MQ=255)
                      gACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGa                                                                   >  2:4831797/1‑93 (MQ=255)
                       accaccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGAt                                                                  <  1:5813727/93‑1 (MQ=255)
                        tcaccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATTGTAAAGATg                                                                 <  2:1046734/92‑1 (MQ=255)
                        ccaccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATg                                                                 <  2:3135614/93‑1 (MQ=255)
                         caccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGa                                                                <  2:4546153/93‑1 (MQ=255)
                          accTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGAc                                                               <  1:4973540/93‑1 (MQ=255)
                          accTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGAc                                                               >  2:1859022/1‑93 (MQ=255)
                           ccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACg                                                              <  2:3759917/93‑1 (MQ=255)
                           ccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACg                                                              <  2:4683798/93‑1 (MQ=255)
                            cTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGt                                                             >  1:4651227/1‑93 (MQ=255)
                             tATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGtt                                                            <  2:2202818/93‑1 (MQ=255)
                             tATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGtt                                                            <  2:130119/93‑1 (MQ=255)
                             tATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGtt                                                            >  2:5667367/1‑93 (MQ=255)
                                   gAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGaa                                                      >  1:5475165/1‑93 (MQ=255)
                                     aGGCGTATCGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACt                                                    <  2:1528568/93‑1 (MQ=255)
                                        cGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCCTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGca                                                 <  1:5846371/93‑1 (MQ=255)
                                         gTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGcac                                                >  1:3791880/1‑93 (MQ=255)
                                         gTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGcac                                                >  2:5567018/1‑93 (MQ=255)
                                            aCGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAg                                             >  1:3987355/1‑93 (MQ=255)
                                            aCGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAg                                             <  1:2907394/93‑1 (MQ=255)
                                             cGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAgc                                            >  2:5890929/1‑93 (MQ=255)
                                              gCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAgcg                                           >  2:2874519/1‑93 (MQ=255)
                                                cAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcc                                         >  2:2478993/1‑93 (MQ=255)
                                                cAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcc                                         >  2:779861/1‑93 (MQ=255)
                                                 aaaTCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcct                                        <  1:5328553/93‑1 (MQ=255)
                                                   aTCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcctgc                                      >  2:358322/1‑93 (MQ=255)
                                                   aTCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcctgc                                      >  1:4714170/1‑93 (MQ=255)
                                                     ctcGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcctgcct                                    <  2:919972/93‑1 (MQ=255)
                                                      tcGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCtt                                   >  1:3752767/1‑93 (MQ=255)
                                                        gACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCt                                 >  2:5066088/1‑93 (MQ=255)
                                                         acacGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTg                                >  1:418916/1‑93 (MQ=255)
                                                         acacGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTg                                >  1:4492735/1‑93 (MQ=255)
                                                          cacGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGa                               <  1:5774012/93‑1 (MQ=255)
                                                          cacGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGa                               >  1:3241867/1‑93 (MQ=255)
                                                          cacGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGa                               >  2:5397842/1‑93 (MQ=255)
                                                           acGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGAc                              <  2:4757287/93‑1 (MQ=255)
                                                           acGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGAc                              <  2:3914647/93‑1 (MQ=255)
                                                            cGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGAca                             <  1:899774/93‑2 (MQ=255)
                                                            cGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACt                             <  1:773641/93‑1 (MQ=255)
                                                            cGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACt                             <  1:1715249/93‑1 (MQ=255)
                                                            cGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACt                             >  1:4294217/1‑93 (MQ=255)
                                                            cGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACt                             >  1:2332896/1‑93 (MQ=255)
                                                             gATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTg                            >  2:1177017/1‑93 (MQ=255)
                                                               ttACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGt                          <  2:2996619/93‑1 (MQ=255)
                                                                  cTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGc                       <  1:5194651/93‑1 (MQ=255)
                                                                   tGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCt                      >  1:2884048/1‑93 (MQ=255)
                                                                   tGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCt                      <  2:2451552/93‑1 (MQ=255)
                                                                      gcTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTc                   <  2:2082543/93‑1 (MQ=255)
                                                                       cTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCt                  >  1:2304632/1‑93 (MQ=255)
                                                                       cTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCt                  >  1:62957/1‑93 (MQ=255)
                                                                        tGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTg                 <  2:5496179/93‑1 (MQ=255)
                                                                         gCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGc                >  2:2283480/1‑93 (MQ=255)
                                                                         gCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGc                <  1:5924840/93‑1 (MQ=255)
                                                                         gCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGc                >  2:2139064/1‑93 (MQ=255)
                                                                           cTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCgg              <  1:1168715/93‑1 (MQ=255)
                                                                           cTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCgg              >  1:1049326/1‑93 (MQ=255)
                                                                           cTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCgg              >  1:3529478/1‑93 (MQ=255)
                                                                            tGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCggg             >  2:3158653/1‑93 (MQ=255)
                                                                            tGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCggg             >  1:2540100/1‑93 (MQ=255)
                                                                              ttCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTc           <  2:951974/93‑1 (MQ=255)
                                                                              ttCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTc           <  1:5605348/93‑1 (MQ=255)
                                                                               tCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCt          <  2:3695455/93‑1 (MQ=255)
                                                                                cAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTg         <  2:6252215/93‑1 (MQ=255)
                                                                                  gggTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGAt       <  2:5248913/93‑1 (MQ=255)
                                                                                  gggTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGAt       <  1:858802/93‑1 (MQ=255)
                                                                                  gggTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGAt       >  1:2415614/1‑93 (MQ=255)
                                                                                     tCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCaa    <  2:5353016/93‑1 (MQ=255)
                                                                                     tCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCaa    >  1:3107040/1‑93 (MQ=255)
                                                                                     tCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCaa    >  1:2173957/1‑93 (MQ=255)
                                                                                     tCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCaa    >  1:1687903/1‑93 (MQ=255)
                                                                                      cTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAAt   <  1:2030292/93‑1 (MQ=255)
                                                                                      cTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAAt   <  2:1919781/93‑1 (MQ=255)
                                                                                      cTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCATTCTGACTGGTTGCTGTCTGCGGGTCTGATCAAt   <  1:2093972/93‑1 (MQ=255)
                                                                                       tGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATg  >  1:5072032/1‑93 (MQ=255)
                                                                                          |                                                                                         
AGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATG  >  NC_000913/4011340‑4011519

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