Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 124,859 G→C 5.2% A615P (GCT→CCT)  aceE → pyruvate dehydrogenase E1 component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913124,8590GC 5.2% 263.1 / 4.4 77A615P (GCT→CCT) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base G (33/40);  new base C (2/2);  total (35/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

CTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCT  >  NC_000913/124712‑125008
                                                                                                                                                   |                                                                                                                                                     
cTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCt                                                                                                                                                     <  2:809322/150‑1 (MQ=255)
 tGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTg                                                                                                                                                    <  2:947824/150‑1 (MQ=255)
 tGCAGGAAGGGAGCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGAGCGCGGTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTg                                                                                                                                                    <  1:630427/150‑1 (MQ=255)
          ggATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGTCGCCAGCCGCCCAGc                                                                                                                                           <  1:983795/150‑1 (MQ=255)
          ggATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGc                                                                                                                                           <  2:1559135/150‑1 (MQ=255)
           gATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCa                                                                                                                                          >  1:1175459/1‑150 (MQ=255)
           gATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCa                                                                                                                                          >  1:1183095/1‑150 (MQ=255)
             tCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGTCGCCAGCCGCCCAGCaca                                                                                                                                        <  2:546324/150‑3 (MQ=255)
               aaCGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAgcg                                                                                                                                      >  2:1806270/1‑150 (MQ=255)
               aaCGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAgcg                                                                                                                                      <  1:663232/150‑1 (MQ=255)
               aaCGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAgcg                                                                                                                                      >  1:449078/1‑150 (MQ=255)
                aCGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAgcgc                                                                                                                                     >  1:1087185/1‑150 (MQ=255)
                 cGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGcgcg                                                                                                                                    <  1:1944308/150‑1 (MQ=255)
                 cGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGcgcg                                                                                                                                    >  1:678335/1‑150 (MQ=255)
                    gCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCGAGCGCGTgg                                                                                                                                 >  1:341975/1‑150 (MQ=255)
                       gggCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCtt                                                                                                                              >  1:1476338/1‑150 (MQ=255)
                          cgcAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCt                                                                                                                           <  2:1987193/150‑1 (MQ=255)
                           gcAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGTCTTCCTg                                                                                                                          <  2:1692340/150‑1 (MQ=255)
                              ggTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATc                                                                                                                       >  2:1517614/1‑150 (MQ=255)
                              ggTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATc                                                                                                                       >  2:668515/1‑150 (MQ=255)
                              ggTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATc                                                                                                                       >  1:685369/1‑150 (MQ=255)
                               gttgttCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATcg                                                                                                                      <  2:1893067/150‑1 (MQ=255)
                                    tCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGt                                                                                                                 >  2:794884/1‑150 (MQ=255)
                                    tCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGt                                                                                                                 <  1:1806270/150‑1 (MQ=255)
                                    tCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGt                                                                                                                 >  1:971059/1‑150 (MQ=255)
                                     ccTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTa                                                                                                                <  1:1318286/150‑1 (MQ=255)
                                         gctggcAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTc                                                                                                            >  2:723360/1‑150 (MQ=255)
                                         gctggcAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTc                                                                                                            >  2:46512/1‑150 (MQ=255)
                                                 cggcgACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTa                                                                                                    <  2:1086458/150‑1 (MQ=255)
                                                   gcgACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTacc                                                                                                  <  1:680240/150‑1 (MQ=255)
                                                    cgACCTCTTACAGCACCAACAATCTGCCGATGATCCCTTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTacca                                                                                                 <  1:323185/150‑1 (MQ=255)
                                                    cgACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTacca                                                                                                 <  2:147688/150‑1 (MQ=255)
                                                         tctTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTg                                                                                            <  2:1474765/150‑1 (MQ=255)
                                                          ctTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGa                                                                                           <  2:1639585/150‑1 (MQ=255)
                                                           ttACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGaa                                                                                          <  1:1525238/150‑1 (MQ=255)
                                                            tACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAAc                                                                                         <  1:38566/150‑1 (MQ=255)
                                                              cAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACgg                                                                                       >  1:419728/1‑150 (MQ=255)
                                                               aGCACCAACAATCTGCCGATGATCCCGTTCTACGTCTATTACTCGTTGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGc                                                                                      <  2:942538/150‑1 (MQ=255)
                                                                gCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCg                                                                                     >  2:455056/1‑150 (MQ=255)
                                                                 cACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGa                                                                                    >  2:1990570/1‑150 (MQ=255)
                                                                     aacaaTCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGt                                                                                >  2:1833140/1‑150 (MQ=255)
                                                                         aTCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTgc                                                                            <  1:668515/150‑1 (MQ=255)
                                                                         aTCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTgc                                                                            <  2:1070671/150‑1 (MQ=255)
                                                                            tGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTgcagc                                                                         >  1:479220/1‑150 (MQ=255)
                                                                                gatgatCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGa                                                                     >  1:969586/1‑150 (MQ=255)
                                                                                gatgatCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGa                                                                     <  2:479220/150‑1 (MQ=255)
                                                                                  tgatCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAg                                                                   <  2:1169974/150‑1 (MQ=255)
                                                                                   gatCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCGAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGa                                                                  <  2:341975/150‑1 (MQ=255)
                                                                                        cGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTc                                                             <  2:449078/150‑1 (MQ=255)
                                                                                        cGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTc                                                             <  1:1470269/150‑1 (MQ=255)
                                                                                                aTCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCcac                                                     <  1:455056/150‑1 (MQ=255)
                                                                                                   tATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACAtt                                                  >  1:733798/1‑150 (MQ=255)
                                                                                                   tATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACAtt                                                  >  1:1743260/1‑150 (MQ=255)
                                                                                                    aTTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTc                                                 <  2:782905/150‑1 (MQ=255)
                                                                                                       aCTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGTCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTTAACTGCGAAGGTCTGGAGCAGTATGATGGTCACAGCGACATTCatt                                              <  2:831099/150‑3 (MQ=255)
                                                                                                        cTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTcaca                                                          >  1:423033/1‑137 (MQ=255)
                                                                                                        cTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTcaca                                                          <  2:423033/137‑1 (MQ=255)
                                                                                                        atCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTc                                             >  1:950142/2‑150 (MQ=255)
                                                                                                         tCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCg                                            <  2:240052/150‑1 (MQ=255)
                                                                                                         tCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCg                                            >  1:420890/1‑150 (MQ=255)
                                                                                                           gATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCt                                          <  1:631133/150‑1 (MQ=255)
                                                                                                              gTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGAc                                       >  1:1957900/1‑150 (MQ=255)
                                                                                                              gTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGAc                                       <  1:1517614/150‑1 (MQ=255)
                                                                                                                  ggCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATc                                   <  1:499258/150‑1 (MQ=255)
                                                                                                                     ttCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCg                                <  1:1163381/150‑1 (MQ=255)
                                                                                                                     ttCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCg                                <  2:1037277/150‑1 (MQ=255)
                                                                                                                          gCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTg                           <  2:1906959/150‑1 (MQ=255)
                                                                                                                          gCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTg                           >  1:197314/1‑150 (MQ=255)
                                                                                                                          gCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTg                           >  1:1395286/1‑150 (MQ=255)
                                                                                                                             tATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTAt                        >  1:1557675/1‑150 (MQ=255)
                                                                                                                                 ggCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCTAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATctct                    <  1:46512/150‑1 (MQ=255)
                                                                                                                                    gATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTAc                 <  2:685369/150‑1 (MQ=255)
                                                                                                                                       gcttgCTGGTCGCCAGCCGCCCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAc              >  1:879023/4‑150 (MQ=255)
                                                                                                                                       gcttgCTGGTCGCCAGCCGCCCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAc              >  1:1981556/4‑150 (MQ=255)
                                                                                                                                         gtgCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAccc            >  2:1992084/1‑150 (MQ=255)
                                                                                                                                         gtgCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCATGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAccc            >  1:1575706/1‑150 (MQ=255)
                                                                                                                                              gggCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgctt       >  2:1472606/1‑150 (MQ=255)
                                                                                                                                               ggcggcTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgctta      <  1:1990570/150‑1 (MQ=255)
                                                                                                                                                 cggcTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgcttacg    >  1:1463814/1‑150 (MQ=255)
                                                                                                                                                   gcTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgcttacgct  <  1:723360/150‑1 (MQ=255)
                                                                                                                                                   gcTGGCGACCAGCAAGCGCGTGGCGTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAAGTGTATCTCTTACGACCCGgcttacgct  <  2:673455/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTGCAGGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCT  >  NC_000913/124712‑125008

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: