Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,738,928 A→C 5.3% intergenic (+23/+20) pheA → / ← tyrA chorismate mutase and prephenate dehydratase, P‑protein/fused chorismate mutase T/prephenate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,738,9280AC 5.3% 247.1 / 3.2 86intergenic (+23/+20)pheA/tyrAchorismate mutase and prephenate dehydratase, P‑protein/fused chorismate mutase T/prephenate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (40/41);  new base C (2/3);  total (42/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

ATTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTGCT  >  NC_000913/2738780‑2739075
                                                                                                                                                    |                                                                                                                                                   
atTCAGGCAAATCTTGAATCCGCGGAAATGCAAAAAGCATTGTAAGAGTTAGGGGCCATCCCCCGTTCAATGCAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGaa                                                                                                                                                    <  2:811901/150‑1 (MQ=255)
   cAGGGCAATCTTGAATCATCGGAAATGCTAAAAGCATTGAAAGCGTTTGGGGAAATCACCCGTTCAATGAAGGTAGTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTTATCCAACCTGATGAAAAGGTGCCGGATGATGTGAatca                                                                                                                                                 <  1:826971/150‑1 (MQ=255)
    aGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAatcat                                                                                                                                                >  1:1152635/1‑150 (MQ=255)
     ggCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAGGGTATTGGGCTGTTAACCAAGTGACACCGTAGTGCCTGTGGATCCACCATGATGAAAAGGTGCCGGATGTTGTGCatcatc                                                                                                                                               >  2:86391/1‑150 (MQ=255)
       ccAATCTTGAATCAGCGGAAATGCAAATATCATTGAAAGAGTTTGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGACCTTAGTGCCTGGTGATCCATCCTGATTAAAAGGTGCCGGATGATTCACATCATCCg                                                                                                                                             <  1:707685/150‑1 (MQ=255)
        cAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCgg                                                                                                                                            <  2:905759/150‑1 (MQ=255)
          aTCTTTAATCAGCTGAAATGAAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATTCACATCATCCGGCa                                                                                                                                          <  2:1263320/150‑1 (MQ=255)
           tCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATTCACATCATCCGGCAc                                                                                                                                         <  1:526896/150‑1 (MQ=255)
                    aGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGattatt                                                                                                                                >  2:947333/1‑150 (MQ=255)
                      cGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGTGCTGTTACACAAGTGAGAACGGAGTGCCTGTTGATCCAACCTGATGTAAAGGTGCCGGATGATGTTAATCAACCGGCACTGGATTATTAc                                                                                                                              >  2:78454/1‑150 (MQ=255)
                       ggAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCGGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACt                                                                                                                             <  1:1085472/150‑1 (MQ=255)
                        gAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTg                                                                                                                            >  1:1037119/1‑150 (MQ=255)
                        gAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTg                                                                                                                            <  2:781414/150‑1 (MQ=255)
                        gAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTCACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTg                                                                                                                            <  1:41742/150‑1 (MQ=255)
                         aaaTGCAAAACGAATTTAAAGAGTTCGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTgg                                                                                                                           <  2:302645/150‑1 (MQ=255)
                          aaTGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGTATGATGTGAATCATCCGGCACTGGATTATTACTGGc                                                                                                                          <  1:966202/150‑1 (MQ=255)
                              cAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGAtt                                                                                                                      >  1:398450/1‑150 (MQ=255)
                              cAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGAtt                                                                                                                      >  1:1240574/1‑150 (MQ=255)
                               aaaaaGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCCCTGGATTATTACTGGCGATTg                                                                                                                     >  2:1145272/1‑150 (MQ=255)
                                aaaaGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATACTCCGGCACTGGATTATTACTGGCGATTGt                                                                                                                    >  1:208126/1‑150 (MQ=255)
                                    gCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCAtt                                                                                                                <  2:606008/150‑1 (MQ=255)
                                    gCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACGGGATTATTACTGGCGATTGTCAtt                                                                                                                >  2:1202846/1‑150 (MQ=255)
                                     cATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTc                                                                                                               <  2:862090/150‑1 (MQ=255)
                                         gAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGTTGCCGGATGATGTGAATCATCCTGCACTGGATTATTTCTGGCGATTGTCATTCGCCt                                                                                                           >  1:223303/1‑150 (MQ=255)
                                         gAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGAACTGGATTATTACTGGCGATTGTCATTCGCCt                                                                                                           >  1:185843/1‑150 (MQ=255)
                                           aaGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCTACCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGa                                                                                                         >  1:649862/1‑150 (MQ=255)
                                             gagTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACg                                                                                                       <  1:356262/150‑1 (MQ=255)
                                             gagTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACg                                                                                                       <  1:723411/150‑1 (MQ=255)
                                              agTTAGGGGAAATCACCCGTTCAATGAAGGAATTGGGCTGTAACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAAGCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGc                                                                                                      <  1:1314041/150‑1 (MQ=255)
                                               gTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCa                                                                                                     <  1:210171/150‑1 (MQ=255)
                                                 tAGGGGAAATCACCCGTTAAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATACGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAAt                                                                                                   >  1:927572/1‑150 (MQ=255)
                                                  tggggAAATAACCCGTTCAATGAAGGTATTTGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGCTCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATa                                                                                                  <  1:48286/149‑1 (MQ=255)
                                                  cggggAACTCTCCCGTTCACTGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAAGCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATa                                                                                                  <  2:660369/149‑1 (MQ=255)
                                                  aGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATa                                                                                                  <  2:548259/150‑1 (MQ=255)
                                                  aGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATa                                                                                                  >  2:925629/1‑150 (MQ=255)
                                                  aGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATa                                                                                                  <  1:919350/150‑1 (MQ=255)
                                                    gggAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAac                                                                                                <  2:5317/150‑1 (MQ=255)
                                                    gggAAAGCACCCGTTCAATGACGGTATTGGGCTTTTACCCAAGTGAGAACGTTGTGCCTGTTGATCCAACCTGATTAAAAGGGGCCGGATGATGTGATTCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAac                                                                                                <  1:1263472/150‑1 (MQ=255)
                                                          tCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTTATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCCTCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCgg                                                                                          >  2:187552/1‑150 (MQ=255)
                                                          tCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCgg                                                                                          <  1:103003/150‑1 (MQ=255)
                                                           cACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGc                                                                                         >  1:1100972/1‑150 (MQ=255)
                                                            aCCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTTTCATTAGCCTGACGCAATAACACGCGGCt                                                                                        >  1:663782/1‑150 (MQ=255)
                                                             cccGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCtt                                                                                       >  1:642493/1‑150 (MQ=255)
                                                             cccGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCtt                                                                                       >  2:797165/1‑150 (MQ=255)
                                                             cccGTTCAATGAAGGTATTGGGCTTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCTCCTGACGCAATAACACGCGGCttt                                                                                      >  1:666758/1‑150 (MQ=255)
                                                              ccGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCttt                                                                                      <  1:51000/150‑1 (MQ=255)
                                                                   cAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATTAAAAGGTGCCGGATGATGTGAATCATACGGCACTGGATTATTACTGGCGATTGTCATTCGCATGACGCAATAACACACTGCTTAAActc                                                                                 >  2:633247/1‑150 (MQ=255)
                                                                   cAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCActc                                                                                 >  2:1043879/1‑150 (MQ=255)
                                                                    aaTGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGAACTGGATTATTAATGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCActat                                                                                >  1:671944/1‑148 (MQ=255)
                                                                      tGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGa                                                                              >  2:926843/1‑150 (MQ=255)
                                                                      tGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGa                                                                              >  2:26866/1‑150 (MQ=255)
                                                                            tATTGGGCTGTTACCCAAGTTAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGc                                                                        <  2:1094571/150‑1 (MQ=255)
                                                                                gggCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCACCCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCCCTGGATTAATACTGGCGATTGTCATTCGCCTGACGCAGTAACACGCGGCTATCACTCTGAAAACGCtgtg                                                                    >  1:596691/1‑150 (MQ=255)
                                                                                    tGTTACCCAAGTGTGAACGTAGTGCCTGTTGATCCAACCTGATCAACAGGTGCCGGAGGAGGTGACTCATCCGGCACGGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTa                                                                <  1:340251/150‑1 (MQ=255)
                                                                                           cAATTGAGAACGTAGTGCCTTTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCg                                                         <  2:1156482/150‑1 (MQ=255)
                                                                                           cAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCg                                                         <  1:455436/150‑1 (MQ=255)
                                                                                           cAAGTGAGAACGACGTGCCTGTTGACCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTGCTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCg                                                         <  2:813433/150‑1 (MQ=255)
                                                                                                    aCGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGc                                                <  2:1286889/150‑1 (MQ=255)
                                                                                                       tAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGCACCAGTGCTcc                                             >  1:329739/1‑150 (MQ=255)
                                                                                                       tAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTcc                                             >  2:182114/1‑150 (MQ=255)
                                                                                                             cTGTTGATACAACCTGATGAAAAGGATCCGGATGATGAGAATCATCCGGCACTGGATTATAACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTATCACTCTGAAAACGCAGTGCGTAATCGCAGAACCAGAGCTACACATTg                                       >  1:49468/1‑150 (MQ=255)
                                                                                                              tGTTGATCCAACCTGATGAAAAGGGGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGGCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGc                                      <  2:632447/150‑1 (MQ=255)
                                                                                                              ggTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGc                                      <  2:628460/149‑1 (MQ=255)
                                                                                                                    aacAACCTGATGAAAAGGTGCCTGATGCTGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAAc                                <  2:799467/148‑1 (MQ=255)
                                                                                                                     ccAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACt                               <  1:761126/150‑1 (MQ=255)
                                                                                                                      cAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTg                              >  1:328573/1‑150 (MQ=255)
                                                                                                                      cAACCTGATGAAAAGGTGCCGGATGATGAGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTg                              <  1:1048982/150‑1 (MQ=255)
                                                                                                                       aaCCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTGTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGt                             >  2:510648/1‑150 (MQ=255)
                                                                                                                        aCCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTc                            <  1:685523/150‑1 (MQ=255)
                                                                                                                         ccTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCa                           >  1:914106/1‑150 (MQ=255)
                                                                                                                           tgatgaAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAAt                         >  1:17191/1‑150 (MQ=255)
                                                                                                                            gatgaAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATa                        <  2:1275605/150‑1 (MQ=255)
                                                                                                                            gatgaAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATa                        >  2:435313/1‑150 (MQ=255)
                                                                                                                               gaAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAAc                     <  1:699968/150‑1 (MQ=255)
                                                                                                                               gaAAAGGTGCCGGATGATGTGAATCATCCGGCACGGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAAc                     <  2:223642/150‑1 (MQ=255)
                                                                                                                               gaAAAGGTGCCGGATGAGGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAAc                     <  2:118488/150‑1 (MQ=255)
                                                                                                                                aaaaGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACg                    <  1:968971/150‑1 (MQ=255)
                                                                                                                                 tccagTGCCGGATGATTCACATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGc                   >  2:922511/5‑150 (MQ=255)
                                                                                                                                 ctaGGTTCCGTATGACGTGAATCAACCGTCTCGTGAATATTCCTTGAGATTGTCATTCTCCTGACGCAATAACAAGAGGCTTTCACTCTTAACAATCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTTCGGAAACTGGAAATAAACGc                   <  1:807188/148‑1 (MQ=255)
                                                                                                                                  taGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGcc                  <  1:1086208/149‑1 (MQ=255)
                                                                                                                                      tGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCt              >  1:892800/1‑150 (MQ=255)
                                                                                                                                      ggCCGGATGATGTTAATCATCCGGCACGGGATTATTACTGGCGAGTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTGGCGGAAACTGTCAATAAACGCCTGCt              <  1:66641/149‑1 (MQ=255)
                                                                                                                                        ccGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTa            >  2:667998/1‑150 (MQ=255)
                                                                                                                                             tgatgTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGcc       >  2:600636/1‑150 (MQ=255)
                                                                                                                                                tgtgAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTg    >  2:793843/1‑150 (MQ=255)
                                                                                                                                                 gtgAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTTTCGCCCTGc   >  2:1209074/1‑150 (MQ=255)
                                                                                                                                                  tgAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTGCt  >  1:816232/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTGCT  >  NC_000913/2738780‑2739075

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