Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,250,644 A→C G65G (GGT→GGG dhaK ← dihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,250,6440AC97.8% 633.3 / 0.4 183G65G (GGT→GGGdhaKdihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base C (78/101);  total (81/102)
Fisher's exact test for biased strand distribution p-value = 3.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 185 total aligned reads displayed.

CTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATA  >  NC_000913/1250547‑1250739
                                                                                                 |                                                                                               
 ttCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCg                                                                                              <  2:2340981/100‑1 (MQ=255)
 ttCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCg                                                                                              <  1:4482819/100‑1 (MQ=255)
   cgccgcCATCATGGTGCAGGGCGCATTCAATGATTTTATCGGGCGTCGGTGAGGTGAAAATCTCGCCCGGACAGGCCCCCGAAAGCACCCCCTGCCCGat                                                                                            <  2:3815469/100‑1 (MQ=255)
   cgccgcCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGat                                                                                            <  1:740689/100‑1 (MQ=255)
   cgccgcCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGat                                                                                            <  2:4144264/100‑1 (MQ=255)
    gccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGatn                                                                                           <  1:4416890/100‑2 (MQ=255)
     ccgccATCAACTTGCATGGCGCATTCAAGGGTTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGatat                                                                                          <  2:3321685/100‑1 (MQ=255)
     ccgccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGatat                                                                                          <  1:1482705/100‑1 (MQ=255)
        ccATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAAc                                                                                       >  1:347798/1‑100 (MQ=255)
         cATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAAcc                                                                                      >  1:4357783/1‑100 (MQ=255)
          aTCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACca                                                                                     >  1:2161846/1‑100 (MQ=255)
           tCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACcac                                                                                    >  2:1951550/1‑100 (MQ=255)
           tCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACcac                                                                                    <  2:5150035/100‑1 (MQ=255)
           tCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACcac                                                                                    <  2:4666017/100‑1 (MQ=255)
            cAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACcaca                                                                                   >  2:4730288/1‑100 (MQ=255)
               cTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTg                                                                                >  1:2654295/1‑100 (MQ=255)
               cTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTg                                                                                <  2:1198885/100‑1 (MQ=255)
               cTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTg                                                                                >  1:3134453/1‑100 (MQ=255)
               cTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAATg                                                                                <  2:4093618/100‑1 (MQ=255)
                   cATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATc                                                                            <  1:5362605/100‑1 (MQ=255)
                   cATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATc                                                                            <  1:1757515/100‑1 (MQ=255)
                    aTGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCg                                                                           <  2:3254898/100‑1 (MQ=255)
                      ggCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCgcc                                                                         >  2:2171458/1‑98 (MQ=255)
                      ggCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGc                                                                         >  1:4745710/1‑100 (MQ=255)
                      ggCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGACATAACCACAGTGCATCGGc                                                                         >  1:2296754/1‑100 (MQ=255)
                        cgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTc                                                                       <  2:929539/100‑1 (MQ=255)
                        cgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTc                                                                       >  2:1840301/1‑100 (MQ=255)
                        cgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTc                                                                       <  1:4806608/100‑1 (MQ=255)
                        cgcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTc                                                                       <  2:5440468/100‑1 (MQ=255)
                         gcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCg                                                                      <  1:2462474/100‑1 (MQ=255)
                         gcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCg                                                                      <  1:5491774/100‑1 (MQ=255)
                         gcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCg                                                                      <  1:1973031/100‑1 (MQ=255)
                         gcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCg                                                                      <  1:242005/100‑1 (MQ=255)
                         gcATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCg                                                                      <  1:843458/100‑1 (MQ=255)
                           aTTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCgtg                                                                    >  1:4409539/1‑100 (MQ=255)
                            ttCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCgtgt                                                                   <  1:5227441/100‑1 (MQ=255)
                              cAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTcc                                                                 >  1:40988/1‑100 (MQ=255)
                              cAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTcc                                                                 <  1:1437105/100‑1 (MQ=255)
                              cAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTcc                                                                 >  2:4530380/1‑100 (MQ=255)
                               aaaGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCg                                                                <  1:4161398/100‑1 (MQ=255)
                               aaaGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCg                                                                <  1:2332273/100‑1 (MQ=255)
                               aaaGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCg                                                                >  1:1577310/1‑100 (MQ=255)
                                aaGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCACGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGCGTCCGc                                                               >  2:820237/1‑100 (MQ=255)
                                   aTTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgc                                                            >  2:3735775/1‑100 (MQ=255)
                                      ttATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                         >  1:4599492/1‑100 (MQ=255)
                                      ttATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgc                                                         <  2:4324973/100‑1 (MQ=255)
                                       tATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTgccgcc                                                        <  1:4693262/100‑1 (MQ=255)
                                           gggCCTCGGTGAGGTGAAAATTTCGCCCTGACAGCCCCCCGAAAGCAGCCCCTGCCCGAAATAACCACATTTCATCGGCTCACGGCCGCTCCCGCCAccc                                                    >  2:4664845/1‑100 (MQ=255)
                                            ggCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCg                                                   >  1:429157/1‑100 (MQ=255)
                                              cGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGAc                                                 <  1:3785587/100‑1 (MQ=255)
                                                tCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGAcag                                               <  1:68309/100‑1 (MQ=255)
                                                  ggtgaggtgaAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGAcagca                                             >  1:2215357/1‑100 (MQ=255)
                                                  ggtgaggtgaAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGAcagca                                             <  2:5287213/100‑1 (MQ=255)
                                                  ggtgaggtgaAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGAcagca                                             >  1:3977041/1‑100 (MQ=255)
                                                  ggtgaggtgaAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGAcagca                                             >  2:2950425/1‑100 (MQ=255)
                                                   gtgaggtgaAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGAcagcag                                            >  2:4563448/1‑100 (MQ=255)
                                                   gtgaggtgaAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGAcagcag                                            <  2:5391890/100‑1 (MQ=255)
                                                   gtgaggtgaAAATTTCGCCCGGACAGGCCCCCGAAAGCAACCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGAcagcag                                            >  1:5355862/1‑100 (MQ=255)
                                                    tgaggtgaAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAgg                                           <  1:2278219/100‑1 (MQ=255)
                                                       ggtgaAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCg                                        >  2:2642866/1‑100 (MQ=255)
                                                       ggtgaAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCg                                        >  2:2278895/1‑100 (MQ=255)
                                                        gtgaAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGa                                       >  2:2903246/1‑100 (MQ=255)
                                                        gtgaAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGa                                       >  1:2995628/1‑100 (MQ=255)
                                                         tgaAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGAc                                      <  2:2961679/100‑1 (MQ=255)
                                                         tgaAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGAc                                      <  1:831798/100‑1 (MQ=255)
                                                          gaAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACt                                     >  1:1085887/1‑100 (MQ=255)
                                                               tttCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTcc                                >  1:3063312/1‑100 (MQ=255)
                                                                ttCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCt                               >  1:327349/1‑100 (MQ=255)
                                                                ttCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCt                               <  1:5536905/100‑1 (MQ=255)
                                                                 tCGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTg                              <  2:1816666/100‑1 (MQ=255)
                                                                  cGCCCGGACAGGCCCCCGAAAGCGCCCCCTACCCGATATCACCACAGTGCATCAGCTCATGTCCGCTGCCGACAAACGACAGCAGGGCGACTTTTCCTGc                             <  2:754082/100‑1 (MQ=255)
                                                                  cGCCCGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGc                             >  2:4871886/1‑100 (MQ=255)
                                                                    cccGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCaa                           >  1:463661/1‑100 (MQ=255)
                                                                     ccGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAAc                          <  1:284672/100‑1 (MQ=255)
                                                                     ccGGACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAAc                          <  1:3436326/100‑1 (MQ=255)
                                                                       ggACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAg                        <  1:5295608/100‑1 (MQ=255)
                                                                        gACAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAgg                       <  1:2608312/100‑1 (MQ=255)
                                                                          cAGGCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAgggg                     <  1:1131125/100‑1 (MQ=255)
                                                                            ggCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCa                   >  1:1018962/1‑100 (MQ=255)
                                                                            ggCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCa                   >  2:4685859/1‑100 (MQ=255)
                                                                            ggCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCa                   >  1:4783078/1‑100 (MQ=255)
                                                                            ggCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCa                   >  2:3460702/1‑100 (MQ=255)
                                                                            ggCCCCCGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCa                   <  2:1022838/100‑1 (MQ=255)
                                                                              cccccGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATc                 <  2:40592/100‑1 (MQ=255)
                                                                              cccccGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATc                 >  1:3651420/1‑100 (MQ=255)
                                                                              cccccGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATc                 <  1:5357651/100‑1 (MQ=255)
                                                                               ccccGAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCa                <  1:4941125/100‑1 (MQ=255)
                                                                                cccGAAAGCATCCCCGGACCGATATAACCACAGTGCATCGGCTCGGGTCCGCTGCCGCCACCCGACACCAGGGCGACTTGTCCTTTAACAGGGTCAACAg               >  2:4154436/1‑100 (MQ=255)
                                                                                   gAAAGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCg           >  2:152235/1‑100 (MQ=255)
                                                                                     aaGCCTCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCgg          <  2:1002403/100‑1 (MQ=255)
                                                                                      aGCATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCggg         <  2:2628500/100‑1 (MQ=255)
                                                                                        cATCCCCTGCCCGATCTAACCACAGTGCATCGGTTCGTGTCCGCTGCCGCCCCCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTg       <  2:442402/100‑1 (MQ=255)
                                                                                        cATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTg       <  2:604072/100‑1 (MQ=255)
                                                                                        cATCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTg       <  2:1919982/100‑1 (MQ=255)
                                                                                         aTCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGCGCCTCAGCTCGGGTGa      >  2:5205847/1‑100 (MQ=255)
                                                                                          tCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGAc     <  1:332865/100‑1 (MQ=255)
                                                                                          tCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGAc     <  1:1861082/100‑1 (MQ=255)
                                                                                          tCCCCTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGAc     <  1:875315/100‑1 (MQ=255)
                                                                                            cccTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACAt   >  2:2007696/1‑100 (MQ=255)
                                                                                             ccTGCCCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATa  >  1:3070480/1‑100 (MQ=255)
                                                                                                 |                                                                                               
CTTCGCCGCCATCAACTTGCATGGCGCATTCAAAGATTTTATCGGGCGTCGGTGAGGTGAAAATTTCGCCCGGACAGGCCCCCGAAAGCATCCCCTGACCGATATAACCACAGTGCATCGGCTCGTGTCCGCTGCCGCCACCCGACAGCAGGGCGACTTTTCCTGCAACAGGGGCATCAGCTCGGGTGACATA  >  NC_000913/1250547‑1250739

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: