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.7% 1097.4 / NA 350H263R (CAT→CGT) metEtetrahydropteroyltriglutamate methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (167/182);  total (167/183)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 366 total aligned reads displayed.

AGGGACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCATGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATG  >  NC_000913/4011340‑4011519
                                                                                          |                                                                                         
aGGTACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTg                                                                                         >  2:2501737/1‑93 (MQ=255)
 gggACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCgtgt                                                                                        >  2:7263676/1‑93 (MQ=255)
 gggACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCgtgt                                                                                        <  2:1091717/93‑1 (MQ=255)
 gggACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCgtgt                                                                                        <  1:7546623/93‑1 (MQ=255)
  ggACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGtt                                                                                       >  2:6069416/1‑93 (MQ=255)
  ggACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGtt                                                                                       <  2:3986828/93‑1 (MQ=255)
   gACAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTg                                                                                      <  1:7679254/93‑1 (MQ=255)
    aCAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGa                                                                                     <  1:4709151/93‑1 (MQ=255)
    aCAGGTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGa                                                                                     >  2:2394236/1‑93 (MQ=255)
    aCAGGTGAAACTGCTGCTGACCACCTATTTCGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGa                                                                                     <  2:3502444/93‑1 (MQ=255)
       ggTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCt                                                                                  >  1:1413298/1‑93 (MQ=255)
       ggTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCt                                                                                  >  2:3365971/1‑93 (MQ=255)
        gTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTc                                                                                 >  2:589966/1‑93 (MQ=255)
        gTGAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTc                                                                                 >  1:64214/1‑93 (MQ=255)
          gAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGt                                                                               >  1:398252/1‑93 (MQ=255)
          gAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGt                                                                               <  1:114281/93‑1 (MQ=255)
          gAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGt                                                                               <  2:3306440/93‑1 (MQ=255)
          gAAACTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGt                                                                               <  2:5533706/93‑1 (MQ=255)
             aCTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACa                                                                            <  1:3684155/93‑1 (MQ=255)
             aCTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACa                                                                            <  2:2028866/93‑1 (MQ=255)
             aCTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACa                                                                            <  1:3029928/93‑1 (MQ=255)
             aCTGCTGCTGACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGGCTGCGTGTTGACCTCGTACa                                                                            >  1:5963205/1‑93 (MQ=255)
              ctgctgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGAc                                                                                    >  1:333428/1‑84 (MQ=255)
              ctgctgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGAc                                                                                    <  2:333429/84‑1 (MQ=255)
               tgctgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATg                                                                          >  2:6150769/1‑93 (MQ=255)
                  tgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCATGGTCTGCGTGTTGACCTCGTACATGGTa                                                                       >  2:5738125/1‑93 (MQ=255)
                  tgctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTa                                                                       <  1:5419104/93‑1 (MQ=255)
                   gctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTaa                                                                      <  2:378261/93‑1 (MQ=255)
                   gctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTaa                                                                      >  2:739025/1‑93 (MQ=255)
                    ctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTaaa                                                                     <  2:2221904/93‑1 (MQ=255)
                    ctgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTaaa                                                                     <  1:351531/93‑1 (MQ=255)
                     tgACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAg                                                                    >  2:4434764/1‑93 (MQ=255)
                      gACCACCTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGa                                                                   >  1:1271755/1‑93 (MQ=255)
                       accaccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGAt                                                                  <  2:3055386/93‑1 (MQ=255)
                        ccaccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATg                                                                 <  1:7304687/93‑1 (MQ=255)
                        ccaccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATg                                                                 >  1:4999374/1‑93 (MQ=255)
                        ccaccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGACTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATg                                                                 >  2:2690873/1‑93 (MQ=255)
                         caccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGa                                                                <  2:288964/93‑1 (MQ=255)
                         caccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGa                                                                <  1:5227821/93‑1 (MQ=255)
                          accTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGAc                                                               <  1:6589700/93‑1 (MQ=255)
                          accTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGAc                                                               >  2:794438/1‑93 (MQ=255)
                           ccTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACg                                                              >  2:4995247/1‑93 (MQ=255)
                            cTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGt                                                             <  2:834615/93‑1 (MQ=255)
                            cTATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGt                                                             >  1:2804111/1‑93 (MQ=255)
                             tATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACgct                                                            <  1:4656337/93‑3 (MQ=255)
                             tATTTTGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGtt                                                            >  1:1636910/1‑93 (MQ=255)
                               ttttGAAGGCGTAACGCCAAATTTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGc                                                          >  1:5671304/1‑93 (MQ=255)
                               ttttGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGc                                                          >  1:5290001/1‑93 (MQ=255)
                               ttttGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGc                                                          <  2:5250029/93‑1 (MQ=255)
                               ttttGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGc                                                          >  1:7368499/1‑93 (MQ=255)
                               ttttGAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGc                                                          >  2:7619867/1‑93 (MQ=255)
                                   gAAGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGaa                                                      >  1:5518304/1‑93 (MQ=255)
                                    aaGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAAc                                                     >  2:6903667/1‑93 (MQ=255)
                                     aGGCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACt                                                    >  1:7220064/1‑93 (MQ=255)
                                      ggCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTg                                                   <  2:3176570/93‑1 (MQ=255)
                                      ggCGTAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTg                                                   <  2:5642168/93‑1 (MQ=255)
                                          tAACGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGcaca                                               >  1:7931937/1‑93 (MQ=255)
                                           aaCGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACaa                                              <  2:3892851/93‑1 (MQ=255)
                                           aaCGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACaa                                              >  1:4898605/1‑93 (MQ=255)
                                            aCGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAg                                             <  2:939154/93‑1 (MQ=255)
                                             cGCCAAATCTCGACACGATTACTGCGCTGCCTGTTCATGGTCTGCGTGTTGCCCGAGTACATGGTAAAGATGACGTTGCTGAACTGCACAAgc                                            <  1:2709684/93‑1 (MQ=255)
                                              gCCAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAgcg                                           >  1:3612130/1‑93 (MQ=255)
                                                cAAATCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcc                                         <  2:6351800/93‑1 (MQ=255)
                                                 aaaTCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcct                                        <  1:3388207/93‑1 (MQ=255)
                                                 aaaTCTCGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcct                                        >  1:4215167/1‑93 (MQ=255)
                                                    tctcGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcctgcc                                     <  2:6515124/93‑1 (MQ=255)
                                                     ctcGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCgcctgcct                                    <  2:7213612/93‑1 (MQ=255)
                                                      tcGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCtt                                   <  1:1569352/93‑1 (MQ=255)
                                                       cGACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTc                                  <  2:2483836/93‑1 (MQ=255)
                                                        gACACGATTACTGCGCTGCCTGTTCAGGGTCTGTGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCt                                 >  1:2994979/1‑93 (MQ=255)
                                                        gACACGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCt                                 <  1:4093667/93‑1 (MQ=255)
                                                          cacGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGa                               >  1:4139303/1‑93 (MQ=255)
                                                          cacGATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGa                               <  1:4471590/93‑1 (MQ=255)
                                                           ccGAGTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGAc                              <  1:6040872/92‑1 (MQ=255)
                                                             gATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTg                            >  1:705469/1‑93 (MQ=255)
                                                             gATTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTg                            >  1:1028349/1‑93 (MQ=255)
                                                             gATTACTGCGCTGCCTGTTCAGGGTATGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTg                            >  2:6708497/1‑93 (MQ=255)
                                                              aTTACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTgg                           >  1:3105237/1‑93 (MQ=255)
                                                               ttACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGt                          >  2:4325994/1‑93 (MQ=255)
                                                               ttACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGt                          <  1:7756078/93‑1 (MQ=255)
                                                                tACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGtt                         <  1:182418/93‑1 (MQ=255)
                                                                tACTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGtt                         >  1:679692/1‑93 (MQ=255)
                                                                 aCTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTg                        >  1:2425249/1‑93 (MQ=255)
                                                                  cTGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGc                       <  1:5894501/93‑1 (MQ=255)
                                                                   tGCGCTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCt                      <  2:7858393/93‑1 (MQ=255)
                                                                    gcgcTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTg                     <  2:7348023/93‑1 (MQ=255)
                                                                      gcTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTc                   <  2:6632633/93‑1 (MQ=255)
                                                                       cTGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACCGGCTGCTGTCt                  >  2:7007866/1‑93 (MQ=255)
                                                                        tGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTg                 >  1:2127701/1‑93 (MQ=255)
                                                                        tGCCTGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTg                 >  2:4123407/1‑93 (MQ=255)
                                                                         gCCTGTTCAGGGTCTGCGTGTTGACCTGGTACATGGTAAAGATGACGTTGCTGAACTGCCCAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGc                <  2:5353114/93‑1 (MQ=255)
                                                                            tGTTCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCggg             <  2:6454317/93‑1 (MQ=255)
                                                                              ttCCGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGCACTGCACACGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTc           <  1:2555030/93‑1 (MQ=255)
                                                                               tCAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCt          <  1:6893891/93‑1 (MQ=255)
                                                                                cAGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTg         >  1:648345/1‑93 (MQ=255)
                                                                                 aGGGTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGa        <  1:1487469/93‑1 (MQ=255)
                                                                                  gggTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGAt       <  2:1818649/93‑1 (MQ=255)
                                                                                    gTCTGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGCCTGATCa     >  1:5074984/1‑93 (MQ=255)
                                                                                       tGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCATAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATg  <  2:4675807/93‑1 (MQ=255)
                                                                                       tGCGTGTTGACCTCGTACATGGTAAAGATGACGTTGCTGAACTGCACAAGCGCCTGCCTTCTGACTGGTTGCTGTCTGCGGGTCTGATCAATg  <  1:2093698/93‑1 (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: