Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,072,807 T→G 5.4% P255P (CCA→CCC)  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,8070TG 5.4% 326.6 / 6.6 112P255P (CCA→CCC) yihU3‑sulpholactaldehyde (SLA) reductase, NADH‑dependent; gamma‑hydroxybutyrate dehydrogenase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base T (53/53);  new base G (2/4);  total (55/57)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

Only 100 of 116 total aligned reads displayed.

ATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCAT  >  NC_000913/4072658‑4072956
                                                                                                                                                     |                                                                                                                                                     
aTTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCAtt                                                                                                                                                       >  1:59063/1‑150 (MQ=255)
aTTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCAtt                                                                                                                                                       >  1:64275/1‑150 (MQ=255)
aTTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCAtt                                                                                                                                                       >  2:140358/1‑150 (MQ=255)
  tCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTgg                                                                                                                                                     >  2:1109247/1‑150 (MQ=255)
    gtttATCCAGTCGTTTAATTTTTTCTTTGGCCGTCATCCCCGCACTGACACGGACCGGTTCCAGAATGGCGGACCAGTCCTGGGGGCCGCGAACCCAGGCGCGCGCCTGGGATTAAACCGACCGTGCGGCGGTACCCAGCGGGAATGGca                                                                                                                                                   >  2:183845/2‑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:1593687/150‑4 (MQ=255)
                  tACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCTGCATGTGCCAATGCCGCTggggg                                                                                                                                     <  1:1034405/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:1158818/150‑1 (MQ=255)
                                 gCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCg                                                                                                                      >  1:289967/1‑150 (MQ=255)
                                  cAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAACCCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGa                                                                                                                     >  2:1581058/1‑150 (MQ=255)
                                  cAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGa                                                                                                                     >  1:1481775/1‑150 (MQ=255)
                                    ttCATCCCGGCACTGACACGGCCCTGTTCCAGAATGGCGGACCAGTCCGGGCGACCGCGACCCGCCGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATg                                                                                                                   <  1:869606/149‑1 (MQ=255)
                                    gTCATCCCGGCACTGACACGGACCTATTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATg                                                                                                                   >  1:381101/1‑150 (MQ=255)
                                    ctCATCCCGGCACGGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCTCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGATGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATg                                                                                                                   <  1:1513588/149‑1 (MQ=255)
                                      cATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGCCCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGcc                                                                                                                 >  2:429960/1‑150 (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)
                                                      ggtccgtGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAAACTCCCGTGGGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGCTCCTTATGGGCa                                                                                                 >  2:217379/7‑150 (MQ=255)
                                                      cGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCCCATGCAGCTGGTTGGCGCCATCCAGGGCGATGCCAAGATCCTTATGGGCa                                                                                                 >  1:83926/1‑150 (MQ=255)
                                                        gACCTGTTCCAGAACGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCGGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAg                                                                                               <  1:1629206/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:753290/150‑1 (MQ=255)
                                                                          ggACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCggg                                                                             <  2:64275/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)
                                                                                         accgcgaccCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATGATGAAGGCGGGAGAAAGATCTCCGcg                                                              >  1:177188/1‑149 (MQ=255)
                                                                                            gcgaccCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAg                                                           >  1:1280642/1‑150 (MQ=255)
                                                                                            gcgaccCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAg                                                           >  1:1704599/1‑150 (MQ=255)
                                                                                               accCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCTGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGTCAAGAGCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGAc                                                        <  1:1249144/150‑1 (MQ=255)
                                                                                                 ccGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGTCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACtt                                                      <  2:46115/150‑1 (MQ=255)
                                                                                                 ccGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACtt                                                      <  2:1463598/150‑1 (MQ=255)
                                                                                                  cGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACttt                                                     <  1:140358/150‑1 (MQ=255)
                                                                                                  cGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACttt                                                     <  2:381101/150‑1 (MQ=255)
                                                                                                     tGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACtttgtt                                                  <  2:1177542/150‑1 (MQ=255)
                                                                                                        gcgcgcCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTgg                                               <  1:572886/150‑1 (MQ=255)
                                                                                                          gcgcCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGcc                                             <  1:611981/150‑1 (MQ=255)
                                                                                                             ccTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAgg                                          >  1:564783/1‑150 (MQ=255)
                                                                                                             ccTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGAAAGATCGATAATGAAGGGGGTAGAAAGATCTCCGATGAGGAATTTGTTTGGCAAgg                                          >  2:933205/1‑150 (MQ=255)
                                                                                                                ggCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAg                                       >  2:1201798/1‑150 (MQ=255)
                                                                                                                 gCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAgt                                      <  2:945316/150‑1 (MQ=255)
                                                                                                                   tataAACCTGCCGTGGGGGGGCCCGCAGCGGCATTGGCACATGCAGCTGGGTGGCGACAGCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGATGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAgtgg                                    <  1:217379/150‑1 (MQ=255)
                                                                                                                     taAACCTCCCGTGAGGCTGCCCCCAGCGGCATTGGCACATGCATCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAgtggtg                                  <  2:1013895/150‑1 (MQ=255)
                                                                                                                     taAACCTCCCGTGACGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAgtggtg                                  <  2:88297/150‑1 (MQ=255)
                                                                                                                         ccTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACCTGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGATGAGGAGTTTGTTTGTCCAGGAAGTGGTGAAGt                              >  2:821764/1‑150 (MQ=255)
                                                                                                                          cTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTg                             <  2:625902/150‑1 (MQ=255)
                                                                                                                            cccGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGGTCGTTTTGGGCTAGATCGATCATGTGGGCGGGTGAATGTTGTGCGCTGAGGGCTTTGTTTGGCGAGGATGTGGTGAAGTgtc                           >  1:236278/1‑148 (MQ=255)
                                                                                                                            cccGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGAc                           >  1:99870/1‑150 (MQ=255)
                                                                                                                             ccGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGAcc                          <  1:1690443/150‑1 (MQ=255)
                                                                                                                                  cggcggcCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTAc                     <  2:438353/149‑1 (MQ=255)
                                                                                                                                    gcggcCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACcg                   >  2:1033501/1‑150 (MQ=255)
                                                                                                                                       ccaaccAGCTGCATGTGCCAATGCCGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACcggcg                >  1:170748/5‑150 (MQ=255)
                                                                                                                                        cccccAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACcggcgg               >  2:782524/1‑150 (MQ=255)
                                                                                                                                         ccccAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACcggcggc              <  1:1109247/150‑1 (MQ=255)
                                                                                                                                          accAGCTGCATGTGCCAATGCCGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACcggcggcg             >  1:1184947/2‑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)
                                                                                                                                                   aTTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCtc    <  1:1263812/150‑1 (MQ=255)
                                                                                                                                                     tGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCtcat  <  1:237110/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCAT  >  NC_000913/4072658‑4072956

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: ‑