Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,579,447 T→C S598P (TCG→CCG)  tktB → transketolase 2 isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,579,4470TC99.6% 797.0 / ‑5.8 261S598P (TCG→CCG) tktBtransketolase 2 isozyme
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (133/126);  total (134/127)
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.84e-01

Only 100 of 270 total aligned reads displayed.

GAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCA  >  NC_000913/2579357‑2579536
                                                                                          |                                                                                         
gagaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCg                                                                                         >  1:1183193/1‑93 (MQ=255)
gagaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCg                                                                                         >  1:5809565/1‑93 (MQ=255)
gagaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCg                                                                                         >  1:7598696/1‑93 (MQ=255)
gagaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCg                                                                                         <  2:77774/93‑1 (MQ=255)
 agaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCgg                                                                                        >  1:1222485/1‑93 (MQ=255)
  gaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGgt                                                                                       <  1:3954485/93‑1 (MQ=255)
  gaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGgt                                                                                       >  1:6825422/1‑93 (MQ=255)
  gaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTAGACGCCCAGGATGAGGAATATCGGGAGCCggg                                                                                       >  2:2322074/1‑92 (MQ=255)
   aaaTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGgtg                                                                                      <  1:6896177/93‑1 (MQ=255)
    aaTTAGCAGGAGAAGGTCGCAATGTACGCGTGGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGgtgt                                                                                     <  1:1824292/93‑1 (MQ=255)
    aaTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGgtgt                                                                                     <  1:2288513/93‑1 (MQ=255)
     aTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGtt                                                                                    <  2:227601/93‑1 (MQ=255)
       tAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGc                                                                                  <  1:7343412/93‑1 (MQ=255)
        aGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGcc                                                                                 <  2:6602817/93‑1 (MQ=255)
        aGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCGGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGCGCCGGTGTTGcc                                                                                 >  2:1510013/1‑93 (MQ=255)
          cAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCtt                                                                               <  1:3047702/93‑1 (MQ=255)
           aGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTc                                                                              >  2:583486/1‑93 (MQ=255)
            ggAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCt                                                                             >  1:4987789/1‑93 (MQ=255)
                 aGGTCGCAATGTGCGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGt                                                                        <  2:1465484/93‑1 (MQ=255)
                 aGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGt                                                                        <  2:6944193/93‑1 (MQ=255)
                  ggTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGtt                                                                       <  1:4125330/93‑1 (MQ=255)
                    tCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGc                                                                     <  2:1134387/93‑1 (MQ=255)
                      gCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCgg                                                                   >  1:3327446/1‑93 (MQ=255)
                       cAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGc                                                                  <  2:571798/93‑1 (MQ=255)
                        aaTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCt                                                                 >  2:268710/1‑93 (MQ=255)
                        aaTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCt                                                                 >  2:2933507/1‑93 (MQ=255)
                         aTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTc                                                                >  2:2131558/1‑93 (MQ=255)
                         aTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTc                                                                <  1:7284714/93‑1 (MQ=255)
                         aTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTc                                                                <  1:706993/93‑1 (MQ=255)
                         aTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTc                                                                <  2:851131/93‑1 (MQ=255)
                          tGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTcg                                                               >  2:5518367/1‑93 (MQ=255)
                           gTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTccc                                                              >  2:7312638/1‑91 (MQ=255)
                            tACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTcgcg                                                             <  2:4584790/93‑1 (MQ=255)
                             aCGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGt                                                            <  1:2918833/93‑1 (MQ=255)
                              cgcgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTg                                                           >  1:3466916/1‑93 (MQ=255)
                              cgcgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTg                                                           <  1:6597759/93‑1 (MQ=255)
                               gcgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTgg                                                          <  2:6311736/93‑1 (MQ=255)
                                cgTAGTTTCCCTGCCCTCGATCGATATTATCGACGCCCAGGATGAGGAATATCGGGAGCCTGTGTTGCCTTCTCACGTTGCGGCTCGCGTGGc                                                         <  2:7936276/93‑1 (MQ=255)
                                cgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGc                                                         >  1:5628864/1‑93 (MQ=255)
                                cgTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGc                                                         <  1:7395970/93‑1 (MQ=255)
                                 gTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCg                                                        >  1:3394455/1‑93 (MQ=255)
                                  tAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGgggggg                                                       >  2:6690552/1‑93 (MQ=255)
                                  tAGTTTCCCTGCCCTCGACAGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGGGGCgg                                                       >  1:448981/1‑93 (MQ=255)
                                    gTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTg                                                     <  2:6187621/93‑1 (MQ=255)
                                    gTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTg                                                     >  1:3241297/1‑93 (MQ=255)
                                     tttGCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTgg                                                    <  1:6275793/93‑1 (MQ=255)
                                     tttCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTgg                                                    <  1:1920047/93‑1 (MQ=255)
                                     tttCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTgg                                                    <  1:2487128/93‑1 (MQ=255)
                                       tCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGaa                                                  <  1:7640555/93‑1 (MQ=255)
                                       tCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGaa                                                  <  2:5447737/93‑1 (MQ=255)
                                        cccTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAg                                                 <  1:5105407/93‑1 (MQ=255)
                                         ccTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCCCGGGGTGGGGGGAGc                                                >  2:3303434/1‑93 (MQ=255)
                                          cTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCa                                               >  1:5514867/1‑93 (MQ=255)
                                           tGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAg                                              <  2:5012990/93‑1 (MQ=255)
                                           tGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAggcg                                              >  2:256525/1‑90 (MQ=255)
                                            gCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAgg                                             <  2:5761132/93‑1 (MQ=255)
                                             gccTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGt                                            <  2:2234736/92‑1 (MQ=255)
                                              ccTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTa                                           <  1:5604764/93‑1 (MQ=255)
                                              ccTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAgggt                                           >  1:3796320/1‑91 (MQ=255)
                                               cTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTAt                                          <  2:5137413/93‑1 (MQ=255)
                                               cTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTAt                                          <  1:4649619/93‑1 (MQ=255)
                                                tCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTAtt                                         <  1:709320/93‑1 (MQ=255)
                                                tCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGAtt                                         >  2:4762242/1‑93 (MQ=255)
                                                tCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAGGGAtt                                         >  2:1771959/1‑93 (MQ=255)
                                                tCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCCCGGGGGGGGGGGAGGAGGGAtt                                         >  1:1514290/1‑93 (MQ=255)
                                                  gACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGc                                       <  1:5855845/93‑1 (MQ=255)
                                                   aCCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGCGTTGCCTTCTAACGATGCGGCTCCCGTCGCGGCGGAAACAGGTGtttcg                                      >  1:4246023/1‑90 (MQ=255)
                                                   aCCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGcc                                      <  1:2072624/93‑1 (MQ=255)
                                                     cGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGGATTGCCGa                                    >  1:1661565/1‑93 (MQ=255)
                                                     cGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGa                                    <  1:1717538/93‑1 (MQ=255)
                                                     cGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGa                                    <  2:2735420/93‑1 (MQ=255)
                                                      gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCCGGGATTTCCGAt                                   >  1:810319/1‑93 (MQ=255)
                                                      gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGAAGCAGGGGTTTCCGAt                                   >  2:1422184/1‑93 (MQ=255)
                                                       atatTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTTCCGAtt                                  >  1:2536448/1‑93 (MQ=255)
                                                        tatTTTCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGTATTGCCGATTa                                 >  2:7325568/1‑93 (MQ=255)
                                                          ttttCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGAAGCAGGTATTGCCGGTTACt                               >  1:1010304/1‑93 (MQ=255)
                                                           tttCGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTTCCGATTACTg                              >  2:5062260/1‑93 (MQ=255)
                                                              cGACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTa                           <  2:6569618/93‑1 (MQ=255)
                                                               gACGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTAc                          <  1:837278/93‑1 (MQ=255)
                                                                 cGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACaa                        <  2:4141499/93‑1 (MQ=255)
                                                                 cGCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACaa                        <  2:7598695/93‑1 (MQ=255)
                                                                  gCCCAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAg                       <  2:4905692/93‑1 (MQ=255)
                                                                     cAGGATGAGGAATATCGGGCGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGGGGTGGAAGCGGGTCTTGCCCGTTACTGGGACAAGTAt                    >  2:3445116/1‑93 (MQ=255)
                                                                     cAGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTGCCCATTACCGGGACAAGTAt                    >  2:1074872/1‑93 (MQ=255)
                                                                      aGGATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTGCCGATTACTGGTACAAGTATg                   >  2:4591561/1‑93 (MQ=255)
                                                                        gATGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTGCCGATTACTGGTACAAGTATGtt                 >  1:883633/1‑93 (MQ=255)
                                                                         aTGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGGTTACTGGGACAAGTATGGTg                >  2:1572664/1‑93 (MQ=255)
                                                                          tGAGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCCTGGCGGTGGAAGCAGGGATTGCCGGTTACTGGTACAAGGATGTTgg               >  2:7073275/1‑93 (MQ=255)
                                                                          tGAGGAATATCAGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTgg               >  1:4350864/1‑93 (MQ=255)
                                                                            aGGAATATCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTc             <  1:2832959/93‑1 (MQ=255)
                                                                                 tatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaaa        <  2:5301266/93‑1 (MQ=255)
                                                                                 tatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaaa        <  1:7753398/93‑1 (MQ=255)
                                                                                 tatCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGGATTGCCGCTTACTGGTACAAGTATGTTGGTCTGaaa        >  1:7454018/1‑93 (MQ=255)
                                                                                  atCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAg       <  1:5975540/93‑1 (MQ=255)
                                                                                   tCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGGATTGCCGATTACTGGTACCAGGATGTTGGTCTGgaaag      >  2:7811760/1‑91 (MQ=255)
                                                                                   tCGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGGGGGAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAgg      >  2:3630194/1‑93 (MQ=38)
                                                                                    cGGGAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAggg     <  1:6099260/93‑1 (MQ=255)
                                                                                     gggAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAgggg    <  2:5912168/93‑1 (MQ=255)
                                                                                       gAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGGGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTTGTCTGAAAGGGGCa  >  2:4475514/1‑93 (MQ=255)
                                                                                       gAGCCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGGAGCAGGTATTGCCGATTACTTGTACAAGTATGTTGGTCTGAAAGGGGCa  >  1:5195108/1‑93 (MQ=255)
                                                                                          |                                                                                         
GAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCA  >  NC_000913/2579357‑2579536

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