Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,738,927 G→A 10.3% intergenic (+22/+21) 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,9270GA10.3% 277.5 / 14.1 99intergenic (+22/+21)pheA/tyrAchorismate mutase and prephenate dehydratase, P‑protein/fused chorismate mutase T/prephenate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (44/43);  new base A (5/5);  total (50/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.31e-01

Only 100 of 103 total aligned reads displayed.

ATATTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTGCTC  >  NC_000913/2738778‑2739076
                                                                                                                                                     |                                                                                                                                                     
atatTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGAtgtg                                                                                                                                                       <  1:511525/150‑1 (MQ=255)
atatTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGAtgtg                                                                                                                                                       <  1:1459430/150‑1 (MQ=255)
   ttCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAat                                                                                                                                                    >  2:1114326/1‑150 (MQ=255)
     cAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAatca                                                                                                                                                  >  1:179730/1‑150 (MQ=255)
        gCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATcc                                                                                                                                               >  2:1238822/1‑150 (MQ=255)
         ccAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCg                                                                                                                                              >  2:449880/1‑150 (MQ=255)
          cAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATTCACATCATCCgg                                                                                                                                             <  2:1345830/150‑1 (MQ=255)
          cAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCgg                                                                                                                                             <  2:814249/150‑1 (MQ=255)
          cAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCgg                                                                                                                                             <  1:550905/150‑1 (MQ=255)
           aaTCTTGAATCAGCGGAAATTCAAAAATCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATTCACATCATCCGGc                                                                                                                                            <  1:669782/150‑1 (MQ=255)
           aaTCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATTCACATCATCCGGc                                                                                                                                            <  1:819293/150‑1 (MQ=255)
           aaTCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATTCACATCATCCGGc                                                                                                                                            <  2:40997/150‑1 (MQ=255)
            aTCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATTCACATCATCCGGCa                                                                                                                                           <  1:956091/150‑1 (MQ=255)
               ttGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTg                                                                                                                                        <  2:1678668/150‑1 (MQ=255)
                tGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTgg                                                                                                                                       >  1:1286206/1‑150 (MQ=255)
                tGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTAACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTgg                                                                                                                                       >  1:1283362/1‑150 (MQ=255)
                     cAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGCTCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGattat                                                                                                                                  >  2:1075942/1‑150 (MQ=255)
                     cAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGattat                                                                                                                                  >  1:469768/1‑150 (MQ=255)
                       gCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGAttatta                                                                                                                                <  1:2976/150‑1 (MQ=255)
                         tgAAATGCAAAAAGAATTGAAAGAGTTAGGTTAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATTATGTGAATCATCCGGCACTGGATTATTACt                                                                                                                              <  1:1719951/149‑1 (MQ=255)
                         ggAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGCAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGGGAATCATCCGGCACTGGATTATTACt                                                                                                                              <  1:1713440/150‑1 (MQ=255)
                         ggAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTCATGAAAAGGTGCCGGATGATGTGAATCATCCGGAACTGGATTATTACt                                                                                                                              >  1:1346816/1‑150 (MQ=255)
                              tGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGa                                                                                                                         >  2:1041861/1‑150 (MQ=255)
                                cAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGAtt                                                                                                                       >  1:979514/1‑150 (MQ=255)
                                cAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGAtt                                                                                                                       <  1:229437/150‑1 (MQ=255)
                                 aaaaaGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTg                                                                                                                      >  2:507116/1‑150 (MQ=255)
                                  aaaaGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGt                                                                                                                     <  1:449880/150‑1 (MQ=255)
                                      gCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCAtt                                                                                                                 >  2:1307119/1‑150 (MQ=255)
                                       cATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTc                                                                                                                >  2:676098/1‑150 (MQ=255)
                                          tGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGATCGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGcc                                                                                                             >  2:379112/1‑150 (MQ=255)
                                           gAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCt                                                                                                            >  2:156926/1‑150 (MQ=255)
                                             aaGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACAGGCGAATGACAATCGccagt                                                                                                          <  1:1496976/150‑4 (MQ=255)
                                              agagTTAGGGGAAATCACCCGTGCAATGACGGCATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGAAGGACAATCGccagta                                                                                                         <  1:528558/150‑5 (MQ=255)
                                               tagTTAGGGGAAATCACGCGTTCAATGAAGGTTTTGGGCTGTTACCCAAGTGAGAACGTTGTGCCTTTTGATCCAACCTGATTAAATGGTGCCGTATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACg                                                                                                        <  2:576967/149‑1 (MQ=255)
                                               gattatggggAAAGCACACGTTCAATGACTGTTTTGGGCGGTTACCGAAGTGAGAACGTTGTGCCTGTTGTTCCAACCTGACGAATTGGTGCGGGTTGATGTGATTCATCCGGAACTGGATTATTACTGGCGATTGTCATTCGCCTGACg                                                                                                        <  2:71677/144‑1 (MQ=255)
                                               gagTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGAATGACAATCGccagtaa                                                                                                        <  1:1468667/150‑6 (MQ=255)
                                                   tatcggAAATCACCCGTTCAATGAAGGAATTGGGCTGTGACCGAAGTGAGACCGTAGTGCCTGTTGATCCAACCTGATGAAAAGGGGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAAt                                                                                                    <  1:863835/146‑1 (MQ=255)
                                                   tAGGGGAAATCACCCGTTCAATTAAGGTAATGTGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTTATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAAt                                                                                                    <  1:869745/150‑1 (MQ=255)
                                                       ggAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAaca                                                                                                >  1:426319/1‑150 (MQ=255)
                                                       ggAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACATGATGAAAGGGTGCCGGATGATGTGAATCATCCGTCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAacc                                                                                                >  2:1554098/1‑149 (MQ=255)
                                                        gAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAacac                                                                                               >  2:4232/1‑150 (MQ=255)
                                                         aaaTCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACAcg                                                                                              >  1:1622292/1‑150 (MQ=255)
                                                          aaTCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACAcgc                                                                                             <  1:1242881/150‑1 (MQ=255)
                                                             cACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGc                                                                                          >  2:1001657/1‑150 (MQ=255)
                                                             aaCCCGTTCAATGTAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGc                                                                                          <  2:874734/149‑1 (MQ=255)
                                                              aCCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTAGATCCAACCTGATGAAAAGGTGCCGGATGATCTCAATCATCCGGCACTGTATTATAACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCt                                                                                         >  2:1449404/1‑150 (MQ=255)
                                                               cccGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCtt                                                                                        <  1:676098/150‑1 (MQ=255)
                                                               cccGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCtt                                                                                        <  1:229404/150‑1 (MQ=255)
                                                                ccGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCttt                                                                                       <  2:426319/150‑1 (MQ=255)
                                                                 cGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTc                                                                                      <  1:543876/150‑1 (MQ=255)
                                                                  gTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGAAGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCa                                                                                     >  2:1554198/1‑150 (MQ=255)
                                                                   ttCAATTAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAAACTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCAc                                                                                    <  2:1231717/150‑1 (MQ=255)
                                                                    tCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACGGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCAct                                                                                   <  2:159347/150‑1 (MQ=255)
                                                                       aTGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTg                                                                                <  2:744335/150‑1 (MQ=255)
                                                                       aTGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTg                                                                                >  1:28255/1‑150 (MQ=255)
                                                                         gAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGaa                                                                              >  1:1316518/1‑150 (MQ=255)
                                                                          aaGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGaaa                                                                             <  2:624266/150‑1 (MQ=255)
                                                                                tggggCGGTTACCCAGTTGTGAACGTAGGGCCTGTTGATCCAACCTGTTGAAATGGGGCCGGAGGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAACCGCtg                                                                       <  1:849077/148‑1 (MQ=255)
                                                                                  gggCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCtgtg                                                                     <  2:447091/150‑1 (MQ=255)
                                                                                  gggCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCtgtg                                                                     <  1:280499/150‑1 (MQ=255)
                                                                                   ggCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGc                                                                    <  1:4232/150‑1 (MQ=255)
                                                                                          aCCCATGTGAGAACGTTGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGTATGATGTTAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCg                                                             <  2:1526112/150‑1 (MQ=255)
                                                                                           cccAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAACGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGc                                                            <  2:1221906/150‑1 (MQ=255)
                                                                                            ccAAGTGAGAACGTAGTGCCTGTTGATCCAACCTTATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGcc                                                           <  1:1130773/150‑1 (MQ=255)
                                                                                             cAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCg                                                          >  1:171513/1‑150 (MQ=255)
                                                                                             cAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCg                                                          >  2:21958/1‑150 (MQ=255)
                                                                                             cAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCAATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCg                                                          >  1:1346733/1‑150 (MQ=255)
                                                                                               aGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGaa                                                        <  1:112814/150‑1 (MQ=255)
                                                                                                 tGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCCGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAAcc                                                      >  2:662919/1‑150 (MQ=255)
                                                                                                   agaACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGCCGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAg                                                    >  1:1290185/1‑150 (MQ=255)
                                                                                                   agaACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAg                                                    >  1:1283528/1‑150 (MQ=255)
                                                                                                      aCGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGc                                                 <  2:527235/150‑1 (MQ=255)
                                                                                                      aCGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGc                                                 <  2:1159781/150‑1 (MQ=255)
                                                                                                          aGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCa                                             >  2:1466511/1‑150 (MQ=255)
                                                                                                                     aTCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGaaa                                  >  2:1082603/1‑150 (MQ=255)
                                                                                                                     aTCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGaaa                                  >  1:1361534/1‑150 (MQ=255)
                                                                                                                        cAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTg                               <  1:1114326/150‑1 (MQ=255)
                                                                                                                         aaCCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGCCGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGt                              >  2:667513/1‑150 (MQ=255)
                                                                                                                           ccTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCa                            <  1:737335/150‑1 (MQ=255)
                                                                                                                           ccTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCa                            >  2:1098213/1‑150 (MQ=255)
                                                                                                                           ccTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCa                            <  1:1596548/150‑1 (MQ=255)
                                                                                                                            cTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCaa                           <  1:298895/150‑1 (MQ=255)
                                                                                                                            atgatgaAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATAGCCGAACAAGTGCTCCACCTTGCGGAAACTGTCaa                           >  2:173155/2‑150 (MQ=255)
                                                                                                                             tgatgaAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAAt                          <  2:171513/150‑1 (MQ=255)
                                                                                                                                  aaaaGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACg                     >  1:536871/1‑150 (MQ=255)
                                                                                                                                  aaaaGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACg                     <  2:469768/150‑1 (MQ=255)
                                                                                                                                   aaaGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGc                    >  2:530873/1‑150 (MQ=255)
                                                                                                                                      agTGCCGGATGATTCACATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTg                 >  1:1670566/2‑150 (MQ=255)
                                                                                                                                        tGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCt               >  2:466402/1‑150 (MQ=255)
                                                                                                                                         gCCGGATGATTCACATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCtt              >  1:402280/1‑150 (MQ=255)
                                                                                                                                          ccGGATGATTCACATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCCCCTTGCGGAAACTGTCAATAAACGCCTGCTTa             >  2:1038396/1‑150 (MQ=255)
                                                                                                                                          ccGGATGATTCACATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTa             >  2:1569470/1‑150 (MQ=255)
                                                                                                                                          ccGGATGATTCACATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTa             >  2:358911/1‑150 (MQ=255)
                                                                                                                                          ccGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTACACTCTGAAAACGCTGTGCGTAATCGCCGACCCAGTGCTCCACCTTGCGGAACCTGTCAATAAACGCCTGCTTa             >  2:319377/1‑150 (MQ=255)
                                                                                                                                           cGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTAt            >  1:187951/1‑150 (MQ=255)
                                                                                                                                            ggATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATc           <  2:1078938/150‑1 (MQ=255)
                                                                                                                                               tgatgTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGATCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGcc        >  1:249992/1‑150 (MQ=255)
                                                                                                                                                gatgTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGccc       <  1:1349745/150‑1 (MQ=255)
                                                                                                                                                   gtgAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTGc    <  1:388563/150‑1 (MQ=255)
                                                                                                                                                     gAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTGCTc  <  1:156926/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATATTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTGCTC  >  NC_000913/2738778‑2739076

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