Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,738,927 G→A 5.4% intergenic (+22/+21) 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,9270GA 5.4% 353.1 / 6.0 110intergenic (+22/+21)pheA/tyrAfused chorismate mutase/prephenate dehydratase/fused chorismate mutase/prephenate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (42/62);  new base A (3/3);  total (45/65)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.19e-01

Only 100 of 120 total aligned reads displayed.

ATATTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTGCTC  >  NC_000913/2738778‑2739076
                                                                                                                                                     |                                                                                                                                                     
atatTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGAtgtg                                                                                                                                                       <  2:427697/150‑1 (MQ=255)
atatTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGAtgtg                                                                                                                                                       <  2:1395172/150‑1 (MQ=255)
atatTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGAtgtg                                                                                                                                                       <  2:1342521/150‑1 (MQ=255)
atatTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGAtgtg                                                                                                                                                       <  2:1138284/150‑1 (MQ=255)
 tatTCAGGCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGa                                                                                                                                                      >  1:665969/1‑150 (MQ=255)
      aGGCCAATCTTGAATCAGCGGAAATGCAATAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAatcat                                                                                                                                                 <  2:766566/150‑1 (MQ=255)
        gCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATcc                                                                                                                                               <  2:925442/150‑1 (MQ=255)
        gCCAATCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATcc                                                                                                                                               <  1:418488/150‑1 (MQ=255)
           aaTCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATTCACATCATCCGGc                                                                                                                                            <  2:1069513/150‑1 (MQ=255)
           aaTCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATTCACATCATCCGGc                                                                                                                                            <  2:284302/150‑1 (MQ=255)
            aTCTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGATTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATTCACATCATCCGGCa                                                                                                                                           <  2:1650612/150‑1 (MQ=255)
              cTTGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACt                                                                                                                                         >  1:138450/1‑150 (MQ=255)
               ttGAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGTGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTg                                                                                                                                        <  1:1168970/150‑1 (MQ=255)
                 gAATCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGa                                                                                                                                      >  1:567718/1‑150 (MQ=255)
                  aaTCAGCTGAAATGCAAAAAGCATTGAAAGAGTTATTGTAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACTTTGTGCCTGTTGATCCATCCTGATTAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGat                                                                                                                                     <  2:201000/150‑1 (MQ=255)
                  aaTCAGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGat                                                                                                                                     <  1:935196/150‑1 (MQ=255)
                     cAGCGGAAATCCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGattat                                                                                                                                  <  1:869211/150‑1 (MQ=255)
                      aGCGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGattatt                                                                                                                                 >  1:466607/1‑150 (MQ=255)
                       gccgAAATGCAAAAAGCATTGAAAGAGTTATGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTCCCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTTATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGAttatta                                                                                                                                <  2:1433569/147‑1 (MQ=255)
                        cGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGc                                                                                                                                            <  2:639881/137‑1 (MQ=255)
                        cGGAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTAc                                                                                                                               >  2:1320357/1‑150 (MQ=255)
                          gAAATGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTg                                                                                                                             <  2:1185320/150‑1 (MQ=255)
                            aaTGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATTAAGGTATTTGTCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTTGc                                                                                                                           <  2:240654/150‑1 (MQ=255)
                            aaTGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGc                                                                                                                           >  2:734610/1‑150 (MQ=255)
                             aTGCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCg                                                                                                                          >  1:1453985/1‑150 (MQ=255)
                               gCAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGAt                                                                                                                        <  2:1600162/150‑1 (MQ=255)
                                cAAAAAGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGAtt                                                                                                                       >  1:482793/1‑150 (MQ=255)
                                 aaaaaGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGAttt                                                                                                                      >  2:493110/1‑149 (MQ=255)
                                 aaaaaGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTg                                                                                                                      >  2:580899/1‑150 (MQ=255)
                                   aaaGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTTATCCAACATGATGAAAAGGTGCCGGATGATGTGAATCATCCGGAACTGGATTATTACTGGCGATTGTc                                                                                                                    <  2:259595/150‑1 (MQ=255)
                                   aaaGCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTc                                                                                                                    <  1:717289/150‑1 (MQ=255)
                                      gCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCAtt                                                                                                                 <  2:1520881/150‑1 (MQ=255)
                                      gCATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCAtt                                                                                                                 <  1:638903/150‑1 (MQ=255)
                                       cATTGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTc                                                                                                                <  1:1046312/150‑1 (MQ=255)
                                         ttGAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGc                                                                                                              >  2:346554/1‑150 (MQ=255)
                                           gAAAGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCt                                                                                                            >  2:423047/1‑150 (MQ=255)
                                            aaaGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTg                                                                                                           >  2:1020565/1‑150 (MQ=255)
                                             aaGAGTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGa                                                                                                          >  2:357970/1‑150 (MQ=255)
                                               gagTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCAATTGTCATTCGCCTGACg                                                                                                        <  2:1241429/150‑1 (MQ=255)
                                                agTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGGTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGc                                                                                                       <  2:380588/150‑1 (MQ=255)
                                                agTTAGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGc                                                                                                       <  2:1403424/150‑1 (MQ=255)
                                                    aGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATa                                                                                                   >  2:29312/1‑150 (MQ=255)
                                                    aGGGGAAATCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATa                                                                                                   >  2:1546740/1‑150 (MQ=255)
                                                         aaaTCACCCGTTCAATTAATGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACAcg                                                                                              <  1:201439/150‑1 (MQ=255)
                                                          aaTCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACAcgc                                                                                             <  1:691197/150‑1 (MQ=255)
                                                          aaTCACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGACTGACGCAATAACAcgc                                                                                             >  1:628664/1‑150 (MQ=255)
                                                             cACCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGc                                                                                          <  1:1546740/150‑1 (MQ=255)
                                                              aCCCGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCt                                                                                         >  2:782159/1‑150 (MQ=255)
                                                                 cGTTCAATGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTAATTCGCCTGACGCAATAACACGCGGCTTTc                                                                                      >  2:1656591/1‑150 (MQ=255)
                                                                     cAATGAAGGTATTTGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCActc                                                                                  <  1:684160/150‑1 (MQ=255)
                                                                       aTGAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTg                                                                                >  2:918164/1‑150 (MQ=255)
                                                                        tGAAGGTATTCGGCTGTTACCCAAGTGAGAACGTAGGGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGa                                                                               <  1:526366/150‑1 (MQ=255)
                                                                         gAAGGTATTGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGaa                                                                              <  2:204093/150‑1 (MQ=255)
                                                                                ttGGGCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCtg                                                                       <  1:883542/150‑1 (MQ=255)
                                                                                    gCTGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCg                                                                   <  2:482793/150‑1 (MQ=255)
                                                                                      tGTTACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTa                                                                 >  1:103457/1‑150 (MQ=255)
                                                                                         tACCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATc                                                              <  2:541850/150‑1 (MQ=255)
                                                                                          aCCCAAGTGAGAACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCg                                                             >  2:442651/1‑150 (MQ=255)
                                                                                                    gaACGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGt                                                   <  1:666883/150‑1 (MQ=255)
                                                                                                     aaCGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAg                                                    <  2:783607/148‑1 (MQ=255)
                                                                                                     aaCGTAGTGCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAg                                                    >  1:783607/1‑148 (MQ=255)
                                                                                                             gCCTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCt                                          >  1:1395316/1‑150 (MQ=255)
                                                                                                              ccTGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCtt                                         <  1:341741/150‑1 (MQ=255)
                                                                                                                tGTTGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGc                                       >  1:1516380/1‑150 (MQ=255)
                                                                                                                   tGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGa                                    <  1:55450/150‑1 (MQ=255)
                                                                                                                   tGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGa                                    <  2:407742/150‑1 (MQ=255)
                                                                                                                   tGATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGa                                    <  1:918164/150‑1 (MQ=255)
                                                                                                                    gATCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGaa                                   >  1:1392720/1‑150 (MQ=255)
                                                                                                                     aTCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGaaa                                  >  1:301549/1‑150 (MQ=255)
                                                                                                                      tCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAAc                                 <  1:275302/150‑1 (MQ=255)
                                                                                                                      tCCAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAAc                                 <  1:357970/150‑1 (MQ=255)
                                                                                                                       ccAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACt                                <  2:138450/150‑1 (MQ=255)
                                                                                                                       ccAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACt                                >  2:1435974/1‑150 (MQ=255)
                                                                                                                       ccAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACt                                <  1:966613/150‑1 (MQ=255)
                                                                                                                        cAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTg                               >  2:49727/1‑150 (MQ=255)
                                                                                                                        cAACCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTg                               >  1:558519/1‑150 (MQ=255)
                                                                                                                          aCCTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTc                             <  2:1157276/150‑1 (MQ=255)
                                                                                                                           ccTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCa                            <  1:249689/150‑1 (MQ=255)
                                                                                                                           ccTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCa                            <  2:1453985/150‑1 (MQ=255)
                                                                                                                            cTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAAc                                 <  1:1022029/144‑1 (MQ=255)
                                                                                                                            cTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAAc                                 >  2:1022029/1‑144 (MQ=255)
                                                                                                                            cTGATGAAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCaa                           <  1:782159/150‑1 (MQ=255)
                                                                                                                              gatgaAAAGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATa                         >  2:370310/1‑150 (MQ=255)
                                                                                                                                tgaAAAGGTGCCGGATGATGTGACTCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATaaa                       <  1:742597/150‑1 (MQ=255)
                                                                                                                                 taaaaGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAAc                      <  1:389494/149‑1 (MQ=255)
                                                                                                                                  aaaaGGTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACg                     <  2:629688/150‑1 (MQ=255)
                                                                                                                                   aaaGGGGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGc                    <  1:408885/150‑1 (MQ=255)
                                                                                                                                      ggTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTg                 <  1:123702/150‑1 (MQ=255)
                                                                                                                                      ggTGCCGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTg                 <  1:1435974/150‑1 (MQ=255)
                                                                                                                                         gCCGGATGATTCACATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCtt              >  2:630254/1‑150 (MQ=255)
                                                                                                                                          ccGGATGATTCACATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCGATAAACGCCTGCTTa             >  2:448860/1‑150 (MQ=255)
                                                                                                                                          ccGGATGATTCACATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTa             >  1:127300/1‑150 (MQ=255)
                                                                                                                                          ccGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTa             <  1:1656591/150‑1 (MQ=255)
                                                                                                                                          ccGGATGATGTGAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTa             >  2:1340960/1‑150 (MQ=255)
                                                                                                                                             gatgatGTGAATCATCCAGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCg          >  2:1090095/1‑150 (MQ=255)
                                                                                                                                                  tgtgAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCAc                                            >  1:517496/1‑111 (MQ=255)
                                                                                                                                                  tgtgAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCAc                                            <  2:517496/111‑1 (MQ=255)
                                                                                                                                                   gtgAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTGc    <  1:1020565/150‑1 (MQ=255)
                                                                                                                                                    tgAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTGCt   >  1:1481576/1‑150 (MQ=255)
                                                                                                                                                     gAATCATCCGGCACTGGATTATTACTGGCGATTGTCATTCGCCTGACGCAATAACACGCGGCTTTCACTCTGAAAACGCTGTGCGTAATCGCCGAACCAGTGCTCCACCTTGCGGAAACTGTCAATAAACGCCTGCTTATCGCCCTGCTc  >  1:515221/1‑150 (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 < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: