Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 375,273 C→T 16.7% intergenic (+392/‑186) mhpE → / → mhpT 4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913375,2730CT16.7% 385.4 / 53.9 177intergenic (+392/‑186)mhpE/mhpT4‑hyroxy‑2‑oxovalerate/4‑hydroxy‑2‑oxopentanoic acid aldolase, class I/3‑hydroxyphenylpropionic transporter
Reads supporting (aligned to +/- strand):  ref base C (62/84);  new base T (5/25);  total (67/110)
Fisher's exact test for biased strand distribution p-value = 7.56e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 183 total aligned reads displayed.

GAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGC  >  NC_000913/375175‑375372
                                                                                                  |                                                                                                   
gagaATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGgcg                                                                                                    >  1:2972445/1‑100 (MQ=38)
gagaATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGACGTAAACGCCTTATCCGGCCTACggcgggcg                                                                                                    >  2:2011363/1‑100 (MQ=255)
     tttGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGa                                                                                               <  1:1581706/100‑1 (MQ=38)
       ccgactgcctgatgcgacgctcgcgcgtcttatcaggcagtcggcactgtTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAt                                                                                             <  1:1243543‑M1/50‑1 (MQ=255)
        cgactgcctgatgcgacgcttgcgcgtcttatcaggcagtcggcaccatTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAAtt                                                                                            <  1:2431427‑M1/51‑1 (MQ=255)
          aGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTg                                                                                          <  2:2653641/100‑1 (MQ=38)
           ggCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGt                                                                                         >  1:2839607/1‑100 (MQ=38)
           ctgcctgatgcgacgcttgcgcgtcttatcaggcagttggcactgtTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGCGGGTGCGGGAATTTGt                                                                                         <  2:1842451‑M1/54‑1 (MQ=255)
           ctgcctgatgcgacgcttgcgcgtcttatcaggcagttggcactgtTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGt                                                                                         <  1:2375050‑M1/54‑1 (MQ=255)
           ctgcctgatgcgacgcttgcgcgtcttatcaggcagttggcactgtTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGt                                                                                         <  2:1251754‑M1/54‑1 (MQ=255)
           ctgcctgatgcgacgcttgcgcgtcttatcaggcagtcggcaccgtTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGt                                                                                         <  1:989228‑M1/54‑1 (MQ=255)
           ctgcctgatgcgacgcttgcgcgtcttatcaggcagtcggcaccgtTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGt                                                                                         <  1:1505131‑M1/54‑1 (MQ=255)
           ctgcctgatgcgacgcttgcgcgtcttatcaggcagtcggcaccgtTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGt                                                                                         <  2:1712940‑M1/54‑1 (MQ=255)
           ctgcctgatgcgacgctcgcgcgtcttatcaggcagtcggcactgtTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGt                                                                                         <  2:1722070‑M1/54‑1 (MQ=255)
           ctgcctgatgcgacgctcgcgcgtcttatcaggcagtcggcactgtTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGt                                                                                         <  1:2970286‑M1/54‑1 (MQ=255)
            tgcctgatgcgacgcttgcccgtcttatcaggcagtcgtcaccgtTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTa                                                                                        <  2:1004970‑M1/55‑1 (MQ=255)
             gcctgatgcgacgcttgcgcgtcttatcaggcagtcggcaccgtTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAg                                                                                       <  1:2876937‑M1/56‑1 (MQ=255)
             gcctgatgcgacgcttgcgcgtcttatcaggcagtcggcaccatTGTCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAg                                                                                       <  2:1730693‑M1/56‑1 (MQ=255)
             gcctgatgcgacgctcgcgcgtcttatcaggcagtcggcactgtTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAg                                                                                       <  1:1250846‑M1/56‑1 (MQ=255)
             ccTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAg                                                                                       >  1:421718/1‑100 (MQ=33)
             ccTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAg                                                                                       <  1:1289233/100‑1 (MQ=35)
              cTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAgg                                                                                      <  1:2087048/100‑1 (MQ=33)
              cTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAgg                                                                                      <  2:111948/100‑1 (MQ=33)
              cTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAg                                                                                       <  2:2345954/99‑1 (MQ=33)
              cTGATAAGACGCGTTAGCGTCGCATCAGGCATCCGCGCACGACTGCCGGAGGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAgg                                                                                      <  1:1853576/100‑1 (MQ=16)
               tGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAGCGCCTTATCCCGCCTACGGATCGCGCCGGCATTTGt                                                                                         >  2:2637589/1‑96 (MQ=11)
                tcttgcgacgcttgcgcgtcttatcaggcagtcggcaccgtTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGTGCGGGAATTTGTAGGcc                                                                                    >  1:2226169‑M1/42‑100 (MQ=255)
                gATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGcc                                                                                    <  2:680144/100‑1 (MQ=33)
                 aTAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCt                                                                                   >  1:892498/1‑100 (MQ=33)
                  tAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTg                                                                                  <  1:2304290/100‑1 (MQ=33)
                  ntgcgacgcttgcgcgtcttatcaggcagttggcactgtTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAggccag                                                                                  <  2:722338‑M1/61‑3 (MQ=255)
                    aGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTCTCCCGCCTACGGATGGCGCGGGAATTTGTATGCCTg                                                                                  >  2:588926/1‑98 (MQ=255)
                    aGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGAt                                                                                >  1:1971608/1‑100 (MQ=33)
                     gACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATa                                                                               >  1:1322395/1‑100 (MQ=33)
                       cgcgTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAg                                                                             >  1:1086100/1‑100 (MQ=33)
                         gggtatgcggctcatcaggcagtcggcactgtTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGGTGCGGGAATTTGTAGGCCTGATa                                                                               >  2:1473414‑M1/33‑96 (MQ=255)
                          gTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGAcg                                                                          <  2:338182/100‑1 (MQ=32)
                              gCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGtt                                                                      >  1:2267468/1‑100 (MQ=33)
                                 tCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGc                                                                   <  1:785922/100‑1 (MQ=33)
                                  cGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCg                                                                  <  2:2403354/100‑1 (MQ=33)
                                  cGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCg                                                                  <  1:401608/100‑1 (MQ=33)
                                        aGGCATCTGCTCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGtn                                                                <  2:368685/96‑2 (MQ=17)
                                        aGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCAt                                                            >  1:2348548/1‑100 (MQ=33)
                                        aGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCAt                                                            >  1:774812/1‑100 (MQ=33)
                                         gggaTCTGCGCACGTCTGCCGGATGTGGCGTAAACGCGTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATc                                                           <  2:1031616/97‑1 (MQ=11)
                                         ggCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATc                                                           >  1:682046/1‑100 (MQ=33)
                                         ggCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATc                                                           <  1:1602584/100‑1 (MQ=33)
                                           cATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAg                                                         <  1:603790/100‑1 (MQ=32)
                                           cATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAg                                                         <  2:1171289/100‑1 (MQ=33)
                                            aTCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAgg                                                        <  1:2236511/100‑1 (MQ=32)
                                             tCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGc                                                       <  1:2554920/100‑1 (MQ=33)
                                              cTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTCCGGATGGCGCGGGAGTTTGTAGGCCTGTTAAGACGCGTTAGCGGCGCATCAGGCa                                                      >  2:2550908/1‑100 (MQ=11)
                                               tGCGCANGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCAt                                                     >  2:1578047/1‑100 (MQ=17)
                                               tGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCAt                                                     >  2:309549/1‑100 (MQ=33)
                                                   ctcGACTGCCGGATGCGGCGTACACGCCTTATCCGATCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATc                                                    <  2:1248653/95‑1 (MQ=11)
                                                    aCGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGca                                                <  1:3096695/100‑1 (MQ=33)
                                                     cGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGcac                                               >  2:2815122/1‑100 (MQ=35)
                                                      gACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGcaca                                              <  2:2439536/100‑1 (MQ=35)
                                                      gACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGcaca                                              <  1:156379/100‑1 (MQ=35)
                                                      gACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGcaca                                              <  1:1917865/100‑1 (MQ=38)
                                                       aCTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCacac                                             >  2:2613223/1‑100 (MQ=38)
                                                       aCTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCacac                                             >  2:2282611/1‑100 (MQ=35)
                                                        cTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACg                                            <  2:2768978/100‑1 (MQ=35)
                                                        cTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACg                                            >  1:167893/1‑100 (MQ=38)
                                                          gCCGGATGCGGCGTAAACGCCCTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACg                                            <  2:2190636/98‑1 (MQ=25)
                                                           ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACt                                         >  1:2245954/1‑100 (MQ=38)
                                                           ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACt                                         >  1:2748651/1‑100 (MQ=38)
                                                           ccGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACt                                         >  1:2852597/1‑100 (MQ=38)
                                                            ccgatctggcgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTg                                        <  1:946952/93‑1 (MQ=255)
                                                               aTGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCg                                     <  2:1362088/100‑1 (MQ=38)
                                                               aTGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCg                                     <  1:1079908/100‑1 (MQ=38)
                                                                  cggcgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGAt                                  <  1:2192185/100‑1 (MQ=38)
                                                                  cggcgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGAt                                  >  1:2579138/1‑100 (MQ=38)
                                                                  cggcgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGCCGCATCAGGCATCTGCACACGACTGCCGGa                                   >  2:2496602/1‑99 (MQ=25)
                                                                   ggcgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTgagatcgg                                 >  2:946948/1‑93 (MQ=255)
                                                                   ggcgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATg                                 <  1:1406867/100‑1 (MQ=38)
                                                                   gccgtactcGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATg                                 <  2:2087881/92‑1 (MQ=255)
                                                                    gcgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGc                                >  2:2035004/1‑100 (MQ=35)
                                                                    gcgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGc                                >  1:861943/1‑100 (MQ=38)
                                                                     cgTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCg                               >  2:1677725/1‑100 (MQ=38)
                                                                      gTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGa                              <  1:1376749/100‑1 (MQ=38)
                                                                      gTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGa                              <  1:184213/100‑1 (MQ=38)
                                                                        aaaCGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATa                            <  1:2265752/100‑1 (MQ=38)
                                                                         aaCGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCg                               >  2:585352/1‑96 (MQ=32)
                                                                            tccTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACg                        <  2:1151749/99‑1 (MQ=255)
                                                                             ccTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACga                       >  1:1151755/1‑99 (MQ=255)
                                                                             ccTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGt                       >  1:2020936/1‑100 (MQ=255)
                                                                              cTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGt                       >  2:224244/1‑99 (MQ=255)
                                                                                   ccGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTc                 >  2:1446151/1‑100 (MQ=255)
                                                                                    cGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTcc                >  2:615271/1‑100 (MQ=255)
                                                                                     ggCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCg               <  2:2972416/100‑1 (MQ=255)
                                                                                       ccTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGATAAACGCCt                     >  2:1899432/1‑92 (MQ=16)
                                                                                        cTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGcn            <  1:57997/100‑2 (MQ=255)
                                                                                           cGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTAc         >  1:3014161/1‑100 (MQ=255)
                                                                                            ggATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACa        <  1:2395363/100‑1 (MQ=255)
                                                                                             gATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACAt       <  2:39568/100‑1 (MQ=255)
                                                                                             gATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACAt       >  2:768501/1‑100 (MQ=255)
                                                                                                  cgcgGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGTAAACGCCCTGTCCGGccccca        >  2:194827/1‑90 (MQ=11)
                                                                                                  cgcgGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGc  <  2:2852568/100‑1 (MQ=255)
                                                                                                  cgcgGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGc  <  1:1467847/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
GAGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCGCACGACTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGATGGCGCGGGAATTTGTAGGCCTGATAAGACGCGTTAGCGTCGCATCAGGCATCTGCACACGACTGCCGGATGCGATAAACGTCTTGTCCGGCCTACATTTCGC  >  NC_000913/375175‑375372

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