Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 257,250 A→C T242P (ACC→CCC)  frsA → fermentation‑respiration switch protein; PTS Enzyme IIA(Glc)‑binding protein; pNP‑butyrate esterase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913257,2500AC98.6% 258.6 / ‑3.9 71T242P (ACC→CCC) frsAfermentation‑respiration switch protein; PTS Enzyme IIA(Glc)‑binding protein; pNP‑butyrate esterase activity
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (38/31);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.85e-01

CGTTCCCGACGGTATTAATGTGTGGTGGTCTGGATGCGATGCAGACGGACTATTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCACCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCTAACGTTGCCGTGCGTCTGGCATACCTTGAATCGCCGCGTCTG  >  NC_000913/257104‑257396
                                                                                                                                                  |                                                                                                                                                  
cGTTCCCGACGGTATTAATGTGTGGTGGTCTGGATGCGATGCAGACGGACTATTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCa                                                                                                                                                <  2:991965/151‑1 (MQ=255)
 gTTCCCGACGGTATTAATGTGTGGTGGTCTGGATGCGATGCAGACGGACTATTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCAccc                                                                                                                                                 >  1:1027233/1‑149 (MQ=255)
  ttCCCGACGGTATTAATGTGTGGTGGTCTGGATGCGATGCAGACGGACTATTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAgg                                                                                                                                              >  1:493721/1‑151 (MQ=255)
    cccGACGGTATTAATGTGTGGTGGTCTGGATGCGATGCAGACGGACTATTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGAc                                                                                                                                            <  1:1300299/151‑1 (MQ=255)
     ccGACGGTATTAATGTGTGGTGGTCTGGATGCGATGCAGACGGACTATTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACt                                                                                                                                           <  1:999726/151‑1 (MQ=255)
      cGACGGTATTAATGTGTGGTGGGCTGGATGCGATGCAGACGGACTATTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTc                                                                                                                                          >  1:1043593/1‑151 (MQ=255)
           gTATTAATGTGTGGTGGTCTGGATGCGATGCAGACGGACTATTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGcc                                                                                                                                     <  1:428228/151‑1 (MQ=255)
             aTTAATGTGTGGTGGTCTGGATGCGATGCAGACGGACTATTACAGGCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGGTTTTCTTCAAAATGGAAGCTCCCCCATGACTCCAGCCt                                                                                                                                    <  1:551298/150‑1 (MQ=255)
               tAATGTGTGGTGGTCTGGATGCGATGCAGACGGACTATTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGtt                                                                                                                                 >  2:826676/1‑151 (MQ=255)
                       ggtggtCTGGATGCGATGCAGACGGACTATTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAg                                                                                                                          <  2:1037005/150‑1 (MQ=255)
                        gtggtCTGGATGCGATGCAGACGGACTATTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCa                                                                                                                        <  1:1103172/151‑1 (MQ=255)
                              tGGATGCGATGCAGACGGACTATTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCCGGACTCCAGCCTGTTGCATCAGCACGTCtt                                                                                                                  >  2:735465/1‑151 (MQ=255)
                               ggATGCGATGCAGACGGACTATTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCtt                                                                                                                  <  1:378111/150‑1 (MQ=255)
                               ggATGCGATGCAGACGGACTATTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCt                                                                                                                   <  1:844673/149‑1 (MQ=255)
                                 atgcgatgcAGACGGACTATTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTaa                                                                                                                >  2:69961/1‑150 (MQ=255)
                                   gcgatgcAGACGGACTATTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAg                                                                                                              <  1:258351/150‑1 (MQ=255)
                                    cgatgcAGACGGACTATTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCATCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGgc                                                                                                            >  2:1026767/1‑151 (MQ=255)
                                         cAGACGGACTATTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGc                                                                                                       >  1:981979/1‑151 (MQ=255)
                                           gACGGACTATTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCt                                                                                                     <  1:275551/151‑1 (MQ=255)
                                               gACTATTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACg                                                                                                 >  1:320506/1‑151 (MQ=255)
                                                   aTTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTa                                                                                               <  1:1219139/149‑1 (MQ=255)
                                                   aTTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTa                                                                                               <  1:179741/149‑1 (MQ=255)
                                                   aTTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTAcc                                                                                             >  1:834272/1‑151 (MQ=255)
                                                   aTTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTAcc                                                                                             >  1:1155871/1‑151 (MQ=255)
                                                   aTTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTAcc                                                                                             >  1:436491/1‑151 (MQ=255)
                                                   aTTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTAcc                                                                                             >  2:664128/1‑151 (MQ=255)
                                                        aGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCgtggg                                                                                        <  2:306993/151‑1 (MQ=255)
                                                          ccTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCgtgggtg                                                                                      <  1:189051/151‑1 (MQ=255)
                                                            tGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCgtgggtgg                                                                                     >  2:456411/1‑150 (MQ=255)
                                                                  aaCGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCacac                                                                              <  1:501122/151‑1 (MQ=255)
                                                                     gTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACtcg                                                                           >  2:107160/1‑151 (MQ=255)
                                                                       tATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACtcgc                                                                          <  2:894745/150‑1 (MQ=255)
                                                                        aTTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACtcgcgt                                                                        >  1:257727/1‑151 (MQ=255)
                                                                          tttGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACtcgcgtc                                                                       <  1:1031679/150‑1 (MQ=255)
                                                                           ttGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACtcgcgtcgc                                                                     >  2:738522/1‑151 (MQ=255)
                                                                           ttGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACtcgcgtcg                                                                      >  2:119649/1‑150 (MQ=255)
                                                                            tGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACtcgcgtcgcg                                                                    >  2:92796/1‑151 (MQ=255)
                                                                              cgccgcGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGc                                                                  <  1:617935/151‑1 (MQ=255)
                                                                                ccgcGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCtt                                                               <  2:1245358/151‑1 (MQ=255)
                                                                                atgcccgGCATTGCGATGCTGACTATTGATATGCTGTCGGTGGGATTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCTTGTTGCATCAGCACGTCTTAAAGGCGTTGCTTAACGTACCGTGGGTGGATCCCACTCGCTTCGCCGCCt                                                                >  1:65061/6‑151 (MQ=255)
                                                                                 cgcgcgGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCtt                                                               >  1:892476/1‑151 (MQ=255)
                                                                                 cgcgcgGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCt                                                                <  1:917488/150‑1 (MQ=255)
                                                                                     cgGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGt                                                           >  1:128395/1‑151 (MQ=255)
                                                                                            gCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGtt                                                    >  1:765477/1‑151 (MQ=255)
                                                                                            gCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGt                                                     <  1:362631/150‑1 (MQ=255)
                                                                                             cGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGttt                                                   >  2:859031/1‑151 (MQ=255)
                                                                                               aTGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCg                                                 >  1:185053/1‑151 (MQ=255)
                                                                                                 gCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCgg                                                >  1:211278/1‑150 (MQ=255)
                                                                                                  cTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGgcg                                              >  1:402588/1‑151 (MQ=255)
                                                                                                     aCTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCt                                            >  1:443141/1‑150 (MQ=255)
                                                                                                     aCTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCt                                            >  2:677645/1‑150 (MQ=255)
                                                                                                          tGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCTAACGt                                       <  1:468616/150‑1 (MQ=255)
                                                                                                            atatGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCTAACGTTGc                                    <  2:955034/151‑1 (MQ=255)
                                                                                                            atatGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCTAACGTTGc                                    >  2:554937/1‑151 (MQ=255)
                                                                                                             tatGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCCTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCTAACGTTGc                                    >  2:293697/1‑150 (MQ=255)
                                                                                                              atGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCTAACGTTGcc                                   <  2:104577/150‑1 (MQ=255)
                                                                                                                gCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCTAACGTTGCCGt                                 <  2:86510/150‑1 (MQ=255)
                                                                                                                    tCGGTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCTAACGTTGCCGTGCGTc                            >  1:525530/1‑151 (MQ=255)
                                                                                                                       gTGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCTAACGTTGCCGTGCGTCTgg                         >  1:633888/1‑151 (MQ=255)
                                                                                                                        tGGGCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCTAACGTTGCCGTGCGTCTgg                         >  1:1083011/1‑150 (MQ=255)
                                                                                                                         gggCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCTAACGTTGCCGTGCGTCTGGc                        >  1:104203/1‑150 (MQ=255)
                                                                                                                          ggCTTTTCTTCAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATAACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGCCGCTAACGTTGCCGTGCGTCTGGCAt                      <  2:1163040/151‑1 (MQ=255)
                                                                                                                               ttcttcAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCTAACGTTGCCGTGCGTCTGGCATACCtt                 >  2:817205/1‑151 (MQ=255)
                                                                                                                               ttcttcAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCTAACGTTGCCGTGCGTCTGGCATACCt                  <  2:287264/150‑1 (MQ=255)
                                                                                                                                tcttcAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCTAACGTTGCCGTGCGTCTGGCATACCtt                 <  2:1155944/150‑1 (MQ=255)
                                                                                                                                  ttcAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCTAACGTTGCCGTGCGTCTGGCATACCTTGa               >  2:502075/1‑150 (MQ=255)
                                                                                                                                   tcAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCTAACGTTGCCGTGCGTCTGGCATACCTTGaa              <  2:942571/150‑1 (MQ=255)
                                                                                                                                   tcAAAATGGAAGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCTAACGTTGCCGTGCGTCTGGCATACCTTGa               <  2:493734/149‑1 (MQ=255)
                                                                                                                                            aaGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCTAACGTTGCCGTGCGTCTGGCATACCTTGAATCGCCGCGTc    >  2:524822/1‑151 (MQ=255)
                                                                                                                                             aGCTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCTAACGTTGCCGTGCGTCTGGCATACCTTGAATCGCCGCGTCt   <  2:587031/151‑1 (MQ=255)
                                                                                                                                               cTCCCCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCTAACGTTGCCGTGCGTCTGGCATACCTTGAATCGCCGCGTCTg  >  1:378327/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGTTCCCGACGGTATTAATGTGTGGTGGTCTGGATGCGATGCAGACGGACTATTACAGCCTGTATGAACGTTATTTTGCGCCGCGCGGCATTGCGATGCTGACTATTGATATGCCGTCGGTGGGCTTTTCTTCAAAATGGAAGCTCACCCAGGACTCCAGCCTGTTGCATCAGCACGTCTTAAAGGCGCTGCCTAACGTACCGTGGGTGGATCACACTCGCGTCGCGGCCTTTGGTTTCCGTTTCGGCGCTAACGTTGCCGTGCGTCTGGCATACCTTGAATCGCCGCGTCTG  >  NC_000913/257104‑257396

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: