Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,738,925 G→T 7.2% intergenic (+20/+23) pheA → / ← tyrA fused chorismate mutase/prephenate dehydratase/fused chorismate mutase/prephenate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,738,9250GT 7.2% 468.6 / 16.0 139intergenic (+20/+23)pheA/tyrAfused chorismate mutase/prephenate dehydratase/fused chorismate mutase/prephenate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (55/74);  new base T (5/5);  total (60/79)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.54e-01

Only 100 of 153 total aligned reads displayed.

GGATATTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTG  >  NC_000913/2738776‑2739073
                                                                                                                                                     |                                                                                                                                                    
ggATATTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGatgatg                                                                                                                                                      >  2:878141/1‑150 (MQ=255)
 gATATTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGAtgt                                                                                                                                                     <  1:2438416/150‑1 (MQ=255)
 gATATTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGAtgt                                                                                                                                                     <  2:2276363/150‑1 (MQ=255)
  atatTCAGGCCATTCTTGAATCAGCGGAAATGCATTAAGCATTGAAAGATTTAGGGGAAATCACGCGTTCAATGAATGTATTGGGCTGTTACCCAAGTGAGATCGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGAtgtg                                                                                                                                                    <  1:1580762/150‑1 (MQ=255)
  atatTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGAtgtg                                                                                                                                                    <  1:1340072/150‑1 (MQ=255)
    atTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGaa                                                                                                                                                  <  1:458828/150‑1 (MQ=255)
         ggCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAatcatc                                                                                                                                             >  1:1211031/1‑150 (MQ=255)
         ggCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAat                                                                                                                                                 <  1:1208945/146‑1 (MQ=255)
         ggCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAat                                                                                                                                                 >  2:1208945/1‑146 (MQ=255)
          gCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATcc                                                                                                                                            <  1:872605/150‑1 (MQ=255)
             aaTCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATTCACATCATCCGGc                                                                                                                                         <  1:1126363/150‑1 (MQ=255)
             aaTCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATTCACATCATCCGGc                                                                                                                                         <  1:1807244/150‑1 (MQ=255)
              aTCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATTCACATCATCCGGCa                                                                                                                                        <  1:497572/150‑1 (MQ=255)
              aTCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCa                                                                                                                                        >  2:519009/1‑150 (MQ=255)
                 ttGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTg                                                                                                                                     <  2:214054/150‑1 (MQ=255)
                     aTCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGatt                                                                                                                                 >  2:804090/1‑150 (MQ=255)
                         gCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGAttatta                                                                                                                             <  2:1192612/150‑1 (MQ=255)
                          cGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGCAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAACGTGCCGGATGATGTGAATCATCCGGCACGGGATTATTAc                                                                                                                            <  2:1353496/150‑1 (MQ=255)
                           ggAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACt                                                                                                                           >  2:1287717/1‑150 (MQ=255)
                              aaTGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGc                                                                                                                        <  2:2315144/150‑1 (MQ=255)
                              aaTGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGc                                                                                                                        <  2:1066774/150‑1 (MQ=255)
                                 gCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGAt                                                                                                                     >  2:2418861/1‑150 (MQ=255)
                                  cAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGAtt                                                                                                                    >  2:2357545/1‑150 (MQ=255)
                                   aaaaaGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTg                                                                                                                   >  1:928848/1‑150 (MQ=255)
                                   aaaaaGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTg                                                                                                                   >  2:30247/1‑150 (MQ=255)
                                    aaaaGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGt                                                                                                                  <  2:949348/150‑1 (MQ=255)
                                        tcATTGAAATAGGTTGGGGTAATCACCCGTTGATTGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCCACGTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCAtt                                                                                                              <  1:2017891/149‑1 (MQ=255)
                                           ttGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGc                                                                                                           >  2:2035565/1‑150 (MQ=255)
                                            tGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGcc                                                                                                          >  2:1276925/1‑150 (MQ=255)
                                            tGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGcc                                                                                                          <  1:387135/150‑1 (MQ=255)
                                              aaaGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTg                                                                                                        <  1:1489873/150‑1 (MQ=255)
                                              aaaGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTg                                                                                                        >  1:1342127/1‑150 (MQ=255)
                                               aaGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTTATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGa                                                                                                       <  2:133591/150‑1 (MQ=255)
                                                agagTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGCTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGAc                                                                                                      <  2:77554/150‑1 (MQ=255)
                                                  agTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGc                                                                                                    <  1:30247/150‑1 (MQ=255)
                                                  agTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGAATGACAATCGccagtaat                                                                                                    <  1:1702061/150‑7 (MQ=255)
                                                  agTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGAATGACAATCGccagtaat                                                                                                    <  1:1598005/150‑7 (MQ=255)
                                                  agTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGAATGACAATCGccagtaat                                                                                                    <  1:2341459/150‑7 (MQ=255)
                                                   gTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCa                                                                                                   >  2:565458/1‑150 (MQ=255)
                                                      aGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATa                                                                                                >  1:731404/1‑150 (MQ=255)
                                                       ggggAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCTCCTGACGCAATaa                                                                                               <  2:1839876/150‑1 (MQ=255)
                                                        gggAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAac                                                                                              <  1:1538353/150‑1 (MQ=255)
                                                          gAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAacac                                                                                            >  2:648585/1‑150 (MQ=255)
                                                             aTCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACAcgcg                                                                                         >  1:1888854/1‑150 (MQ=255)
                                                              tCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCgg                                                                                        >  2:1620454/1‑150 (MQ=255)
                                                               cACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACAcgcg                                                                                         <  1:1959108/148‑1 (MQ=255)
                                                               cACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACAcgcg                                                                                         >  2:1959108/1‑148 (MQ=255)
                                                               cACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGc                                                                                       <  2:225938/150‑1 (MQ=255)
                                                                aCCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCt                                                                                      >  2:645422/1‑150 (MQ=255)
                                                                 cccGTTCAATGAAGGTATTGGGCTGTTACCCAATTGAGAACGTATTGCCTGTTGATCCAACCTGATGAAAATGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGAAGCAAAAAAACGCGGCtt                                                                                     <  2:1900234/150‑1 (MQ=255)
                                                                     ttCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCAc                                                                                 <  1:1409592/150‑1 (MQ=255)
                                                                          tGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGa                                                                            <  2:1800458/150‑1 (MQ=255)
                                                                                tATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAAt                                                                                                 <  2:1468962/123‑1 (MQ=255)
                                                                                tATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGc                                                                      <  2:1335282/150‑1 (MQ=255)
                                                                                  ttGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCtg                                                                    <  1:746869/150‑1 (MQ=255)
                                                                                  ttGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCtg                                                                    >  2:780795/1‑150 (MQ=255)
                                                                                   tGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCtgt                                                                   >  1:263813/1‑150 (MQ=255)
                                                                                    gggCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCtgtg                                                                  >  1:1216287/1‑150 (MQ=255)
                                                                                     ggCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGc                                                                 >  2:1208363/1‑150 (MQ=255)
                                                                                       cTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGt                                                               <  1:539626/150‑1 (MQ=255)
                                                                                       cTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGt                                                               <  2:2399420/150‑1 (MQ=255)
                                                                                              ccAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGcc                                                        <  2:671212/150‑1 (MQ=255)
                                                                                              ccAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGcc                                                        <  2:1635417/150‑1 (MQ=255)
                                                                                              ccAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGcc                                                        <  1:80528/150‑1 (MQ=255)
                                                                                              ccAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGAGGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGcc                                                        <  1:2321280/150‑1 (MQ=255)
                                                                                               cAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCg                                                       >  1:140557/1‑150 (MQ=255)
                                                                                                  gTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAAc                                                    >  2:2396302/1‑150 (MQ=255)
                                                                                                      gaACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGt                                                >  1:935456/1‑150 (MQ=255)
                                                                                                          gTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTc                                            <  1:1620454/150‑1 (MQ=255)
                                                                                                           tAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTcc                                           >  2:796925/1‑150 (MQ=255)
                                                                                                               gCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCt                                       >  2:1411532/1‑150 (MQ=255)
                                                                                                                  tGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGc                                    <  2:1211031/150‑1 (MQ=255)
                                                                                                                  tGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCAGAAAACGCTGTCCGTAATCGCCGAACCAGTGCTCCACCTTGc                                    >  1:978143/1‑150 (MQ=255)
                                                                                                                  tGTTGATCCAACCTGATGAAAAGGTGCAGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGc                                    >  1:986418/1‑150 (MQ=255)
                                                                                                                   gTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCg                                   <  1:1718914/150‑1 (MQ=255)
                                                                                                                     tGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGa                                 <  2:452142/150‑1 (MQ=255)
                                                                                                                      gATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGaa                                <  2:1342127/150‑1 (MQ=255)
                                                                                                                        tCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAAc                              >  1:1463612/1‑150 (MQ=255)
                                                                                                                        tCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAAc                              <  1:1360127/150‑1 (MQ=255)
                                                                                                                         ccAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACt                             >  2:344713/1‑150 (MQ=255)
                                                                                                                          cAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTg                            <  1:648585/150‑1 (MQ=255)
                                                                                                                            aCCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCTCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTc                          <  1:2213884/150‑1 (MQ=255)
                                                                                                                            aCCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTc                          <  1:344709/150‑1 (MQ=255)
                                                                                                                            aCCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTc                          <  1:2127869/150‑1 (MQ=255)
                                                                                                                               tgatgaAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAAt                       >  1:1056164/1‑150 (MQ=255)
                                                                                                                                  tgaAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATaaa                    >  2:1370708/1‑150 (MQ=255)
                                                                                                                                      aaGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGcc                <  1:1276925/150‑1 (MQ=255)
                                                                                                                                         gTGCCGGATGATTCACATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGc             >  1:2406425/1‑150 (MQ=255)
                                                                                                                                          tGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCt            <  2:2076062/150‑1 (MQ=255)
                                                                                                                                          tGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCt            <  2:1380480/150‑1 (MQ=255)
                                                                                                                                           gCCGGATGATTCACATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCtt           >  2:1500102/1‑150 (MQ=255)
                                                                                                                                            ccGGATGATTCACATCATCCGGCACTGGATTATTACTGGCGATTGTCAtt                                                                                                              >  1:143700/1‑50 (MQ=18)
                                                                                                                                            ccGGATGATTCACATCATCCGGCACTGGATTATTACTGGCGATTGTCAtt                                                                                                              <  2:143700/50‑1 (MQ=18)
                                                                                                                                               gatgatGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCGGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCg       <  2:2212217/150‑1 (MQ=255)
                                                                                                                                                atgatgTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGc      <  2:1577823/150‑1 (MQ=255)
                                                                                                                                                  gatgTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGccc    <  2:464458/150‑1 (MQ=255)
                                                                                                                                                   atgTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCt   >  1:374279/1‑150 (MQ=255)
                                                                                                                                                   atgTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCt   <  1:645422/150‑1 (MQ=255)
                                                                                                                                                    tgtgAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTg  <  1:1234297/150‑1 (MQ=255)
                                                                                                                                                    tgtgAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTg  >  1:451359/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGATATTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTG  >  NC_000913/2738776‑2739073

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