Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,739,791 T→G 5.5% T93T (ACA→ACC tyrA ← fused chorismate mutase T/prephenate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,739,7910TG 5.5% 329.5 / 8.2 110T93T (ACA→ACCtyrAfused chorismate mutase T/prephenate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (61/43);  new base G (4/2);  total (65/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.19e-01

Only 100 of 111 total aligned reads displayed.

CATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGACGCGAGGCCAACATAGATGCCTCGCGCTCCGGAACATA  >  NC_000913/2739643‑2739937
                                                                                                                                                    |                                                                                                                                                  
cATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAgtg                                                                                                                                                   <  1:915454/150‑1 (MQ=255)
        acacTATCCGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACCGGCGGCCCATCTGACCGCCACCGCCGACGCTAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAAt                                                                                                                                           <  1:842225/150‑1 (MQ=255)
        acaaTATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGTGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCTCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAAt                                                                                                                                           <  1:1400109/150‑1 (MQ=255)
              tCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTg                                                                                                                                     >  1:899724/1‑150 (MQ=255)
              tCAGCCGCTCGATCCCAGTCATGGTGCGCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTg                                                                                                                                     <  1:998044/150‑1 (MQ=255)
                aGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTTATAACACGAGAGGGTCAGCATCTTCTCGAACAGGCTTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTc                                                                                                                                   <  1:1638336/150‑1 (MQ=255)
                 gCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCCCCTGATAACCCGAGAGGGTCAGCATCTGCTCGCACCGGCGTCCCATCTGCCCGCCGCCGCCGACGATTACCACCGGACGCAGTGACGGAAAAAGTGTTTTAAATCCTTTGTCg                                                                                                                                  >  1:996700/1‑150 (MQ=255)
                   cGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGtt                                                                                                                                <  1:969112/150‑1 (MQ=255)
                      tCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTc                                                                                                                             <  2:430371/150‑1 (MQ=255)
                       cGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCa                                                                                                                            >  1:1125167/1‑150 (MQ=255)
                       cGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGGCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCa                                                                                                                            <  2:1383292/150‑1 (MQ=255)
                       cGATCCCAGTCATGTTGCTACAGAATCCGCACATGATAAACCGAGAGGTTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACAAACGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCa                                                                                                                            >  1:1816780/1‑150 (MQ=255)
                        gATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCGTTGTCGTTTTCAc                                                                                                                           >  2:501784/1‑150 (MQ=255)
                         aTCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGCCGATAACCACCGGACGCAGTGCCGGACAAAGTGTTTTAAATCCTTTGTCGGTTTCCCt                                                                                                                          >  2:1229112/1‑150 (MQ=255)
                           cccAGTCATGTTGCTCGAGCATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGGCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTgg                                                                                                                        <  1:41958/150‑1 (MQ=255)
                             cAGTAATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAg                                                                                                                      <  2:402768/150‑1 (MQ=255)
                               gTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTa                                                                                                                    <  1:147011/150‑1 (MQ=255)
                                  aTGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGa                                                                                                                 <  1:1549789/150‑1 (MQ=255)
                                  aTGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGa                                                                                                                 <  1:1757561/150‑1 (MQ=255)
                                  aTGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGa                                                                                                                 <  1:1291767/150‑1 (MQ=255)
                                      tGCTCCAGAATCCGCACCTTATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCa                                                                                                             <  1:1789101/150‑1 (MQ=255)
                                       gCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCAc                                                                                                            >  2:692657/1‑150 (MQ=255)
                                        cTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACg                                                                                                           <  2:573476/150‑1 (MQ=255)
                                        cTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACg                                                                                                           >  1:820508/1‑150 (MQ=255)
                                              attCCGCACCAGATAACCCTTGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGCTAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCAc                                                                                                     <  2:78330/148‑1 (MQ=255)
                                               aTCCGCCCCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCAcc                                                                                                    <  1:2002165/150‑1 (MQ=255)
                                                  cGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGa                                                                                                 >  1:1104503/1‑150 (MQ=255)
                                                   gCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGAc                                                                                                <  1:933354/150‑1 (MQ=255)
                                                   gCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGGGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGAc                                                                                                <  1:636845/150‑1 (MQ=255)
                                                       cTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCTGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCaa                                                                                            >  1:636649/1‑150 (MQ=255)
                                                        tGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCCCCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGCTTCACGCATCACCCGACGCaaa                                                                                           >  2:487681/1‑150 (MQ=255)
                                                          aTAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAAc                                                                                         <  2:1125167/150‑1 (MQ=255)
                                                           tAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACa                                                                                        >  1:614432/1‑150 (MQ=255)
                                                             aCCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATc                                                                                      >  2:1089393/1‑150 (MQ=255)
                                                              cccGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATcc                                                                                     >  1:1274259/1‑150 (MQ=255)
                                                                caagagGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTGTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTc                                                                                   >  1:132211/3‑150 (MQ=255)
                                                                cGATAGTGTCAGCATCTTATCGAACAGGCGTCCCATCTTACCGCCACCGCCGACGATAACCACCGGACTCAGTGACTGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATAACCCGACGCAAAACATCCTc                                                                                   <  2:870079/150‑1 (MQ=255)
                                                                    agGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAAACCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATc                                                                               <  2:1376249/150‑1 (MQ=255)
                                                                      ggTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAg                                                                             >  1:1150831/1‑150 (MQ=255)
                                                                      ggTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAg                                                                             >  1:986049/1‑150 (MQ=255)
                                                                        tCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGAt                                                                           >  1:1531504/1‑150 (MQ=255)
                                                                         cAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCTCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATc                                                                          <  1:1171228/150‑1 (MQ=255)
                                                                         cAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGCCGCAGTGACGGACAAAGTGTGTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGGATCCCCCGACGGATACCCTCCTGACTCAGATc                                                                          >  2:358628/1‑150 (MQ=255)
                                                                         cAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATc                                                                          <  1:1826091/150‑1 (MQ=255)
                                                                           gCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTg                                                                        >  2:1680340/1‑150 (MQ=255)
                                                                             aTCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTATAAATCCTTTGTCGATTTCACTGGAGTAAGAATCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGc                                                                      >  1:1919221/1‑150 (MQ=255)
                                                                             aTCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGGATTTAAAACACTTTGTccg                                                                                                                                 >  2:1827301/1‑89 (MQ=38)
                                                                             aTCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGGATTTAAAACACTTTGTccg                                                                                                                                 <  1:1827301/91‑3 (MQ=38)
                                                                             aTCTTCGCGAACAGGCGACCCATCTGACCGCCCCCGACGACGATAACCACCGGCCGCTGTGACGGACAACGTGTGTTAAATCCGTTGTCGTTTTCACTGGAGTAAGATTCTCGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGc                                                                      <  2:1252691/150‑1 (MQ=255)
                                                                              tcttctCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCg                                                                     <  1:2004240/150‑1 (MQ=255)
                                                                               cttctCGACCAGGCGGCCCAGCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTATGTCGTTTTCACTGGAGTAAGATTCACGCCTCACCCGACGCAAAACATCCTCAATCAGATCTGGCgg                                                                    <  2:205950/150‑1 (MQ=255)
                                                                               cttctCGAACAGGCGTCCCATCTGACAGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAATATTCACGCATCAACCGACGCAAAACATCCTCAATCAGATCTGGCgg                                                                    >  2:871001/1‑150 (MQ=255)
                                                                                 tctcGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTGGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTa                                                                  <  2:479514/150‑1 (MQ=255)
                                                                                 tctcGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGCCGCAAAACCTCCTCAATCAGATCTGGCGGTa                                                                  >  1:1445565/1‑150 (MQ=255)
                                                                                  ctcGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATCCACGCGTCTCCCGACGCAAAACATCCTCAATCAGATCTGGCGGTac                                                                 >  1:1724685/1‑150 (MQ=255)
                                                                                   tcGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCCCCCGACGCAAAACATCCTCAATCAGATCTGGCGGTacc                                                                >  1:1124536/1‑149 (MQ=255)
                                                                                    cGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCAGACGCAAAACATCCTCAATCAGATCTGGCGGTacac                                                               >  2:909539/1‑150 (MQ=255)
                                                                                     gAACAGGCGTCCCATCTGACCGCCAACGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACCTCCTCAATCAGATCTGGCGGTACAcc                                                              >  2:887649/1‑150 (MQ=255)
                                                                                            cGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCCCTGGAGTAAGATTCACGCCTCACCCGCCGCAAAAGATCCTCAATCAGATCTGGCGGTACACCCAGAGCt                                                       >  1:369816/1‑150 (MQ=255)
                                                                                               cccATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTcc                                                    >  2:154332/1‑150 (MQ=255)
                                                                                                ccATCTGACCGCCACCTCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTATAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCg                                                   <  2:967973/150‑1 (MQ=255)
                                                                                                ccATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTACATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTCCACCCAGAGCTTCCg                                                   >  1:498402/1‑150 (MQ=255)
                                                                                                 cATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGc                                                  >  2:992603/1‑150 (MQ=255)
                                                                                                  aTCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGcc                                                 <  2:537952/150‑1 (MQ=255)
                                                                                                  aTCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATAAGATCTGGCGGTACACCCAGAGCTTCCGcc                                                 <  1:636634/150‑1 (MQ=255)
                                                                                                   tCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGccg                                                >  2:188429/1‑149 (MQ=255)
                                                                                                     tGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCctct                                              >  2:1397398/1‑150 (MQ=255)
                                                                                                     tGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCctct                                              >  1:2005974/1‑150 (MQ=255)
                                                                                                     tGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCcact                                              >  2:1881789/1‑150 (MQ=255)
                                                                                                      gACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGCGCTTCCGCCTCTg                                             >  2:1936010/1‑150 (MQ=255)
                                                                                                      gACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTg                                             >  1:437402/1‑150 (MQ=255)
                                                                                                      gACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTg                                             >  2:1574873/1‑150 (MQ=255)
                                                                                                      gACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTg                                             >  1:628048/1‑150 (MQ=255)
                                                                                                       accgccaccgccCACGATAACCACCGGCCGGACGGACGGACAAACTTTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCTTCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGc                                            <  1:150645/150‑1 (MQ=255)
                                                                                                        ccgccaccgccGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCa                                           <  2:1417161/150‑1 (MQ=255)
                                                                                                        ccgccaccgccGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCa                                           >  2:1278905/1‑150 (MQ=255)
                                                                                                          gccaccgccGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCacg                                         >  2:1334908/1‑150 (MQ=255)
                                                                                                          gccaccgccGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCacg                                         >  1:1884421/1‑150 (MQ=255)
                                                                                                           ccaccgccGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCacga                                        >  1:1252165/1‑150 (MQ=255)
                                                                                                              ccgccgACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGAcgc                                     >  2:144885/1‑150 (MQ=255)
                                                                                                                        tAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGACGCGAGGCCAACa                           >  2:2001239/1‑150 (MQ=255)
                                                                                                                        tAACCACCGGACGCAGTGAAGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGACGCGAGGCCAACa                           >  1:28991/1‑150 (MQ=255)
                                                                                                                           ccaccGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGACGCGAGGCCAACATAg                        >  2:1145875/1‑150 (MQ=255)
                                                                                                                             accGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGACGCGAGGCCAACATAGAt                      <  2:1376943/150‑1 (MQ=255)
                                                                                                                             accGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGACGCGAGGCCAACATAGAt                      >  2:212809/1‑150 (MQ=255)
                                                                                                                              ccGGACGCAGTGACGGACATAGTGTTTTAAATCCTTTGTCGTTTTCACTTGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGACGCGAGGCCATCATAGATg                     <  1:1636603/150‑1 (MQ=255)
                                                                                                                               cGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGACGCGAGGCCAACATAGATGc                    >  1:911068/1‑150 (MQ=255)
                                                                                                                                   cGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGACGCGAGGCCAACATAGATGCCTcg                >  2:1118180/1‑150 (MQ=255)
                                                                                                                                    tctgTGACTGACAAATTGTTTTTAATCCTTTGTCTTTTTCACTGGAGTAAGATTCACGCATCAGCCGACGCAAAACATCCTCATTCAGATCTGGCGGTACACCCAGATCTTCCGCCTCTGCACGACGCGAGGCCAACATAGATGCCTcgc               <  2:1936320/147‑1 (MQ=255)
                                                                                                                                     cAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGACGCGAGGCCAACATAGATGCCTcgcg              >  2:580934/1‑150 (MQ=255)
                                                                                                                                      tgaaaacgACAAAGGATTTAAAACACTTTGTCGTTTTCACTGGAGTAAGATACACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCCGAGCTTCCGCGTCTGCACGCCGCGACGCCACCATTGCTGCCTCgcgc             >  1:1685823/8‑150 (MQ=255)
                                                                                                                                       gTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGACGCGAGGCCAACATAGATGCCTCgcgcg            >  1:68212/1‑149 (MQ=255)
                                                                                                                                       gTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGACGCGAGGCCAACATAGATGCCTCGCGCt            >  1:775072/1‑150 (MQ=255)
                                                                                                                                        aaaacgACAAAGGATTTAAAACACTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGAACACCCAGAGCTTCAGCCTCTGCACGACGCGAGGCCAACATAGATGCCTCGAGCTc           >  1:1460440/6‑150 (MQ=255)
                                                                                                                                         gACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGACGCGAGGCCAACATAGATGCCTCGCGCTcc          >  1:495535/1‑150 (MQ=255)
                                                                                                                                         gACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGACGCGAGGCCAACATAGATGCCTCGCGCTcc          >  2:1093447/1‑150 (MQ=255)
                                                                                                                                         aaacgACAAAGGATTTAAAACACTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGCGCTTCCGCCTCTGCACGACGCGAGGCCAACATAGATGCCTCGCGCTcc          >  1:878585/5‑150 (MQ=255)
                                                                                                                                          aCGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGACGCGAGGCCAACATAGATGCCTCGCGCTCCg         <  2:1953225/150‑1 (MQ=255)
                                                                                                                                            ggACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGACGCGAGGCCAACATAGATGCCTCGCGCTCCGGa       <  2:877402/150‑1 (MQ=255)
                                                                                                                                                 aaGTGTTTTAATTCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGACGCGAGGCCAACATAGATGCCTCGCGCTCCGGAACATa  <  1:1296078/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CATCGGCAACAATATCAGCCGCTCGATCCCAGTCATGTTGCTCCAGAATCCGCACCTGATAACCCGAGAGGGTCAGCATCTTCTCGAACAGGCGTCCCATCTGACCGCCACCGCCGACGATAACCACCGGACGCAGTGACGGACAAAGTGTTTTAAATCCTTTGTCGTTTTCACTGGAGTAAGATTCACGCATCACCCGACGCAAAACATCCTCAATCAGATCTGGCGGTACACCCAGAGCTTCCGCCTCTGCACGACGCGAGGCCAACATAGATGCCTCGCGCTCCGGAACATA  >  NC_000913/2739643‑2739937

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