Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,270,146 T→C 15.1% intergenic (+19/‑92) 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,1460TC15.1% 494.0 / 73.5 205intergenic (+19/‑92)aphA/yjbQacid phosphatase/phosphotransferase, class B, non‑specific/thiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (81/93);  new base C (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 = 8.40e-01

Only 100 of 208 total aligned reads displayed.

CGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGT  >  NC_000913/4270047‑4270244
                                                                                                   |                                                                                                  
cgcgcCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGAt                                                                                                    >  2:1185475/1‑100 (MQ=255)
cgcgcCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGAt                                                                                                    <  2:311049/100‑1 (MQ=255)
 gcgcCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATg                                                                                                   <  2:1467771/100‑1 (MQ=255)
  cgcCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCTTCAATTCAGAATACTGACAGAGCGCGAGAGCGTGATGc                                                                                                  <  1:2731112/100‑1 (MQ=255)
    ccTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGctc                                                                                                >  2:3194681/1‑100 (MQ=255)
      tCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCtctc                                                                                              <  1:1029369/100‑1 (MQ=255)
        cAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATCCTGACAGAGCGGGAGAGCGTGATGctctgc                                                                                             >  2:2938789/1‑97 (MQ=255)
        cAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTccc                                                                                            >  2:3228346/1‑100 (MQ=255)
        cAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTccc                                                                                            >  2:3628136/1‑100 (MQ=255)
        cAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTccc                                                                                            <  1:3444392/100‑1 (MQ=255)
         aaCTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCg                                                                                           <  2:4291996/100‑1 (MQ=255)
         aaCTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCg                                                                                           <  1:1587267/100‑1 (MQ=255)
         aaCTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCg                                                                                           <  2:723068/100‑1 (MQ=255)
         aaCTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCg                                                                                           <  1:3329668/100‑1 (MQ=255)
         aaCTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCg                                                                                           <  2:1257722/100‑1 (MQ=255)
         aaCTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCg                                                                                           <  1:3798924/100‑1 (MQ=255)
         aaCTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCg                                                                                           <  2:3717536/100‑1 (MQ=255)
          aCTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGa                                                                                          >  1:4319267/1‑100 (MQ=255)
          aCTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCgc                                                                                          <  1:1435845/100‑2 (MQ=255)
            tctACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAAt                                                                                        <  2:214523/100‑1 (MQ=255)
            tctACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAAt                                                                                        <  2:3438837/100‑1 (MQ=255)
            tctACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCgctc                                                                                        <  1:2833710/100‑4 (MQ=38)
            tctACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCgctc                                                                                        <  2:2065145/100‑4 (MQ=38)
             ctacctacATACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGCCAGAGCGGGAGAGCGTGATGCTCTCCCGAAt                                                                                        <  2:564111/99‑1 (MQ=255)
             ctacctacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATg                                                                                       >  1:235941/1‑100 (MQ=255)
                cctacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTg                                                                                    <  2:3142042/100‑1 (MQ=255)
                  tacAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCATCACGCTCTCCCgctctgtcag                                                                                  <  2:2452841/100‑10 (MQ=37)
                   acAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGttt                                                                                 <  2:1955909/100‑1 (MQ=255)
                   acAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCGCTCCCGAATGCTGttt                                                                                 <  1:4201376/100‑1 (MQ=255)
                    cAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGtttt                                                                                >  2:3081083/1‑100 (MQ=255)
                      aaCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGtttttt                                                                              >  2:1504237/1‑100 (MQ=255)
                      aaCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGtttttt                                                                              <  2:2292193/100‑1 (MQ=255)
                      aaCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGtttttt                                                                              <  2:808500/100‑1 (MQ=255)
                      aaCCCTTCCCATAAGCGTGTGCGTTTGGGTAAGCGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGATCTCCCGAATGCTGtttttt                                                                              <  1:997681/100‑1 (MQ=255)
                       aCCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGttttttt                                                                             <  2:1826342/100‑1 (MQ=255)
                         ccTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCGATTCAGAATACTGACAGGGCGGGAGAGCGTGCTGCTCTCCCGAATGCTGTTTTTTTaa                                                                           >  2:3273192/1‑100 (MQ=255)
                         ccTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTaa                                                                           <  1:4084384/100‑1 (MQ=255)
                         ccTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTaa                                                                           <  1:2232798/100‑1 (MQ=255)
                          cTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAAt                                                                          <  2:2824224/100‑1 (MQ=255)
                          cTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAAt                                                                          <  2:2728926/100‑1 (MQ=255)
                          cTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAAt                                                                          >  2:156326/1‑100 (MQ=255)
                           ttGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATc                                                                         <  1:2916684/100‑1 (MQ=255)
                             gCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATcac                                                                       <  2:1114776/100‑1 (MQ=255)
                               cacaAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCacac                                                                     >  2:463558/1‑100 (MQ=255)
                                acaAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACAcc                                                                    <  2:2482294/100‑1 (MQ=255)
                                acaAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACAcc                                                                    <  1:3236015/100‑1 (MQ=255)
                                 caAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGtttttttttt                                                                          >  1:2823155/1‑90 (MQ=255)
                                 caAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACAcc                                                                    >  1:2846373/1‑99 (MQ=255)
                                   aGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGGGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCttt                                                                 >  2:199097/1‑100 (MQ=255)
                                     cGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTAt                                                               >  2:3469180/1‑100 (MQ=255)
                                      gcgTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATc                                                              >  2:4096320/3‑100 (MQ=255)
                                            gTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCg                                                        >  1:2042643/1‑100 (MQ=255)
                                            gTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCg                                                        >  2:4139412/1‑100 (MQ=255)
                                              ttGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCt                                                      >  2:2132436/1‑100 (MQ=255)
                                                 gTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTc                                                   >  2:2620806/1‑100 (MQ=255)
                                                  tGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATACTGTTATTTTAATCACACCTTTATCCTTTCGCTGTCt                                                  >  1:2971381/1‑100 (MQ=255)
                                                   gAAGAGGTGATCGTCAATTCAGAATACTGACTGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTc                                                   <  1:1260292/98‑1 (MQ=255)
                                                   gAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCGCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCtt                                                 <  2:128584/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:1523973/100‑1 (MQ=255)
                                                        ggTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgctgc                                            >  1:696535/1‑100 (MQ=255)
                                                        ggTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgctgc                                            <  1:4032113/100‑1 (MQ=255)
                                                        ggTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTtgctg                                             <  1:3624044/99‑1 (MQ=255)
                                                          tGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCaa                                          >  2:498014/1‑100 (MQ=255)
                                                           gATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCaaa                                         >  2:2926608/1‑100 (MQ=255)
                                                           gATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCaaa                                         >  2:4163401/1‑100 (MQ=255)
                                                              cGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAAc                                        <  2:3352698/98‑1 (MQ=255)
                                                              cGTCAAGTCAGAATACTGACAGAGCGGGAGAGCGTTATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACt                                       <  1:183265/99‑1 (MQ=255)
                                                               gTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGa                                     <  2:2529459/100‑1 (MQ=255)
                                                               gTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGa                                     >  1:3414747/1‑100 (MQ=255)
                                                                   aTTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTaa                                 >  1:2456134/1‑100 (MQ=255)
                                                                    ttCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAag                                >  2:4019968/1‑100 (MQ=255)
                                                                     tCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCCGACCTTTATCCTTTCGTTGTCCTTCTGCTAACTGATTCag                                >  2:3579062/1‑99 (MQ=255)
                                                                       aGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAagaga                             >  1:1648647/1‑100 (MQ=255)
                                                                        gAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAagagag                            >  1:784454/1‑100 (MQ=255)
                                                                           tACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGttt                         >  1:2575706/1‑100 (MQ=255)
                                                                            aCTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGtttt                        <  1:1298827/100‑1 (MQ=255)
                                                                              tctgcattcGGGAGAGCATCACGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTAt                      <  1:454926/92‑1 (MQ=38)
                                                                               gACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATc                     >  2:2347376/1‑100 (MQ=255)
                                                                               gACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATc                     <  2:1515396/100‑1 (MQ=255)
                                                                                 cAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCATACCTTTATCCTTTCGCTGTCTTACTTCAAACTGATTAAGAGAGTTTTATCaa                   >  1:3310255/1‑100 (MQ=255)
                                                                                 cAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCaa                   >  1:1147870/1‑100 (MQ=255)
                                                                                  agagCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAg                  >  2:1237530/1‑100 (MQ=255)
                                                                                  agagCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAg                  <  1:1197701/100‑1 (MQ=255)
                                                                                   gagCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAgg                 <  2:3839303/100‑1 (MQ=255)
                                                                                   gagCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAgg                 <  1:934164/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)
                                                                                     gCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGag               <  1:3080317/100‑1 (MQ=255)
                                                                                      cGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGagc              <  2:1224211/100‑1 (MQ=255)
                                                                                        ggAGAGCATCACGCTCTCCCGCATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGagcag            >  1:2337489/1‑100 (MQ=255)
                                                                                        ggAGAGCATCACGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCCCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAgta                >  1:3119368/1‑94 (MQ=255)
                                                                                          agagCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAgcagca          <  2:4273801/100‑1 (MQ=255)
                                                                                             gCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAACGAGCAGcaca        >  1:399935/1‑99 (MQ=255)
                                                                                               gTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACCtgt     >  2:3064647/1‑100 (MQ=255)
                                                                                                tGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACAtgtg    <  2:173103/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
CGCGCCTCCAACTCTACCTACAAACCCTTGCCACAAGCGGGTGCGTTTGGTGAAGAGGTGATCGTCAATTCAGAATACTGACAGAGCGGGAGAGCGTGATGCTCTCCCGAATGCTGTTTTTTTAATCACACCTTTATCCTTTCGCTGTCTTGCTGCAAACTGATTAAGAGAGTTTTATCAAGGAGCAGCACATGTGGT  >  NC_000913/4270047‑4270244

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