Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 467,839 T→C 21.3% V143A (GTC→GCC)  cof → thiamine pyrimidine pyrophosphate hydrolase; HMP‑PP phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913467,8390TC21.3% 341.0 / 95.3 183V143A (GTC→GCC) cofthiamine pyrimidine pyrophosphate hydrolase; HMP‑PP phosphatase
Reads supporting (aligned to +/- strand):  ref base T (77/67);  new base C (21/18);  total (98/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

Only 100 of 183 total aligned reads displayed.

GTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTTGTGTTTTTCC  >  NC_000913/467741‑467936
                                                                                                  |                                                                                                 
gttTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTc                                                                                                  >  1:1639233/1‑100 (MQ=255)
 tttACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCa                                                                                                 >  2:2703387/1‑100 (MQ=255)
  ttACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCAc                                                                                                >  2:3276942/1‑100 (MQ=255)
   tACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCAcc                                                                                               >  1:2462601/1‑100 (MQ=255)
    aCCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCa                                                                                              >  2:1590082/1‑100 (MQ=255)
     ccGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCaa                                                                                             >  1:4258330/1‑100 (MQ=255)
        ggAAAGAGATCCATGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGAt                                                                                          >  2:2898849/1‑100 (MQ=255)
            agagATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGc                                                                                      <  2:1280399/100‑1 (MQ=255)
             gagaTCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAGAATGCCACTCGGCAGCGTCACCAAGATCTg                                                                                       >  1:4166702/1‑98 (MQ=255)
             gagaTCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCt                                                                                     <  2:3417939/100‑1 (MQ=255)
             gagaTCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCt                                                                                     >  2:1895128/1‑100 (MQ=255)
             gagaTCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCt                                                                                     >  2:3455325/1‑100 (MQ=255)
              agaTCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCtt                                                                                    >  2:3604456/1‑100 (MQ=255)
               gaTCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTc                                                                                   >  2:3400582/1‑100 (MQ=255)
                 tCCCTGCGTTGTTGCAGGCATTTGTCTCTAGCGGTTTTCGTTATCAGATAATGGATGTCAAAAAAATGCCACTCGGCCGCGCCACCAAGATCTGCTTCt                                                                                  >  2:1967038/1‑99 (MQ=255)
                  cccTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCtgt                                                                                <  2:747795/100‑1 (MQ=255)
                  cccTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCtgt                                                                                <  1:2324965/100‑1 (MQ=255)
                   ccTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCtgtg                                                                               >  2:1426977/1‑100 (MQ=255)
                     tGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGc                                                                             >  2:3307898/1‑100 (MQ=255)
                     tGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGc                                                                             <  2:3510268/100‑1 (MQ=255)
                     tGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGc                                                                             <  1:1988041/100‑1 (MQ=255)
                     tGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGc                                                                             <  1:756400/100‑1 (MQ=255)
                       cGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGa                                                                           <  1:3283777/100‑1 (MQ=255)
                        gttgttGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGAt                                                                          <  2:3261371/100‑1 (MQ=255)
                         ttgttgCAGTCATTTGTCTATAGTGGTTTTCGTTATCAGAGAAAAAAAGAAAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGAt                                                                          <  1:3796717/99‑1 (MQ=255)
                         ttgttgCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATc                                                                         <  2:685008/100‑1 (MQ=255)
                           gttgCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCac                                                                       >  1:2997309/1‑100 (MQ=255)
                            ttgCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCacg                                                                      >  1:3866488/1‑100 (MQ=255)
                             tgCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCacga                                                                     >  1:3680879/1‑100 (MQ=255)
                              gCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCacgac                                                                    >  2:402666/1‑100 (MQ=255)
                               cAGGTACTTTTCTATATCGGTTTTCGTTATCAGATCACCGATGTCAAAAAAATGCCACTCGGCCGCGTCACCAAGATCTGCTTCTGTGGCGATCacgac                                                                    <  1:1647264/99‑1 (MQ=255)
                               cAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCacgacg                                                                   >  2:3035005/1‑100 (MQ=255)
                                aGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAAAAAAAAAAAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCacgac                                                                    <  1:1498103/98‑1 (MQ=255)
                                 ggCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATGACGACGAt                                                                 >  1:1790662/1‑100 (MQ=255)
                                  gCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGCTCACGACGATc                                                                >  1:631693/1‑100 (MQ=255)
                                   cATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCt                                                               <  2:4086796/100‑1 (MQ=255)
                                    aTTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCtt                                                              <  1:429622/100‑1 (MQ=255)
                                       tGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTTGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTaca                                                           <  2:1767947/100‑1 (MQ=255)
                                       tGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTaca                                                           <  2:3797648/100‑1 (MQ=255)
                                        gTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTacac                                                          >  2:2511085/1‑100 (MQ=255)
                                        gTCTATAGCGGTTTGCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTaca                                                           <  2:2969550/99‑1 (MQ=255)
                                          cTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGTGCCACCAAGATCTGCTTCTGTGGCGATCCCGACGATCCTACACGc                                                        >  1:3934685/1‑100 (MQ=255)
                                          cTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGc                                                        <  1:3786757/100‑1 (MQ=255)
                                           tataGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCt                                                       <  2:267439/100‑1 (MQ=255)
                                           tataGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCt                                                       >  1:2208958/1‑100 (MQ=255)
                                           tataGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCt                                                       <  1:239353/100‑1 (MQ=255)
                                           tataGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCt                                                       <  2:210235/100‑1 (MQ=255)
                                           tataGCGGTTTTCGTTATCAGATATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCtt                                                      >  1:3391805/1‑100 (MQ=255)
                                            ataGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCTACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGCCGATCTTACACACtt                                                      >  1:548584/1‑100 (MQ=255)
                                             taGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTg                                                     >  2:2059800/1‑100 (MQ=255)
                                               gCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGc                                                    <  2:777271/99‑1 (MQ=255)
                                               gCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCa                                                   <  2:1761167/100‑1 (MQ=255)
                                                  nttttCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACAGGCTTGCAGAt                                                >  1:2896748/2‑100 (MQ=255)
                                                    tttCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATcc                                              >  2:2484431/1‑100 (MQ=255)
                                                     ttCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTCTGGCGATCACGACGATCTTACACGCTTGCAGATCCa                                             <  2:3107129/100‑1 (MQ=255)
                                                      tCGTCATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATcc                                              <  1:3719484/98‑1 (MQ=255)
                                                       cGTTATCAGATAATCGATGTCAAAAAATTGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGc                                           <  1:3318913/100‑1 (MQ=255)
                                                       cGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGc                                           <  2:1798252/100‑1 (MQ=255)
                                                        gTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGCTCTTACCCGCTTGCAGATCCAGc                                           >  1:830437/1‑99 (MQ=255)
                                                        gTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCt                                          >  1:1051454/1‑100 (MQ=255)
                                                         ttATCAGATAATCGATGTCAAAAAAATGTCACTGGGCGGCGTCGCCAAGATCTGCTTGTGCGGCGATGAGGACGATCTTACACGCTTGCAGGTCCAGCta                                         >  1:3091077/1‑100 (MQ=255)
                                                         ttATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCta                                         <  2:3680904/100‑1 (MQ=255)
                                                          tATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCtat                                        >  1:230094/1‑100 (MQ=255)
                                                           aTCAGATAATCGATGTGAAAAAAATCCCACTCGTCAGCGTCACCAAGATCTGCTTCTGTGGCGATCAAGACGACAATACAGGCTTGCATATCCAGCtata                                       >  1:579535/1‑100 (MQ=255)
                                                           aTCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCtata                                       <  1:1590075/100‑1 (MQ=255)
                                                            tCAGATTATCGATTTCAAAAAACTTCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACTACGATCTTACACGCTTGCAGATCCAGCta                                         <  1:1204611/97‑1 (MQ=255)
                                                            tCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATAc                                      >  1:1903637/1‑100 (MQ=255)
                                                             cAGATAATCGATGTCAAAAAAAGGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACg                                     >  1:2559236/1‑100 (MQ=255)
                                                              aGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGa                                    >  1:4036826/1‑100 (MQ=255)
                                                               gctAATCGATGTCAAAAAAATGCCACTCTGCAGCGCCACCAAGATCTACTTCTGTGGCGATCACGACGATCCTACACCTTTGCAGATCCAGCTATAcaa                                    >  1:3351451/3‑97 (MQ=255)
                                                               gATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGaa                                   <  2:869158/100‑1 (MQ=255)
                                                               gATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGaa                                   <  2:2369007/100‑1 (MQ=255)
                                                                 tAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGc                                 >  2:834768/1‑100 (MQ=255)
                                                                 tAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGc                                 <  2:1660417/100‑1 (MQ=255)
                                                                 tAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGc                                 <  2:1214261/100‑1 (MQ=255)
                                                                  aaTCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCa                                >  1:4382904/1‑100 (MQ=255)
                                                                     cGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTa                             <  2:2103210/100‑1 (MQ=255)
                                                                       aTGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAgg                           <  2:1503688/100‑1 (MQ=255)
                                                                        tGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAg                            >  1:2763654/1‑98 (MQ=255)
                                                                        tGTCAAAAAAATGCCACTCGGCAGCGCCACCACGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGc                          >  2:593312/1‑100 (MQ=255)
                                                                         gTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCg                         <  2:2858559/100‑1 (MQ=255)
                                                                         gTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCg                         <  2:3881722/100‑1 (MQ=255)
                                                                          tCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGa                        >  1:3551484/1‑100 (MQ=255)
                                                                           cAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGGTCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAg                            >  2:10506/1‑95 (MQ=255)
                                                                           cAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAg                       >  1:2503964/1‑100 (MQ=255)
                                                                           cAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAg                       >  2:529684/1‑100 (MQ=255)
                                                                             aaaaaaTGCCTCTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCg                     <  1:743109/100‑1 (MQ=255)
                                                                             aaaaaaTGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCg                     >  2:293841/1‑100 (MQ=255)
                                                                              aaaaaTGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGt                    <  2:28071/100‑1 (MQ=255)
                                                                              aaaaaTGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACCACGATCTTACACGCTTGCAGATCCAGCTATACGCAGCATTAGGCGAGCGt                    >  1:924692/1‑100 (MQ=255)
                                                                               aaaaTGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTg                   >  2:3205870/1‑100 (MQ=255)
                                                                               aaaaTGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTg                   <  1:1093612/100‑1 (MQ=255)
                                                                                     ccACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACAtt             <  2:3952/100‑1 (MQ=255)
                                                                                      cACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACAttt            >  2:722011/1‑100 (MQ=255)
                                                                                         tCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTtgtg         <  2:328512/100‑1 (MQ=255)
                                                                                          cGGCAGCTTCACCAAAATCTGCTTCTGTTGCGATCACGGCTATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGTGACCGTGCACAt              >  1:264055/1‑94 (MQ=255)
                                                                                           ggCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTTGTGtt       <  2:4191905/100‑1 (MQ=255)
                                                                                           ggCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTTGTGtt       >  2:2553996/1‑100 (MQ=255)
                                                                                               gCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTTGTGTTTTTc   <  2:3936747/100‑1 (MQ=255)
                                                                                                cGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTTGTGTTTTTcc  >  2:2229637/1‑100 (MQ=255)
                                                                                                  |                                                                                                 
GTTTACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTTGTGTTTTTCC  >  NC_000913/467741‑467936

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