Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,971,075 2 bp→GC coding (942‑943/1068 nt) rffG → dTDP‑glucose 4,6‑dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,971,0750CG100.0% 801.5 / NA 246C314W (TGC→TGG) rffGdTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (120/126);  total (120/126)
*NC_0009133,971,0760GC99.6% 801.7 / NA 246V315L (GTG→CTG) rffGdTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base C (119/126);  total (120/126)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 255 total aligned reads displayed.

TGACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGCGTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGAC  >  NC_000913/3970988‑3971162
                                                                                       ||                                                                                      
   ccTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCa                                                                                 <  1:1312511/93‑1 (MQ=255)
   ccTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCa                                                                                 >  2:3637679/1‑93 (MQ=255)
    cTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAg                                                                                <  2:6186629/93‑1 (MQ=255)
    cTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAg                                                                                >  2:2768447/1‑93 (MQ=255)
      gATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGa                                                                              >  1:852073/1‑93 (MQ=255)
      gATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGa                                                                              >  1:5644079/1‑93 (MQ=255)
      gATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGa                                                                              >  1:4778983/1‑93 (MQ=255)
        tCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGaaa                                                                            >  2:678669/1‑93 (MQ=255)
          aCCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAccc                                                                          >  1:5612783/1‑93 (MQ=255)
          aCCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAAcc                                                                          <  2:5928549/93‑1 (MQ=255)
           ccTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCt                                                                         <  2:2419207/93‑1 (MQ=255)
            cTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCtt                                                                        <  2:5764671/93‑1 (MQ=255)
             tttGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCttt                                                                       <  1:2177657/93‑1 (MQ=255)
             tttGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCttt                                                                       <  2:5062565/93‑1 (MQ=255)
             tttGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCttt                                                                       >  2:4036604/1‑93 (MQ=255)
             tttGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCttt                                                                       <  2:3947293/93‑1 (MQ=255)
             tttGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCttt                                                                       >  1:2510092/1‑93 (MQ=255)
             tttGTCGCAGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCttt                                                                       >  1:6523522/1‑93 (MQ=255)
              ttGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTg                                                                      <  1:2073987/93‑1 (MQ=255)
               tGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGa                                                                     >  1:5406322/1‑93 (MQ=255)
                gccGCTGACCGCCCGGGGCATGATCTGCGTTATTCCATTGATGCTACGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGaa                                                                    <  2:6474135/91‑1 (MQ=255)
                gTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGaa                                                                    <  1:4237633/93‑1 (MQ=255)
                 tCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGaaa                                                                   <  2:4904742/93‑1 (MQ=255)
                   gCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGt                                                                 >  1:5839409/1‑93 (MQ=255)
                    cTGAGCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTg                                                                <  1:3743707/93‑1 (MQ=255)
                    cTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTg                                                                <  2:1578848/93‑1 (MQ=255)
                     tGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTgg                                                               <  1:634339/93‑1 (MQ=255)
                      gACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGa                                                              <  1:2672398/93‑1 (MQ=255)
                      gACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGa                                                              >  1:1701125/1‑93 (MQ=255)
                       gccGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGaa                                                             <  2:3158393/92‑1 (MQ=255)
                       aCCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGaa                                                             <  2:4560898/93‑1 (MQ=255)
                        ccGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAAt                                                            >  1:5790723/1‑93 (MQ=255)
                          gTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGc                                                          <  1:5920896/93‑1 (MQ=255)
                          gTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGc                                                          <  1:2568238/93‑1 (MQ=255)
                           tCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCg                                                         >  1:4316222/1‑93 (MQ=255)
                             cGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTa                                                       <  2:3623297/93‑1 (MQ=255)
                                ggCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTaaaa                                                    <  1:6257034/93‑1 (MQ=255)
                                 gCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAAc                                                   >  1:5508205/1‑93 (MQ=255)
                                 gCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAAc                                                   <  2:5651163/93‑1 (MQ=255)
                                 gCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAAc                                                   <  1:4535506/93‑1 (MQ=255)
                                 gCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAAc                                                   <  2:3524800/93‑1 (MQ=255)
                                  cATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACg                                                  >  1:1413051/1‑93 (MQ=255)
                                   aTGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACgg                                                 >  1:3877673/1‑93 (MQ=255)
                                    tGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGt                                                >  1:6490193/1‑93 (MQ=255)
                                     gATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGATATTGGAATGCGTTAAACGGtt                                               <  1:6565746/93‑1 (MQ=255)
                                     gATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGtt                                               >  2:1884279/1‑93 (MQ=255)
                                     gATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGtt                                               >  1:3610592/1‑93 (MQ=255)
                                      aTCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTc                                              <  2:5812546/93‑1 (MQ=255)
                                       tCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCa                                             >  2:5881073/1‑93 (MQ=255)
                                         tGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGt                                           >  1:6640367/1‑93 (MQ=255)
                                         tGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGt                                           <  2:4413454/93‑1 (MQ=255)
                                          gCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTg                                          <  1:190122/93‑1 (MQ=255)
                                            gTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGt                                        <  2:3144979/93‑1 (MQ=255)
                                            gTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGt                                        >  1:3797861/1‑93 (MQ=255)
                                             ttATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTa                                       <  2:159104/93‑1 (MQ=255)
                                             ttATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTa                                       <  1:3012238/93‑1 (MQ=255)
                                             ttATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTa                                       >  2:5304006/1‑93 (MQ=255)
                                               aTGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATc                                     >  2:1432033/1‑93 (MQ=255)
                                               aTGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATc                                     >  2:1196353/1‑93 (MQ=255)
                                               aTGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATc                                     <  2:19715/93‑1 (MQ=255)
                                               aTGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATc                                     >  1:6201603/1‑93 (MQ=255)
                                               aTGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATc                                     <  2:2588341/93‑1 (MQ=255)
                                               aTGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATc                                     <  1:1877002/93‑1 (MQ=255)
                                                 gCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTg                                   <  1:286344/93‑1 (MQ=255)
                                                 gCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGACGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTg                                   >  2:3729401/1‑93 (MQ=255)
                                                  ccATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTgg                                  >  2:4485668/1‑93 (MQ=255)
                                                  ccATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTgg                                  >  2:1480741/1‑93 (MQ=255)
                                                  ccATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTgg                                  <  1:3395085/93‑1 (MQ=255)
                                                   cATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGc                                 <  2:3478868/93‑1 (MQ=255)
                                                   cATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGc                                 <  2:1823718/93‑1 (MQ=255)
                                                   cATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGc                                 <  1:4415531/93‑1 (MQ=255)
                                                    aTTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGGTTCTGGcg                                >  1:4012537/1‑92 (MQ=255)
                                                     ttGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTa                               >  2:108377/1‑93 (MQ=255)
                                                       gATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAAt                             >  1:213122/1‑93 (MQ=255)
                                                       gATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAAt                             <  2:2808785/93‑1 (MQ=255)
                                                        aTGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATg                            <  2:6620026/93‑1 (MQ=255)
                                                          gCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGaa                          >  1:5938553/1‑93 (MQ=255)
                                                           cTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGaaa                         <  1:643690/93‑1 (MQ=255)
                                                            ttCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAg                        >  1:5067396/1‑93 (MQ=255)
                                                             tCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGc                       <  1:484043/93‑1 (MQ=255)
                                                               gAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtg                     >  1:1102733/1‑93 (MQ=255)
                                                               gAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtg                     <  2:1738094/93‑1 (MQ=255)
                                                                 aaaTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtggt                   >  2:5610322/1‑93 (MQ=255)
                                                                   aTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtggtgg                 >  2:5382666/1‑93 (MQ=255)
                                                                     tGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGTGTAACACGGTTCAGTGGCATCTTGCTAAAGAAAGCTGGCGGaa               >  1:2012355/1‑93 (MQ=255)
                                                                      gCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAg              >  2:618411/1‑93 (MQ=255)
                                                                        ccGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCa            >  1:4737943/1‑93 (MQ=255)
                                                                        ccGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCa            >  2:608629/1‑93 (MQ=255)
                                                                        ccGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCa            >  2:4945518/1‑93 (MQ=255)
                                                                        ccGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCa            <  2:2030704/93‑1 (MQ=255)
                                                                        ccGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCa            >  2:3419865/1‑93 (MQ=255)
                                                                         cGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAg           >  1:4163032/1‑93 (MQ=255)
                                                                         cGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTAGCTAATGAAAGCTGGTGGAAGCAg           <  2:5195671/93‑1 (MQ=255)
                                                                          gTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAgg          >  1:6671421/1‑93 (MQ=255)
                                                                          gTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCATAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAgg          <  1:752780/93‑1 (MQ=255)
                                                                           tGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGt         <  2:4222130/93‑1 (MQ=255)
                                                                           tGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGt         >  2:381254/1‑93 (MQ=255)
                                                                               cTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAg     <  1:5153645/93‑1 (MQ=255)
                                                                                ttGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAgg    >  2:1343630/1‑93 (MQ=255)
                                                                                 tggctggcTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAGACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGa   <  2:4721758/93‑1 (MQ=255)
                                                                                       ||                                                                                      
TGACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGCGTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGAC  >  NC_000913/3970988‑3971162

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