Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,738,926 T→C 6.9% intergenic (+21/+22) 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,9260TC 6.9% 426.3 / 15.7 144intergenic (+21/+22)pheA/tyrAfused chorismate mutase/prephenate dehydratase/fused chorismate mutase/prephenate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (57/77);  new base C (5/5);  total (62/82)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

Only 100 of 152 total aligned reads displayed.

GATATTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTGCT  >  NC_000913/2738777‑2739075
                                                                                                                                                     |                                                                                                                                                     
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)
 atatTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGAtgtg                                                                                                                                                      <  2:2313988/150‑1 (MQ=255)
   atTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGaa                                                                                                                                                    <  1:458828/150‑1 (MQ=255)
    ttCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAat                                                                                                                                                   <  1:1800213/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)
           cAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATTCACATCATCCgg                                                                                                                                            <  2:778852/150‑1 (MQ=255)
            aaTCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATTCACATCATCCGGc                                                                                                                                           <  1:1126363/150‑1 (MQ=255)
             aTCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCa                                                                                                                                          >  2:519009/1‑150 (MQ=255)
             aTCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCa                                                                                                                                          >  1:421562/1‑150 (MQ=255)
                    aTCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGatt                                                                                                                                   >  2:804090/1‑150 (MQ=255)
                      cAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCTTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGattat                                                                                                                                 <  2:1404735/150‑1 (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:1066774/150‑1 (MQ=255)
                                gCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGAt                                                                                                                       >  1:1635417/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)
                                           tGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGcc                                                                                                            >  2:1276925/1‑150 (MQ=255)
                                           tGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGcc                                                                                                            <  1:387135/150‑1 (MQ=255)
                                             aaaGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTg                                                                                                          >  1:1342127/1‑150 (MQ=255)
                                             aaaGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTg                                                                                                          <  1:1489873/150‑1 (MQ=255)
                                              aaGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTTATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGa                                                                                                         <  2:133591/150‑1 (MQ=255)
                                              aaGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGAATGACAATCGccagt                                                                                                         <  2:1455878/150‑4 (MQ=255)
                                                 agTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGc                                                                                                      <  1:30247/150‑1 (MQ=255)
                                                 agTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGc                                                                                                      >  2:654222/1‑150 (MQ=255)
                                                 agTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGAATGACAATCGccagtaat                                                                                                      <  1:2341459/150‑7 (MQ=255)
                                                  gTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGAATGACAATCGccagtaata                                                                                                     <  2:2251597/150‑8 (MQ=255)
                                                     aGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATa                                                                                                  >  1:731404/1‑150 (MQ=255)
                                                     aGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATa                                                                                                  <  2:1580126/150‑1 (MQ=255)
                                                       gggAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAac                                                                                                <  1:1538353/150‑1 (MQ=255)
                                                        ggAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAaca                                                                                               >  2:2033034/1‑150 (MQ=255)
                                                            aTCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACAcgcg                                                                                           >  2:1535699/1‑150 (MQ=255)
                                                            aTCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACAcgcg                                                                                           >  1:1888854/1‑150 (MQ=255)
                                                             tCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCgg                                                                                          >  2:1620454/1‑150 (MQ=255)
                                                               aCCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCt                                                                                        >  2:645422/1‑150 (MQ=255)
                                                                cccGTTCAATGAAGGTATTGGGCTGTTACCCAATTGAGAACGTATTGCCTGTTGATCCAACCTGATGAAAATGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGAAGCAAAAAAACGCGGCtt                                                                                       <  2:1900234/150‑1 (MQ=255)
                                                                   gTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCa                                                                                    >  2:1718914/1‑150 (MQ=255)
                                                                    ttCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCAc                                                                                   <  1:1409592/150‑1 (MQ=255)
                                                                         tGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGa                                                                              <  2:1800458/150‑1 (MQ=255)
                                                                               tATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGc                                                                        >  2:5055/1‑150 (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)
                                                                                  tGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCtgt                                                                     <  1:352923/150‑1 (MQ=255)
                                                                                   gggCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCtgtg                                                                    >  1:1216287/1‑150 (MQ=255)
                                                                                    ggCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGc                                                                   >  2:1208363/1‑150 (MQ=255)
                                                                                     gCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCg                                                                  <  2:264551/150‑1 (MQ=255)
                                                                                      cTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGt                                                                 <  2:2399420/150‑1 (MQ=255)
                                                                                      cTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGt                                                                 <  1:539626/150‑1 (MQ=255)
                                                                                             ccAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGcc                                                          <  2:671212/150‑1 (MQ=255)
                                                                                             ccAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGcc                                                          >  1:1922746/1‑150 (MQ=255)
                                                                                             ccAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGcc                                                          <  1:80528/150‑1 (MQ=255)
                                                                                             ccAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGAGGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGcc                                                          <  1:2321280/150‑1 (MQ=255)
                                                                                              cAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCg                                                         <  2:2368037/150‑1 (MQ=255)
                                                                                              cAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCg                                                         >  1:140557/1‑150 (MQ=255)
                                                                                              cAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCg                                                         <  2:803104/150‑1 (MQ=255)
                                                                                                 gTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAAc                                                      >  1:459157/1‑150 (MQ=255)
                                                                                                     gaACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGt                                                  >  1:935456/1‑150 (MQ=255)
                                                                                                         gTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTc                                              <  1:1620454/150‑1 (MQ=255)
                                                                                                          tAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTcc                                             <  1:1595967/150‑1 (MQ=255)
                                                                                                                 tGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGc                                      <  2:1211031/150‑1 (MQ=255)
                                                                                                                 tGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCAGAAAACGCTGTCCGTAATCGCCGAACCAGTGCTCCACCTTGc                                      >  1:978143/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                                <  2:421562/150‑1 (MQ=255)
                                                                                                                       tCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAAc                                <  2:576234/150‑1 (MQ=255)
                                                                                                                       tCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAAc                                <  1:1360127/150‑1 (MQ=255)
                                                                                                                       tCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAAc                                >  1:1463612/1‑150 (MQ=255)
                                                                                                                        ccAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCTCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACt                               >  2:2128648/1‑150 (MQ=255)
                                                                                                                         cAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTg                              <  1:648585/150‑1 (MQ=255)
                                                                                                                         cAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACACAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTg                              >  2:2413473/1‑150 (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)
                                                                                                                                     aaGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGcc                  <  1:1276925/150‑1 (MQ=255)
                                                                                                                                        gTGCCGGATGATTCACATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGc               >  2:1043806/1‑150 (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:186689/1‑150 (MQ=255)
                                                                                                                                          gCCGGATGATTCACATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCtt             >  2:1500102/1‑150 (MQ=255)
                                                                                                                                           ccGGATGATTCACATCATCCGGCACTGGATTATTACTGGCGATTGTCAtt                                                                                                                <  2:143700/50‑1 (MQ=18)
                                                                                                                                           ccGGATGATTCACATCATCCGGCACTGGATTATTACTGGCGATTGTCAtt                                                                                                                >  1:143700/1‑50 (MQ=18)
                                                                                                                                              gatgatGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCGGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCg         <  2:2212217/150‑1 (MQ=255)
                                                                                                                                                  atgTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCt     <  1:645422/150‑1 (MQ=255)
                                                                                                                                                   tgtgAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTg    <  1:1234297/150‑1 (MQ=255)
                                                                                                                                                     tgAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTGCt  >  1:521456/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATATTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTGCT  >  NC_000913/2738777‑2739075

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: