Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,260,901 C→T 16.1% intergenic (‑98/+27) dauA ← / ← prs C4‑dicarboxylic acid transporter/phosphoribosylpyrophosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,260,9010CT16.1% 330.0 / 45.9 136intergenic (‑98/+27)dauA/prsC4‑dicarboxylic acid transporter/phosphoribosylpyrophosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (50/64);  new base T (11/11);  total (61/75)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

Only 100 of 151 total aligned reads displayed.

ACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGAC  >  NC_000913/1260802‑1261000
                                                                                                   |                                                                                                   
acaATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAATGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACccgcc                                                                                                     <  2:2844681/100‑1 (MQ=255)
acaATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACccgcc                                                                                                     <  1:1883182/100‑1 (MQ=255)
 caATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACccgccg                                                                                                    >  2:2301623/1‑100 (MQ=255)
   ataatGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATGCAAATGGATAGTACAGACAAGAAAACCcgccgc                                                                                                   >  2:2501176/1‑99 (MQ=255)
     aatGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGc                                                                                                <  2:303967/100‑1 (MQ=255)
        gtgtgaTCCGTGAGGGCATCATCTGTCAACTAAAAGGCAAGAAAAGGAGGCATAGGTCATACATATGTATAAAAAAAAAAAAAAAACCCGCCGCAGCgg                                                                                              <  2:1263654/94‑1 (MQ=255)
        gTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCggg                                                                                             <  2:1662295/100‑1 (MQ=255)
        gTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCggg                                                                                             <  2:47325/100‑1 (MQ=255)
         tttCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGt                                                                                            <  1:1527902/100‑1 (MQ=255)
           tCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTc                                                                                           >  2:1142440/1‑99 (MQ=255)
           tCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCt                                                                                          <  1:2872617/100‑1 (MQ=255)
             aTCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAGACCCGCTGCGGCGGGTtttt                                                                                        <  1:1981765/100‑1 (MQ=255)
              tCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCTGCGGCGGGTttttt                                                                                       <  1:3395818/100‑2 (MQ=255)
              tCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTg                                                                                       <  1:688336/100‑1 (MQ=255)
              tCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTg                                                                                       <  2:3791032/100‑1 (MQ=255)
              tCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTg                                                                                       <  1:2074501/100‑1 (MQ=255)
               ccGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCTGCAGCGGGTtttttt                                                                                      <  1:2761494/100‑3 (MQ=255)
               ccGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGa                                                                                      <  1:309075/100‑1 (MQ=255)
                cGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAg                                                                                     >  1:2393324/1‑100 (MQ=255)
                cGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAg                                                                                     >  1:3311630/1‑100 (MQ=255)
                    aGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTCCAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCggg                                                                                 >  2:1095925/1‑100 (MQ=255)
                    aGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCgg                                                                                 >  1:2143549/1‑100 (MQ=255)
                      cgcATTATCTGACAACTAAATTGCAGGTCAAGGAGGCATAAGTCATGGAAATGGATATTACAGACAAAAAGACCCGCTGCGGCGGGTttttttgtctgta                                                                               <  2:2789762/100‑10 (MQ=16)
                      cgcATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCgg                                                                                 >  2:1985696/1‑98 (MQ=255)
                      cgcATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGt                                                                               <  2:534916/100‑1 (MQ=255)
                       gcATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGtt                                                                              >  1:905559/1‑100 (MQ=255)
                        catcatCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGAAAAAAAAAAAAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTc                                                                             <  2:1090554/100‑1 (MQ=255)
                         atcatcTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCg                                                                            <  2:3380261/100‑1 (MQ=255)
                            atcTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGAtt                                                                         <  1:2330498/100‑1 (MQ=255)
                             tcTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTa                                                                        >  1:3517689/1‑100 (MQ=255)
                              cTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAg                                                                       <  1:333596/100‑1 (MQ=255)
                              cTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAg                                                                       <  1:1862409/100‑1 (MQ=255)
                              cTGCCAACTAAATGGCAGGTGAAGGAGGCCTAGGTCATACAAATGGATATTACAGACAAAAAAACCTGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAg                                                                       <  2:2518875/100‑1 (MQ=255)
                                gTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAgtg                                                                     <  2:628714/100‑1 (MQ=255)
                                  cAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTGCAGACAAAAAAACCCGCCGCAGAGGATATTTGAGCCGGGATCGATTAGTGtt                                                                   >  2:3524332/1‑100 (MQ=255)
                                  cAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGtt                                                                   >  2:1002705/1‑100 (MQ=255)
                                  cAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGtt                                                                   >  2:296277/1‑100 (MQ=255)
                                   aaCTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTc                                                                  <  2:2413030/100‑1 (MQ=255)
                                     cTAAATGGCAGGTGGAGGAGGCATAGCTCATACAAATGGATATAACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCgc                                                                <  2:830289/100‑2 (MQ=255)
                                     cTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGa                                                                <  1:4102974/100‑1 (MQ=255)
                                       aaaTGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAAc                                                              <  2:392819/100‑1 (MQ=255)
                                       aaaTGGCAGGTGAAGGAGGCATAGGGCATAAAAATGGGTATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGa                                                                <  1:1095925/98‑1 (MQ=255)
                                         aTGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACAt                                                            <  2:3934875/100‑1 (MQ=255)
                                           ggCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATgg                                                          <  1:1825720/100‑1 (MQ=255)
                                            gCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATgg                                                          >  1:2943320/1‑99 (MQ=255)
                                            gCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGc                                                         <  2:2143546/100‑1 (MQ=255)
                                               ggTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCaga                                                      >  2:3732180/1‑100 (MQ=255)
                                                     ggaggCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGa                                                >  2:3008150/1‑100 (MQ=255)
                                                     ggaggCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGa                                                <  2:3026029/100‑1 (MQ=255)
                                                     ggaggCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGa                                                <  1:2476797/100‑1 (MQ=255)
                                                        ggCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAttc                                             >  2:2684386/1‑100 (MQ=255)
                                                         gCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGAttct                                            <  2:4011954/100‑1 (MQ=255)
                                                              ggTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGtt                                       <  2:1292401/100‑1 (MQ=255)
                                                              ggTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGtt                                       >  2:1926504/1‑100 (MQ=255)
                                                              ggTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGtt                                       >  1:307984/1‑100 (MQ=255)
                                                                tCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGc                                     <  1:1297500/100‑1 (MQ=255)
                                                                 cATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTTGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCt                                    <  1:2956007/100‑1 (MQ=255)
                                                                  aTACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTg                                   <  2:3270130/100‑1 (MQ=255)
                                                                   tACAAATGGATATTACAGACAAAACACCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGc                                     <  2:2997818/97‑1 (MQ=255)
                                                                   tACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGa                                  >  2:2765133/1‑100 (MQ=255)
                                                                     cAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATa                                >  1:2222659/1‑100 (MQ=255)
                                                                     cAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATa                                >  1:242419/1‑100 (MQ=255)
                                                                       aaTGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATacg                              >  1:2041837/1‑100 (MQ=255)
                                                                        aTGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATacga                             <  1:1000631/100‑1 (MQ=255)
                                                                         tGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATacgac                            <  2:1290046/100‑1 (MQ=255)
                                                                         tGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATacgac                            <  1:3582239/100‑1 (MQ=255)
                                                                         tGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATacgac                            <  1:3642510/100‑1 (MQ=255)
                                                                          ggATATTACAGACAAAAAAACCCGCCGCAGCGGGGCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATacgacg                           >  2:617784/1‑100 (MQ=255)
                                                                           gATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATAcgacga                          >  1:3896679/1‑100 (MQ=255)
                                                                           gATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATAcgacga                          >  2:12203/1‑100 (MQ=255)
                                                                            atatTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGaa                         <  1:3281422/100‑1 (MQ=255)
                                                                               ttACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCg                      >  2:1686626/1‑100 (MQ=255)
                                                                               ttACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCg                      <  1:1719042/100‑1 (MQ=255)
                                                                                   aGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGCTTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCtt                   <  2:2691759/99‑1 (MQ=255)
                                                                                   aGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTc                  <  1:50550/100‑1 (MQ=255)
                                                                                   aGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTc                  <  1:1383847/100‑1 (MQ=255)
                                                                                    gACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCg                 >  1:2533548/1‑100 (MQ=255)
                                                                                     gctcaaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCgg                >  1:368426/9‑100 (MQ=38)
                                                                                     aCAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCgg                <  1:1726147/100‑1 (MQ=255)
                                                                                      ctcaaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGc               >  2:3813377/8‑100 (MQ=38)
                                                                                      ctcaaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAAATCGATTATTCGTTGCTGATACGACGAATCGCTTCGGc               >  2:889539/4‑100 (MQ=37)
                                                                                      cAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGTTTAGTGTTCGAACATGGTAGATATCGATTCTTCGTTGCTGATACGACGAATCGCTTCgg                >  2:2951888/1‑99 (MQ=255)
                                                                                      cAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGc               <  2:2413866/100‑1 (MQ=255)
                                                                                      cAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATCAGTGTTCGAACATGGCAGAGATCGATTCTTTGTTGCTGATACGACGAATCGCTTCGGc               >  2:3619464/1‑100 (MQ=255)
                                                                                       aaaaaaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGcc              <  1:3910187/100‑1 (MQ=255)
                                                                                        caaagaCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCa             >  1:1571037/2‑100 (MQ=255)
                                                                                        aaaaaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCa             <  1:783617/100‑1 (MQ=255)
                                                                                        aaaaaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCa             >  1:216369/1‑100 (MQ=255)
                                                                                        aaaaaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCa             >  1:2458761/1‑100 (MQ=255)
                                                                                        aaaaaaCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCa             >  1:2572184/1‑100 (MQ=255)
                                                                                         aaagaCCCGCTGCGGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAg            >  1:3670230/1‑100 (MQ=255)
                                                                                          aagaCCCGCTGCGGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGc           >  1:1047065/4‑100 (MQ=255)
                                                                                             aCCCGCTGCGGCGGGTTTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATa        >  1:4137063/1‑100 (MQ=255)
                                                                                              cccGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATa        <  2:3316131/99‑1 (MQ=255)
                                                                                               ccgccgCAGCGTGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATAcc      <  1:1392074/100‑1 (MQ=255)
                                                                                               ccgccgCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATAcc      >  2:1258494/1‑100 (MQ=255)
                                                                                               ccgccgCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATAcc      <  2:2817437/100‑1 (MQ=255)
                                                                                                 gccgcAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTg    <  1:3732201/100‑1 (MQ=255)
                                                                                                 gccgcAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTg    <  2:735627/100‑1 (MQ=255)
                                                                                                   cgcAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGAc  >  2:4141120/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATACCTGAC  >  NC_000913/1260802‑1261000

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