Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 467,839 T→C 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,8390TC98.3% 400.7 / ‑2.7 116V143A (GTC→GCC) cofthiamine pyrimidine pyrophosphate hydrolase; HMP‑PP phosphatase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (61/53);  total (62/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 117 total aligned reads displayed.

TACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTTGTGTTTTT  >  NC_000913/467744‑467934
                                                                                               |                                                                                               
tACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCAcc                                                                                             <  1:1320216/100‑1 (MQ=255)
 aCCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGCTAGACGATTTCTAAAAAATGCCACTCGGCAGCGCCAcc                                                                                             >  2:2409794/1‑99 (MQ=255)
 aCCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGccacca                                                                                            >  2:2555903/1‑100 (MQ=255)
  ccGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCaa                                                                                           <  1:703910/100‑1 (MQ=255)
    gggAAAGAGATCCCTGCGTTGTTGTAGTTAATTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAATATAATGCCACTCGGCAGCGCCACCAAGa                                                                                         <  2:676813/100‑1 (MQ=255)
    gggAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGa                                                                                         >  1:806949/1‑100 (MQ=255)
    gggAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAg                                                                                          >  1:1145606/1‑99 (MQ=255)
    gggAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGa                                                                                         >  1:2861495/1‑100 (MQ=255)
     ggAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGAt                                                                                        <  1:1457510/100‑1 (MQ=255)
      gAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATc                                                                                       >  2:97788/1‑100 (MQ=255)
      gAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATc                                                                                       <  2:96158/100‑1 (MQ=255)
      gAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATc                                                                                       <  1:1835998/100‑1 (MQ=255)
        aaGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTTGGCAGCGCCACCAAGATCTg                                                                                     <  1:2519009/100‑1 (MQ=255)
        aaGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTg                                                                                     <  1:706925/100‑1 (MQ=255)
          gagaTCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCt                                                                                   <  1:2214219/100‑1 (MQ=255)
          gagaTCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCt                                                                                   >  2:986617/1‑100 (MQ=255)
          gagaTCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCt                                                                                   <  1:2219024/100‑1 (MQ=255)
          gagaTCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCt                                                                                   >  2:1692706/1‑100 (MQ=255)
           agaTCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTTGGCAGCGCCACCAAGATCTGCtt                                                                                  >  2:557102/1‑100 (MQ=255)
              tCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCtg                                                                               <  2:1509610/100‑1 (MQ=255)
               cccTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCtgtg                                                                             <  2:736375/100‑1 (MQ=255)
                ccTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCtgtg                                                                             >  1:1087484/1‑100 (MQ=255)
                ccTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCtgtg                                                                             >  1:635235/1‑100 (MQ=255)
                  tGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCTAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCtgtg                                                                             <  2:2309177/98‑1 (MQ=255)
                  tGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGc                                                                           <  2:1268285/100‑1 (MQ=255)
                   gCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCg                                                                          >  2:1590917/1‑100 (MQ=255)
                    cGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGa                                                                         <  2:745898/100‑1 (MQ=255)
                    cGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGa                                                                         <  2:745886/100‑1 (MQ=255)
                    cGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGa                                                                         >  1:2542194/1‑100 (MQ=255)
                     gttgttGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCCGCTTCTGTGGCGAt                                                                        >  2:2070950/1‑100 (MQ=255)
                     gttgttGCAGGCACTTGTATATAGCGGTTTTCGTTATCAGATAAAAAAAAAAAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGAt                                                                        <  2:2826810/100‑1 (MQ=255)
                     gttgttGAAGGCATTTGTCTATAGCGGTTTTCGTTATCCGATAACCTATGTCAAAAAAATGCCAATCGGCAGCGCCACCAAGATCTGCATCTATGACgct                                                                        >  2:565316/1‑98 (MQ=255)
                      ttgttgCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATc                                                                       <  1:364418/100‑1 (MQ=255)
                      ttgttgCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATc                                                                       >  1:924476/1‑100 (MQ=255)
                       tgttgCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCa                                                                      >  2:1870866/1‑100 (MQ=255)
                        gttgCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGCTCac                                                                     >  2:764717/1‑100 (MQ=255)
                        gttgCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCac                                                                     <  1:371771/100‑1 (MQ=255)
                        gttgCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCac                                                                     >  2:2316885/1‑100 (MQ=255)
                          tgCAGGCATTTGGCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCacga                                                                   >  1:1896475/1‑100 (MQ=255)
                           gCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCacgac                                                                  >  2:1732746/1‑100 (MQ=255)
                            cAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCacgacg                                                                 >  2:106744/1‑100 (MQ=255)
                                cATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCt                                                             >  2:1519257/1‑100 (MQ=255)
                                cATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCt                                                             <  1:753854/100‑1 (MQ=255)
                                cATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCt                                                             >  2:2165160/1‑100 (MQ=255)
                                   ttGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTac                                                          <  2:322471/100‑1 (MQ=255)
                                     gTCTATAGCGGTTTTCGCTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTacac                                                        >  2:2126346/1‑100 (MQ=255)
                                       cgataGCGGTTTTAGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGc                                                      >  1:1959329/3‑100 (MQ=255)
                                       cTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGc                                                      <  2:1517588/100‑1 (MQ=255)
                                        tataGCGGTTTTCGTTATCAGATACTCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCt                                                     >  2:1286131/1‑100 (MQ=255)
                                         ataGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCtt                                                    >  2:2342876/1‑100 (MQ=255)
                                          taGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTg                                                   <  2:1119059/100‑1 (MQ=255)
                                          taGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTg                                                   >  1:766053/1‑100 (MQ=255)
                                          taGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTg                                                   <  1:2072491/100‑1 (MQ=255)
                                           aGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGc                                                  >  1:1678304/1‑100 (MQ=255)
                                           aGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGc                                                  <  2:924474/100‑1 (MQ=255)
                                            gCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCa                                                 >  2:1674132/1‑100 (MQ=255)
                                            gCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCa                                                 >  2:419702/1‑100 (MQ=255)
                                              ggTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGa                                               <  2:529508/100‑1 (MQ=255)
                                               gTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGAt                                              >  1:819769/1‑100 (MQ=255)
                                               gTTTTCGTTATCAGAGAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTTCTTCTGTGGCGATCACGACGATCTTACACGCTTGCatat                                              >  2:2305665/1‑100 (MQ=255)
                                                 tttCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATcc                                            >  1:818745/1‑100 (MQ=255)
                                                  ttCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCa                                           <  2:2419419/100‑1 (MQ=255)
                                                   tCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAg                                          <  1:2185154/100‑1 (MQ=255)
                                                   tCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAg                                          <  1:2387916/100‑1 (MQ=255)
                                                    cGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGc                                         <  2:2518784/100‑1 (MQ=255)
                                                    cGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGc                                         <  1:315176/100‑1 (MQ=255)
                                                     gTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCt                                        >  1:2208518/1‑100 (MQ=255)
                                                     gTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCt                                        <  1:394881/100‑1 (MQ=255)
                                                       tATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCtat                                      >  2:1924098/1‑100 (MQ=255)
                                                        aTCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCtata                                     >  1:519967/1‑100 (MQ=255)
                                                        aTCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCtata                                     >  1:2771077/1‑100 (MQ=255)
                                                         tCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATAc                                    >  1:1315976/1‑100 (MQ=255)
                                                         tCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATAc                                    >  1:1456591/1‑100 (MQ=255)
                                                          cAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACTCGCTTGCAGATCCCGCTCTACg                                   >  2:222913/1‑100 (MQ=255)
                                                          cAGATAATCGATGACAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACg                                   >  1:1959735/1‑100 (MQ=255)
                                                            gATAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGaa                                 <  1:1272878/100‑1 (MQ=255)
                                                              tAATCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGc                               <  2:1484149/100‑1 (MQ=255)
                                                                 tCGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCAtt                            <  2:614660/100‑1 (MQ=255)
                                                                  cGATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTa                           <  2:2494614/100‑1 (MQ=255)
                                                                   gATGTCAATAAAATGCCACTCGGCAGCGCCACCAAGAGCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAg                          <  2:2888615/100‑1 (MQ=255)
                                                                   gATGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAg                          >  1:2504235/1‑100 (MQ=255)
                                                                    aTGTCAAAAACATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACTCTTGCAGATCCAGCTATACGAAGCATTAgg                         <  2:1896462/100‑1 (MQ=255)
                                                                    aTGTCAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAgg                         >  1:2070904/1‑100 (MQ=255)
                                                                        cAAAAAAATGCCACTCGGCAGCGCTACCAAGGTCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAg                     <  1:1579132/100‑1 (MQ=255)
                                                                        cAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAg                     >  1:1877001/1‑100 (MQ=255)
                                                                        cAAAAAAATGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAg                     >  1:1220385/1‑100 (MQ=255)
                                                                         aaaaaaaTGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGc                    >  2:2173048/1‑100 (MQ=255)
                                                                          aaaaaaTGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCg                   <  2:336019/100‑1 (MQ=255)
                                                                          aaaaaaTGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCg                   <  2:2771437/100‑1 (MQ=255)
                                                                          aaaaaaTGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCg                   >  1:436047/1‑100 (MQ=255)
                                                                           aaaaaTGCCACTCGGCCGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCg                   <  1:2305679/99‑1 (MQ=255)
                                                                                tGCCACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGcaca             <  2:2189840/100‑1 (MQ=255)
                                                                                  ccACTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCATCTATACGAAGCATTAGGCGAGCGTGCACAtt           >  1:2409235/1‑100 (MQ=255)
                                                                                     cTCGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTtgt        >  1:1570253/1‑100 (MQ=255)
                                                                                       cGGCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTTgtgt      <  2:1087481/100‑1 (MQ=255)
                                                                                        ggCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTTgtgt      >  2:1602082/1‑99 (MQ=255)
                                                                                        ggCAGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTTGTGtt     <  1:298550/100‑1 (MQ=255)
                                                                                           aGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTTGTGttttt  <  2:818744/100‑1 (MQ=255)
                                                                                           aGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTTGTGttttt  <  1:1590926/100‑1 (MQ=255)
                                                                                           aGCGCCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTTGTGttttt  <  1:53028/100‑1 (MQ=255)
                                                                                               |                                                                                               
TACCGGGAAAGAGATCCCTGCGTTGTTGCAGGCATTTGTCTATAGCGGTTTTCGTTATCAGATAATCGATGTCAAAAAAATGCCACTCGGCAGCGTCACCAAGATCTGCTTCTGTGGCGATCACGACGATCTTACACGCTTGCAGATCCAGCTATACGAAGCATTAGGCGAGCGTGCACATTTGTGTTTTT  >  NC_000913/467744‑467934

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