Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,270,142 G→A 15.1% intergenic (+15/‑96) aphA → / → yjbQ acid phosphatase/phosphotransferase, class B, non‑specific/thiamine phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,270,1420GA15.1% 523.2 / 68.8 205intergenic (+15/‑96)aphA/yjbQacid phosphatase/phosphotransferase, class B, non‑specific/thiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (81/93);  new base A (11/20);  total (92/113)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.52e-01

Only 100 of 207 total aligned reads displayed.

TCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGT  >  NC_000913/4270043‑4270241
                                                                                                   |                                                                                                   
 cTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGt                                                                                                    <  2:1565855/100‑1 (MQ=255)
 cTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGt                                                                                                    <  2:2310475/100‑1 (MQ=255)
   gcgcgcCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGa                                                                                                  <  1:2934039/100‑1 (MQ=255)
    cgcgcCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGAt                                                                                                 >  2:1197705/1‑100 (MQ=255)
     gcgcCTCCAACTCTACCTACAAACCCTTGCCACAAGCGTGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATg                                                                                                <  1:2358964/100‑1 (MQ=255)
     gcgcCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATg                                                                                                <  2:1467771/100‑1 (MQ=255)
       gcCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGct                                                                                              >  2:3127426/1‑100 (MQ=255)
       gcCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGct                                                                                              >  1:3142024/1‑100 (MQ=255)
          tCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCtctc                                                                                           >  2:1862693/1‑100 (MQ=255)
          tCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCtctc                                                                                           <  1:1029369/100‑1 (MQ=255)
            cAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATCCTGACAGAGCGGGAGAGCGTGATGctctgc                                                                                          >  2:2938789/1‑97 (MQ=255)
            cAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTccc                                                                                         >  2:3628136/1‑100 (MQ=255)
             aaCTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCg                                                                                        <  2:4291996/100‑1 (MQ=255)
             aaCTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCg                                                                                        <  2:723068/100‑1 (MQ=255)
             aaCTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCg                                                                                        <  1:1587267/100‑1 (MQ=255)
             aaCTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCg                                                                                        <  1:3329668/100‑1 (MQ=255)
              aCTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCgc                                                                                       <  2:1357527/100‑2 (MQ=255)
              aCTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCgc                                                                                       <  1:1435845/100‑2 (MQ=255)
                tctACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAAt                                                                                     <  2:214523/100‑1 (MQ=255)
                tctACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAAt                                                                                     >  1:3981512/1‑100 (MQ=255)
                tctACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCgctc                                                                                     <  2:2065145/100‑4 (MQ=38)
                 ctacctacATACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGCCAGAGCGGGAGAGCGTGATGCTCTCCCGAAt                                                                                     <  2:564111/99‑1 (MQ=255)
                  tacctacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGc                                                                                   >  2:3170556/1‑100 (MQ=255)
                      tacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGtt                                                                               >  2:3313868/1‑100 (MQ=255)
                      tacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCgctctgtcag                                                                               <  2:2452841/100‑10 (MQ=37)
                          aaCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGtttttt                                                                           <  2:808500/100‑1 (MQ=255)
                          aaCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGtttttt                                                                           <  1:2736573/100‑1 (MQ=255)
                            cccTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTa                                                                         >  2:3462465/1‑100 (MQ=255)
                             ccTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTaa                                                                        <  1:2232798/100‑1 (MQ=255)
                             ccTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTa                                                                         >  1:424472/1‑99 (MQ=255)
                              cTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAAt                                                                       <  2:3229193/100‑1 (MQ=255)
                                 gCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATcac                                                                    <  2:336073/100‑1 (MQ=255)
                                  ncacaAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATcaca                                                                   >  2:2282197/2‑100 (MQ=255)
                                  ccACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATcaca                                                                   <  1:4335816/100‑1 (MQ=255)
                                   cacaAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCacac                                                                  >  2:463558/1‑100 (MQ=255)
                                    acaAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACAcc                                                                 <  1:3236015/100‑1 (MQ=255)
                                    acaAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACAcc                                                                 <  2:2482294/100‑1 (MQ=255)
                                      aaGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCtt                                                               <  1:2842554/100‑1 (MQ=255)
                                       aGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGGGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCttt                                                              >  2:199097/1‑100 (MQ=255)
                                        gCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTa                                                             <  1:2870354/100‑1 (MQ=255)
                                          gggTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATc                                                           <  1:3685152/100‑1 (MQ=255)
                                          gcgTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATc                                                           >  2:4096320/3‑100 (MQ=255)
                                            gTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCt                                                         <  1:4001274/100‑1 (MQ=255)
                                                gTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCg                                                     >  2:4139412/1‑100 (MQ=255)
                                                gTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCg                                                     <  2:1179529/100‑1 (MQ=255)
                                                gTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCg                                                     >  1:2042643/1‑100 (MQ=255)
                                                  ttGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCt                                                   >  2:2132436/1‑100 (MQ=255)
                                                   tGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAGTGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTg                                                  <  2:3105925/100‑1 (MQ=255)
                                                   tGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTg                                                  <  1:465051/100‑1 (MQ=255)
                                                     gTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTc                                                >  2:2620806/1‑100 (MQ=255)
                                                       gAAGAGGTGATCGTCAATTCAGAATACTGACTGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTc                                                <  1:1260292/98‑1 (MQ=255)
                                                       gAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCtt                                              <  1:3510976/100‑1 (MQ=255)
                                                         agagGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgc                                            >  1:989062/1‑100 (MQ=255)
                                                         agagGTGATCGTCAAGTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgc                                            <  1:846397/100‑1 (MQ=255)
                                                          gagGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgct                                           <  1:1504231/100‑1 (MQ=255)
                                                           agGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgctg                                          <  2:1628407/100‑1 (MQ=255)
                                                            ggTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgctgc                                         >  2:1487699/1‑100 (MQ=255)
                                                            ggTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgctgc                                         >  1:696535/1‑100 (MQ=255)
                                                            ggTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgctgc                                         >  1:4069306/1‑100 (MQ=255)
                                                            ggTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgctgc                                         >  2:1225611/1‑100 (MQ=255)
                                                              tGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCaa                                       >  2:498014/1‑100 (MQ=255)
                                                               gATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCaaa                                      >  1:311862/1‑100 (MQ=255)
                                                               gATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCaaa                                      >  2:4163401/1‑100 (MQ=255)
                                                                aTCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAAc                                     <  2:3582683/100‑1 (MQ=255)
                                                                  cGTCAAGTCAGAATACTGACAGAGCGGGAGAGCGTTATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACt                                    <  1:183265/99‑1 (MQ=255)
                                                                   gTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGa                                  >  1:3414747/1‑100 (MQ=255)
                                                                   gTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGa                                  <  2:2529459/100‑1 (MQ=255)
                                                                      aaTTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTa                               >  2:3227059/1‑100 (MQ=255)
                                                                       aTTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTaa                              >  1:2456134/1‑100 (MQ=255)
                                                                        ttCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAag                             >  2:4019968/1‑100 (MQ=255)
                                                                          cAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAagag                           <  2:3770716/100‑1 (MQ=255)
                                                                           aGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAagaga                          >  1:1648647/1‑100 (MQ=255)
                                                                            gAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAagagag                         >  1:784454/1‑100 (MQ=255)
                                                                               tACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGttt                      >  1:2575706/1‑100 (MQ=255)
                                                                                aCTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGtttt                     >  2:3573/1‑100 (MQ=255)
                                                                                aCTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGtttt                     <  1:1298827/100‑1 (MQ=255)
                                                                                   gACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATc                  >  2:2347376/1‑100 (MQ=255)
                                                                                   gACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATc                  >  2:1519926/1‑100 (MQ=255)
                                                                                   cagcattcGGGAGAGCATCACGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATc                  >  2:1742779/8‑100 (MQ=37)
                                                                                     cAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCATACCTTTATCCTTTCGCTGTCTTACTTCAAACTGATTAAGAGAGTTTTATCaa                >  1:3310255/1‑100 (MQ=255)
                                                                                     cAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCaa                >  1:1147870/1‑100 (MQ=255)
                                                                                     cAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCaa                <  1:2831859/100‑1 (MQ=255)
                                                                                      cattcGGGAGAGCATCACGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAg               >  2:3073416/5‑100 (MQ=38)
                                                                                      agagCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAg               <  1:1197701/100‑1 (MQ=255)
                                                                                       gagCGTGAGAGCGCGCTGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAGCTGATTAAGAGAGTTTTATCAAgg              <  1:299056/100‑1 (MQ=255)
                                                                                       gagCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAgg              <  1:934164/100‑1 (MQ=255)
                                                                                       gagCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAgg              >  2:1246098/1‑100 (MQ=255)
                                                                                       gagCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAgg              <  1:227627/100‑1 (MQ=255)
                                                                                       gagCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAgg              <  2:2830519/100‑1 (MQ=255)
                                                                                       attcGGGAGAGCATCACGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAgg              >  2:4407790/4‑100 (MQ=38)
                                                                                        agCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGa             >  2:2395388/1‑100 (MQ=255)
                                                                                         gCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGag            <  1:1487693/100‑1 (MQ=255)
                                                                                           gggAGAGCATCACGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGagca          >  2:1514723/1‑100 (MQ=255)
                                                                                            ggAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTTTTTTTATTTTTATCCTTTAGCTGTCTTGCTGTAAACTGATTAATAGAGTTTTATCAAGGagca          >  1:4192243/1‑99 (MQ=255)
                                                                                             gagaGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGagcagc        <  2:1527152/100‑1 (MQ=255)
                                                                                              agagCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAgcagca       <  2:3941996/100‑1 (MQ=255)
                                                                                              agagCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAgcagca       <  2:4273801/100‑1 (MQ=255)
                                                                                                 gCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATCAAGAGAGTTTTATCAAGGAGCAGcacct    >  2:126637/1‑98 (MQ=255)
                                                                                                  cGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACAtg   >  2:2775009/1‑100 (MQ=255)
                                                                                                   gTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACCtgt  >  2:3064647/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TCTGCGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGT  >  NC_000913/4270043‑4270241

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