Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 124,866 A→C 5.7% D617A (GAC→GCC)  aceE → pyruvate dehydrogenase E1 component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913124,8660AC 5.7% 215.0 / 4.9 70D617A (GAC→GCC) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base A (37/29);  new base C (2/2);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

GGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTA  >  NC_000913/124717‑125010
                                                                                                                                                     |                                                                                                                                                
ggAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGa                                                                                                                                                  >  1:375029/1‑150 (MQ=255)
    gggATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGTCGCCAGCCGCCCAg                                                                                                                                              <  2:1121162/150‑1 (MQ=255)
      gATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCa                                                                                                                                            >  1:52782/1‑150 (MQ=255)
        tCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGTCGCCAGCCGCCCAGCaca                                                                                                                                          <  1:976462/150‑3 (MQ=255)
              aGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCACGTTCTACATCTATTACTAGATGTTCGGCTTCCAGCGTATTGGCTATCTGTGCTGGTCGGCTGGCGACCAGCAAGCGCGTg                                                                                                                                    <  2:1388310/150‑1 (MQ=255)
                     cgcAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCCTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCt                                                                                                                             >  1:518969/1‑150 (MQ=255)
                         ggTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATc                                                                                                                         <  1:1326531/150‑1 (MQ=255)
                             gttCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATcggcg                                                                                                                     <  2:91155/150‑1 (MQ=255)
                             gttCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATcggcg                                                                                                                     >  1:917610/1‑150 (MQ=255)
                               tCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGt                                                                                                                   <  2:375029/150‑1 (MQ=255)
                               tCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGt                                                                                                                   >  1:512934/1‑150 (MQ=255)
                                ccTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTa                                                                                                                  >  1:73272/1‑150 (MQ=255)
                                 ctggctggCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTAc                                                                                                                 <  1:212200/150‑1 (MQ=255)
                                  tggctggcAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACCTCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACt                                                                                                                >  2:1212334/1‑150 (MQ=255)
                                  tggctggcAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACt                                                                                                                <  2:892261/150‑1 (MQ=255)
                                      tggcAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCg                                                                                                            <  1:589647/150‑1 (MQ=255)
                                          aGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCATCGTATTGTCTATCTGTGCTGGGAGTCTTTCGAGCATAAAGAGCGTGGCTTCATTATAGGCTGTATTTCCTGTCg                                                                                                        <  2:683248/150‑1 (MQ=255)
                                          aGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCg                                                                                                        <  1:652477/150‑1 (MQ=255)
                                            cggcgACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTa                                                                                                      <  1:32631/150‑1 (MQ=255)
                                            cggcgACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTa                                                                                                      <  2:556034/150‑1 (MQ=255)
                                            cggcgACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTa                                                                                                      <  2:121987/150‑1 (MQ=255)
                                                gACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTaccac                                                                                                  >  2:1141327/1‑150 (MQ=255)
                                                     ctTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGa                                                                                             >  1:695452/1‑150 (MQ=255)
                                                      ttACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGaa                                                                                            >  1:203960/1‑150 (MQ=255)
                                                           gCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCg                                                                                       >  2:396661/1‑150 (MQ=255)
                                                           gCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCg                                                                                       >  1:1040788/1‑150 (MQ=255)
                                                           gCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCCAACAGCAAGCGCGTGCCTTCCTGATCGGCGGTACTTCCGCTCGTACCACCCTGAACGGca                                                                                       >  2:642620/1‑149 (MQ=255)
                                                              ccAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTGCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAg                                                                                    <  2:169673/150‑1 (MQ=255)
                                                              ccAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAg                                                                                    >  1:124292/1‑150 (MQ=255)
                                                               caacaaTCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAgg                                                                                   >  1:1211596/1‑150 (MQ=255)
                                                                aacaaTCTGCCGATGATACCGTTCTACATCAATTACTCGATGTTCGGATTCCAGCGTATTGGCGATCTGTGCAGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGt                                                                                  >  2:989705/1‑150 (MQ=255)
                                                                    aTCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTgc                                                                              <  2:512934/150‑1 (MQ=255)
                                                                               atCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCTGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGAt                                                                   <  1:185070/150‑1 (MQ=255)
                                                                                tCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATg                                                                  >  1:1270154/1‑150 (MQ=255)
                                                                                 cccGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATgg                                                                 <  2:510074/150‑1 (MQ=255)
                                                                                   cGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTc                                                               >  1:125816/1‑150 (MQ=255)
                                                                                     ccgaaCATCGAGTAATAGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTcac                                                             >  1:808227/5‑150 (MQ=255)
                                                                                       cTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATcggcgg                                                                                                                    >  2:1239345/1‑93 (MQ=255)
                                                                                       cTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATcggcgg                                                                                                                    <  1:1239345/93‑1 (MQ=255)
                                                                                         aCATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGcc                                                         >  1:240305/1‑150 (MQ=255)
                                                                                           aTCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCcac                                                       <  1:1141327/150‑1 (MQ=255)
                                                                                           aTCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGAGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCcac                                                       >  1:1197836/1‑150 (MQ=255)
                                                                                            tCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCcaca                                                      >  2:266991/1‑150 (MQ=255)
                                                                                             cTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACAt                                                     >  2:1038734/1‑150 (MQ=255)
                                                                                               aTTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTc                                                   <  2:242787/150‑1 (MQ=255)
                                                                                                ttACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCa                                                  <  1:299475/150‑1 (MQ=255)
                                                                                                  aCTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGt                                                <  2:518969/150‑1 (MQ=255)
                                                                                                     cGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGc                                             >  2:455328/1‑150 (MQ=255)
                                                                                                     cGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGc                                             >  2:593707/1‑150 (MQ=255)
                                                                                                       aTGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTg                                           <  2:1064752/150‑1 (MQ=255)
                                                                                                          ttCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACt                                        >  2:1001851/1‑150 (MQ=255)
                                                                                                           tCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTa                                       >  2:78139/1‑150 (MQ=255)
                                                                                                            cGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTAt                                      >  1:1007774/1‑150 (MQ=255)
                                                                                                              gCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATcc                                    <  1:937537/150‑1 (MQ=255)
                                                                                                                 tCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGa                                 <  1:396661/150‑1 (MQ=255)
                                                                                                                  ccAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGACGCTCACTAACCCGaa                                >  2:140322/1‑150 (MQ=255)
                                                                                                                    aGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACt                              >  2:1360127/1‑150 (MQ=255)
                                                                                                                         aTTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGCACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTAtc                         >  2:882413/1‑150 (MQ=255)
                                                                                                                         aTTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTAtc                         <  2:1007774/150‑1 (MQ=255)
                                                                                                                          ttGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCg                                  >  2:88333/1‑140 (MQ=255)
                                                                                                                          ttGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCg                                  <  1:88333/140‑1 (MQ=255)
                                                                                                                                  gcttgCTGGTCGCCAGCCGCCCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAc                >  1:421008/4‑150 (MQ=255)
                                                                                                                                   cttgCTGGTCGCCAGCCGCCCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGAcc               >  1:784717/3‑150 (MQ=255)
                                                                                                                                      gCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg            >  2:1296064/1‑150 (MQ=255)
                                                                                                                                      acTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCgg            >  1:1127106/2‑150 (MQ=255)
                                                                                                                                       cTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgc           >  2:864454/1‑150 (MQ=255)
                                                                                                                                        tGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgct          <  2:125816/150‑1 (MQ=255)
                                                                                                                                         gggCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTAc                   >  2:1295189/1‑140 (MQ=255)
                                                                                                                                         gggCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTAc                   <  1:1295189/140‑1 (MQ=255)
                                                                                                                                           gcggcTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCTAAGGTCTGCAGCACGAAGATGGTCACAGCCTCATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgcttac       <  1:989705/150‑1 (MQ=255)
                                                                                                                                                tGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGgcttacgctta  >  2:991496/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GGAAGGGATCAACGAGCTGGGCGCAGGTTGTTCCTGGCTGGCAGCGGCGACCTCTTACAGCACCAACAATCTGCCGATGATCCCGTTCTACATCTATTACTCGATGTTCGGCTTCCAGCGTATTGGCGATCTGTGCTGGGCGGCTGGCGACCAGCAAGCGCGTGGCTTCCTGATCGGCGGTACTTCCGGTCGTACCACCCTGAACGGCGAAGGTCTGCAGCACGAAGATGGTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTA  >  NC_000913/124717‑125010

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