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.0% 1210.4 / ‑2.1 380H263R (CAT→CGT) metEtetrahydropteroyltriglutamate methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/3);  new base G (192/184);  total (193/187)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 399 total aligned reads displayed.

AGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATG  >  NC_000913/4011340‑4011519
                                                                                          |                                                                                         
aGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTg                                                                                         >  1:7193149/1‑93 (MQ=255)
 gggACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCgtgt                                                                                        <  2:9833923/93‑1 (MQ=255)
 gggACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCgtgt                                                                                        >  2:7375911/1‑93 (MQ=255)
 gggACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCgtgt                                                                                        >  2:5020742/1‑93 (MQ=255)
 gggACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCgtgt                                                                                        >  1:2406744/1‑93 (MQ=255)
  ggACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGtt                                                                                       <  1:2891486/93‑1 (MQ=255)
     cAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGAc                                                                                    <  2:3164244/93‑1 (MQ=255)
     cAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGAc                                                                                    >  1:8788763/1‑93 (MQ=255)
     cAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGAc                                                                                    <  1:3332971/93‑1 (MQ=255)
      aGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGAcc                                                                                   >  1:5707377/1‑93 (MQ=255)
       ggTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCt                                                                                  <  2:4732581/93‑1 (MQ=255)
       ggTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCt                                                                                  >  2:4115791/1‑93 (MQ=255)
       ggTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCt                                                                                  <  2:3891835/93‑1 (MQ=255)
        gTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTc                                                                                 <  1:9171575/93‑1 (MQ=255)
        gTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTc                                                                                 <  1:2859407/93‑1 (MQ=255)
         tGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCg                                                                                >  1:1628822/1‑93 (MQ=255)
         tGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCg                                                                                >  2:7034248/1‑93 (MQ=255)
          gAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGt                                                                               >  1:208659/1‑93 (MQ=255)
           aaaCTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGACTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTa                                                                              <  2:9195245/93‑1 (MQ=255)
            aaCTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTAc                                                                             >  2:3095597/1‑93 (MQ=255)
            aaCTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTAc                                                                             >  2:967129/1‑93 (MQ=255)
              ctggtgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACAt                                                                           <  1:6396194/93‑1 (MQ=255)
              ctgctgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACAt                                                                           <  2:7186370/93‑1 (MQ=255)
              ctgctgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACAt                                                                           <  1:4191101/93‑1 (MQ=255)
                 ctgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGt                                                                        >  1:1533491/1‑93 (MQ=255)
                  tgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTa                                                                       <  1:6201191/93‑1 (MQ=255)
                   gctgAACAATTATTTTGAAGGCGTAACGCCAAATCACGACGCGATTACTGTGCTTCCTGTTCAGGGTCTGCGTGTTGACTTCGTACATGGTaa                                                                      <  2:6305603/93‑1 (MQ=255)
                     tgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAg                                                                    >  1:7107102/1‑93 (MQ=255)
                     tgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAg                                                                    >  2:7663593/1‑93 (MQ=255)
                     tgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAg                                                                    >  2:897494/1‑93 (MQ=255)
                      gACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGa                                                                   >  2:10044986/1‑93 (MQ=255)
                      gACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGa                                                                   >  2:357280/1‑93 (MQ=255)
                      gACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGa                                                                   >  1:3770367/1‑93 (MQ=255)
                        ccaccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATg                                                                 >  1:4934119/1‑93 (MQ=255)
                        ccaccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATg                                                                 >  1:2226605/1‑93 (MQ=255)
                        ccaccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATg                                                                 <  2:1234337/93‑1 (MQ=255)
                           ccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACg                                                              <  2:6637272/93‑1 (MQ=255)
                            cTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGt                                                             >  1:6121992/1‑93 (MQ=255)
                            cTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGt                                                             <  1:3969542/93‑1 (MQ=255)
                            cTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGt                                                             >  1:7562440/1‑93 (MQ=255)
                             tATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGtt                                                            <  1:2196894/93‑1 (MQ=255)
                              aTTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTg                                                           >  1:3499889/1‑93 (MQ=255)
                               ttttGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGc                                                          <  2:4683961/93‑1 (MQ=255)
                               ttttGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGc                                                          <  1:5375475/93‑1 (MQ=255)
                                 ttGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTg                                                        <  2:8816628/93‑1 (MQ=255)
                                  tGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGa                                                       >  2:3517767/1‑93 (MQ=255)
                                   gAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGaa                                                      <  2:2876984/93‑1 (MQ=255)
                                    aaGGCGTAACTCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTTACCTCGTACATGGTAAAGATGACGTTGCTGAAc                                                     <  1:1364806/93‑1 (MQ=255)
                                     aGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACt                                                    <  1:4141095/93‑1 (MQ=255)
                                      ggCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTg                                                   >  1:6988947/1‑93 (MQ=255)
                                        cGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGca                                                 <  1:251153/93‑1 (MQ=255)
                                          tAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGcaca                                               >  1:9401634/1‑93 (MQ=255)
                                           aaCGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACaa                                              <  2:9552962/93‑1 (MQ=255)
                                           aaCGCCAAATCTCGACACGATCACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACaa                                              <  2:331284/93‑1 (MQ=255)
                                            aCGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAg                                             <  1:47475/93‑1 (MQ=255)
                                             cGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAgc                                            <  1:129432/93‑1 (MQ=255)
                                             cGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAgc                                            <  1:9428430/93‑1 (MQ=255)
                                             cGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAgc                                            <  2:652580/93‑1 (MQ=255)
                                             cGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAgc                                            >  2:8465041/1‑93 (MQ=255)
                                              gCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAgcg                                           >  1:9711058/1‑93 (MQ=255)
                                               ccAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAgcgc                                          <  2:5569056/93‑1 (MQ=255)
                                               ccAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAgcgc                                          <  1:6492769/93‑1 (MQ=255)
                                                cAAATCTCGACACGAGTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcc                                         <  1:2755318/93‑1 (MQ=255)
                                                 aaaTCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcct                                        >  2:1820664/1‑93 (MQ=255)
                                                 aaaTCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcct                                        >  1:121102/1‑93 (MQ=255)
                                                 aaaTCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcct                                        >  1:7498590/1‑93 (MQ=255)
                                                 aaaTCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcct                                        <  1:7913707/93‑1 (MQ=255)
                                                   aTCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcctgc                                      <  2:6109333/93‑1 (MQ=255)
                                                    tctcGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcctgcc                                     >  2:1998161/1‑93 (MQ=255)
                                                     ctcGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcctgcct                                    <  2:7269294/93‑1 (MQ=255)
                                                      tcGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCtt                                   <  1:2584317/93‑1 (MQ=255)
                                                       cGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTc                                  <  1:5220233/93‑1 (MQ=255)
                                                       cGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGAACTCGTACATGGTAAAGATGAAGTTTCTGAACTCCCCAAGAGCCCGCAttt                                  >  1:8112533/1‑92 (MQ=255)
                                                          cacGATTACTGCGCTGCCTGTTCCGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGa                               <  2:2307497/93‑1 (MQ=255)
                                                          cacGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGa                               <  1:3675246/93‑1 (MQ=255)
                                                            cGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACt                             >  2:597365/1‑93 (MQ=255)
                                                              aTTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTgg                           <  1:6834429/93‑1 (MQ=255)
                                                               ttACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGt                          >  2:1056295/1‑93 (MQ=255)
                                                                tACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGtt                         >  2:4350441/1‑93 (MQ=255)
                                                                tACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGtt                         <  2:7846309/93‑1 (MQ=255)
                                                                tACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGtt                         >  2:6755891/1‑93 (MQ=255)
                                                                tACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGtt                         >  2:4515858/1‑93 (MQ=255)
                                                                    gcgcTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTg                     >  1:8383010/1‑93 (MQ=255)
                                                                     cgcTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGt                    <  1:5506137/93‑1 (MQ=255)
                                                                      gcTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTc                   <  2:4606303/93‑1 (MQ=255)
                                                                       cTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCt                  <  1:171100/93‑1 (MQ=255)
                                                                        tGCCTGTTCTGGGTCTGCGTGTTGCCCTCGTACATGGTAAAGCTGCCGTTGCTGACCTGCACCAGCGCCTGCCTTCTGACTGGTTGCTGTCTg                 <  1:7693154/93‑1 (MQ=255)
                                                                        tGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTg                 >  1:911775/1‑93 (MQ=255)
                                                                         gCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGc                <  2:8644628/93‑1 (MQ=255)
                                                                            tGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCggg             >  2:5242487/1‑93 (MQ=255)
                                                                               tCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCt          >  1:1032341/1‑93 (MQ=255)
                                                                               tCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCt          <  1:4453056/93‑1 (MQ=255)
                                                                               tCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCt          <  1:2273160/93‑1 (MQ=255)
                                                                               tCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCt          >  2:7993269/1‑93 (MQ=255)
                                                                                 aGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGa        >  2:8351318/1‑93 (MQ=255)
                                                                                 aGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGa        <  2:5452167/93‑1 (MQ=255)
                                                                                  gggTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGAt       <  2:6057314/93‑1 (MQ=255)
                                                                                    gTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCa     <  1:7327080/93‑1 (MQ=255)
                                                                                    gTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCa     >  1:3380968/1‑93 (MQ=255)
                                                                                     tCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCaa    >  1:8664799/1‑93 (MQ=255)
                                                                                          |                                                                                         
AGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATG  >  NC_000913/4011340‑4011519

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