Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,072,811 A→C 6.0% V254G (GTG→GGG)  yihU ← 3‑sulpholactaldehyde (SLA) reductase, NADH‑dependent; gamma‑hydroxybutyrate dehydrogenase, NADH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,072,8110AC 6.0% 284.1 / 6.9 111V254G (GTG→GGG) yihU3‑sulpholactaldehyde (SLA) reductase, NADH‑dependent; gamma‑hydroxybutyrate dehydrogenase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base A (50/54);  new base C (3/4);  total (53/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

Only 100 of 114 total aligned reads displayed.

CTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCAC  >  NC_000913/4072662‑4072959
                                                                                                                                                     |                                                                                                                                                    
gtttATCCAGTCGTTTAATTTTTTCTTTGGCCGTCATCCCCGCACTGACACGGACCGGTTCCAGAATGGCGGACCAGTCCTGGGGGCCGCGAACCCAGGCGCGCGCCTGGGATTAAACCGACCGTGCGGCGGTACCCAGCGGGAATGGca                                                                                                                                                      >  2:183845/2‑150 (MQ=255)
    aTCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATg                                                                                                                                                  >  1:1463598/1‑150 (MQ=255)
    aTCCAGTCGTTACATTTTTACTTTGGCAGGCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATg                                                                                                                                                  <  1:1281327/150‑1 (MQ=255)
      ccAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCa                                                                                                                                                >  1:251386/1‑150 (MQ=255)
        cgTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGc                                                                                                                                              <  1:899318/149‑1 (MQ=255)
         gTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCt                                                                                                                                             <  1:21266/150‑1 (MQ=255)
         gTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCt                                                                                                                                             <  1:1670106/150‑1 (MQ=255)
            gTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTCCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGt                                                                                                                                          <  2:669550/150‑1 (MQ=255)
            gTTACATTTTTACTTTGGCAGGCATCCCTGAACTGACACGGACCTGTTACAGAATGGATGACCAGTCCTGGCGACAGCGACCAGATGCGCGCGCCTGTCTATAAACCTCCCGTGAGGCGGCCCCCAGCTGCATGTGCCAATGCCGCTggg                                                                                                                                          <  2:1142004/150‑2 (MQ=255)
             ttACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTTGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCACCCGCGGCATTGGCACATGCAGCTGGtt                                                                                                                                         >  2:68795/1‑150 (MQ=255)
             ttACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGtt                                                                                                                                         >  1:794415/1‑150 (MQ=255)
             ttACATTTTTACTTTGGCAGTCATACCGGGACTGACACGGACCTGTTCGAGAATGGCGGACCAGTCCTGGCGACCGCGTCCCGCTGCGCGCGCCTGGCTATAATCCTCCCGTGAGGCGGCCCCCAGCTGCATGTGCCAATGCCGCTgggg                                                                                                                                         <  2:778857/150‑3 (MQ=255)
              tACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCTGCATGTGCCAATGCCGCTggggg                                                                                                                                        <  1:1034405/150‑4 (MQ=255)
              tACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCTGCATGTGCCAATGCCGCTggggg                                                                                                                                        <  1:1593687/150‑4 (MQ=255)
               aCATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTgg                                                                                                                                       >  1:1215175/1‑150 (MQ=255)
                cATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCCGCGGCATTGGCACATGCAGCTGGTTGGc                                                                                                                                      >  1:281756/1‑150 (MQ=255)
                cATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGc                                                                                                                                      >  1:401072/1‑150 (MQ=255)
                cATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACATCCCGTGAGGCGGCCCCCAGCGGCATTGGCACAAGCAGCTGGTTGGc                                                                                                                                      >  2:286129/1‑150 (MQ=255)
                 aTTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAACCCTCCCGTGAGGCAGCCCCCCGCGGCAATGGCACATGCAGCTGGTTGGCg                                                                                                                                     >  2:887774/1‑150 (MQ=255)
                 aTTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCTGCCCCCAGCGGCCTTGGCACATGCAGCTGGTTGGCg                                                                                                                                     >  2:1329663/1‑150 (MQ=255)
                 aTTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATCGCGGACCCGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCg                                                                                                                                     <  2:208333/150‑1 (MQ=255)
                    tttACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACa                                                                                                                                  <  1:788134/150‑1 (MQ=255)
                      tACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATc                                                                                                                                >  2:1463227/1‑150 (MQ=255)
                           tGGCAGTCATCCCGGCACTTACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAggg                                                                                                                           <  2:385325/150‑1 (MQ=255)
                             gCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCg                                                                                                                         >  1:289967/1‑150 (MQ=255)
                             gCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCg                                                                                                                         <  1:1158818/150‑1 (MQ=255)
                              cAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAACCCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGa                                                                                                                        >  2:1581058/1‑150 (MQ=255)
                              cAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGa                                                                                                                        >  1:1481775/1‑150 (MQ=255)
                                ttCATCCCGGCACTGACACGGCCCTGTTCCAGAATGGCGGACCAGTCCGGGCGACCGCGACCCGCCGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATg                                                                                                                      <  1:869606/149‑1 (MQ=255)
                                   aTCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAACCCTCCCGTGAGGCGGCCCCCCGCGGCATTGGCACATGCAGCTGGTTGGCGCGATCAATGGCGATGCCa                                                                                                                   >  2:37651/1‑150 (MQ=255)
                                      ccGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACAGCAAGGGCGATGCCAAGa                                                                                                                >  1:43388/1‑150 (MQ=255)
                                         gcccTGACACGGACCTGATCCAGCATGGCGGACCAGGCCGGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATcc                                                                                                             <  1:327097/147‑1 (MQ=255)
                                          aattGACAAGGACCTGTGCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCt                                                                                                            <  2:120205/147‑1 (MQ=255)
                                           ccTGCCACGGACCTGTTCGAGAATGGCGGACGAGTCCTGGCGAACGCGACGCGCTTCGCGCGCCTGGCTATAAACCTCCCGTGAGGCTGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGATACATCAAGGGGGATGCCAAGATCCtt                                                                                                           <  2:196874/149‑1 (MQ=255)
                                             tGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGAGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTAt                                                                                                         <  1:1721397/150‑1 (MQ=255)
                                              tacacGGACCTGTTCCAGTATTGCGGAGCAGTCCTGGCGACCGCGACGCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCTGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAATGGCGATGCCAAGTTCCTTATg                                                                                                        <  2:1126653/149‑1 (MQ=255)
                                              gACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATg                                                                                                        >  2:943610/1‑150 (MQ=255)
                                                cacGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGGGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATggg                                                                                                      <  1:1431209/150‑1 (MQ=255)
                                                cacGGACCTGTGCCAGAATGGCCGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATggg                                                                                                      <  2:1184837/150‑1 (MQ=255)
                                                 acGGACCTTTTCCAGAATGTCTGACCAGTACTGGCGACCGCTACCCGATTCTCTCGCCTGGCTATAAACCTCCAGTGAGGATGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGAGACATCAAGTTCGATGAAAAGATCCATATGGGc                                                                                                     >  2:158763/1‑150 (MQ=255)
                                                  cGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCCCATGCAGCTGGTTGGCGCCATCCAGGGCGATGCCAAGATCCTTATGGGCa                                                                                                    >  1:83926/1‑150 (MQ=255)
                                                    gACCTGTTCCAGAACGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCGGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAg                                                                                                  <  1:1629206/150‑1 (MQ=255)
                                                      ccTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAgat                                                                                                <  2:794415/150‑1 (MQ=255)
                                                           tCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAgagcgatc                                                                                           <  2:1252546/150‑1 (MQ=255)
                                                            ccAGAATGGCGGACCAGTCCTGGCGACCGCGCCCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCa                                                                                          >  2:391243/1‑150 (MQ=255)
                                                            ccAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCa                                                                                          <  1:164384/150‑1 (MQ=255)
                                                             cAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACCTGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCAt                                                                                         >  1:195724/1‑150 (MQ=255)
                                                              aGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGCCATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATg                                                                                        >  2:1066399/1‑150 (MQ=255)
                                                              aGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATg                                                                                        <  2:147598/150‑1 (MQ=255)
                                                                 aTGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAg                                                                                     >  2:702347/1‑150 (MQ=255)
                                                                      ggACCAGTCCTGGCGAGCGGGACGCGCTGCGCGCGCCTGGCTATAATCCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCggg                                                                                <  2:135172/150‑1 (MQ=255)
                                                                      ggACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCggg                                                                                <  2:59063/150‑1 (MQ=255)
                                                                      ggACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCggg                                                                                <  2:64275/150‑1 (MQ=255)
                                                                      ggACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCggg                                                                                <  2:753290/150‑1 (MQ=255)
                                                                        aCCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACTTCAAGGGCGATGACAAGATCCTTATGGGCAAGATCGATGATGATGGCGGgag                                                                              >  1:592132/1‑150 (MQ=255)
                                                                        aCCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGgag                                                                              >  1:1029877/1‑150 (MQ=255)
                                                                        aCCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGAGGCCAAGCTCCTGATGGGCAAGATCGATCATGAAGGCGGgag                                                                              >  1:592129/1‑150 (MQ=255)
                                                                          cAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCTTGAGGCGGCCCCCAGCGGAATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGaa                                                                            <  2:1255187/150‑1 (MQ=255)
                                                                            gTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAg                                                                          <  2:353899/150‑1 (MQ=255)
                                                                                tGGCGACCGCGTCCCGCTGCGCGCGCCTGGCTATAAACCTCCCGGGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGAtct                                                                      <  2:272256/150‑1 (MQ=255)
                                                                                tGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGCTCGATCATGAAGGCGGGAGAAAGAtct                                                                      >  1:570231/1‑150 (MQ=255)
                                                                                     accgcgaccCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATGATGAAGGCGGGAGAAAGATCTCCGcg                                                                 >  1:177188/1‑149 (MQ=255)
                                                                                     accgcgaccCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCt                                                                 >  1:182943/1‑150 (MQ=255)
                                                                                        gcgaccCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAg                                                              >  1:1704599/1‑150 (MQ=255)
                                                                                        gcgaccCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAg                                                              >  1:1280642/1‑150 (MQ=255)
                                                                                        gcgaccCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAg                                                              >  1:679181/1‑150 (MQ=255)
                                                                                           accCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCTGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGTCAAGAGCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGAc                                                           <  1:1249144/150‑1 (MQ=255)
                                                                                             ccGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGTCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACtt                                                         <  2:46115/150‑1 (MQ=255)
                                                                                              cGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACttt                                                        <  2:381101/150‑1 (MQ=255)
                                                                                              cGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACttt                                                        <  1:140358/150‑1 (MQ=255)
                                                                                                 tGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGTGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACtttgtt                                                     <  1:286129/150‑1 (MQ=255)
                                                                                                 tGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACtttgtt                                                     <  2:1177542/150‑1 (MQ=255)
                                                                                                    gcgcgcCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTgg                                                  <  1:572886/150‑1 (MQ=255)
                                                                                                         ccTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAgg                                             <  2:1505532/150‑1 (MQ=255)
                                                                                                         ccTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAgg                                             >  1:564783/1‑150 (MQ=255)
                                                                                                         ccTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGAAAGATCGATAATGAAGGGGGTAGAAAGATCTCCGATGAGGAATTTGTTTGGCAAgg                                             >  2:933205/1‑150 (MQ=255)
                                                                                                             gCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAgt                                         <  2:945316/150‑1 (MQ=255)
                                                                                                                 taAACCTCCCGTGAGGCTGCCCCCAGCGGCATTGGCACATGCATCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAgtggtg                                     <  2:1013895/150‑1 (MQ=255)
                                                                                                                 taAACCTCCCGTGACGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAgtggtg                                     <  2:88297/150‑1 (MQ=255)
                                                                                                                  aaaCCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGa                                    >  1:1260361/1‑150 (MQ=255)
                                                                                                                  aaaCCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGa                                    >  2:651108/1‑150 (MQ=255)
                                                                                                                   aaCCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGaa                                   >  2:1281647/1‑150 (MQ=255)
                                                                                                                     ccTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACCTGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGATGAGGAGTTTGTTTGTCCAGGAAGTGGTGAAGt                                 >  2:821764/1‑150 (MQ=255)
                                                                                                                        cccGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGGTCGTTTTGGGCTAGATCGATCATGTGGGCGGGTGAATGTTGTGCGCTGAGGGCTTTGTTTGGCGAGGATGTGGTGAAGTgtc                              >  1:236278/1‑148 (MQ=255)
                                                                                                                        cccGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGAc                              >  1:99870/1‑150 (MQ=255)
                                                                                                                              cggcggcCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTAc                        <  2:438353/149‑1 (MQ=255)
                                                                                                                              aGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGTGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTAc                        >  2:437396/1‑150 (MQ=255)
                                                                                                                                gcggcCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACcg                      >  2:1033501/1‑150 (MQ=255)
                                                                                                                                   ccaaccAGCTGCATGTGCCAATGCCGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACcggcg                   >  1:170748/5‑150 (MQ=255)
                                                                                                                                     ccccAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACcggcggc                 <  1:1109247/150‑1 (MQ=255)
                                                                                                                                      accAGCTGCATGTGCCAATGCCGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACcggcggcg                >  1:1184947/2‑150 (MQ=255)
                                                                                                                                        cAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGt              >  1:1651765/1‑150 (MQ=255)
                                                                                                                                         aGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTG‑TTGGGCAGGAAGTGGTGAAGTGACCATTACCGGCGGCGGTAc            >  1:1542445/1‑150 (MQ=255)
                                                                                                                                          gCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTAc            >  1:13549/1‑150 (MQ=255)
                                                                                                                                             gCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGc         >  1:822927/1‑150 (MQ=255)
                                                                                                                                              cATTGGCACAGGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGCCGCCGGTaaaggt        <  2:575814/150‑6 (MQ=255)
                                                                                                                                                 tGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCtcat     <  1:237110/150‑1 (MQ=255)
                                                                                                                                                  ggCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCtcatc    <  1:1581058/150‑1 (MQ=255)
                                                                                                                                                    cacaTGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCAc  >  2:589067/1‑150 (MQ=255)
                                                                                                                                                    cacaTGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCAc  >  2:482033/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCAC  >  NC_000913/4072662‑4072959

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