Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,072,817 A→C 5.2% L252R (CTG→CGG)  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,8170AC 5.2% 324.1 / 7.3 117L252R (CTG→CGG) yihU3‑sulpholactaldehyde (SLA) reductase, NADH‑dependent; gamma‑hydroxybutyrate dehydrogenase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base A (56/55);  new base C (2/4);  total (58/59)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

Only 100 of 119 total aligned reads displayed.

CCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACG  >  NC_000913/4072668‑4072966
                                                                                                                                                     |                                                                                                                                                     
ccAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCa                                                                                                                                                       >  1:251386/1‑150 (MQ=255)
  cgTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGc                                                                                                                                                     <  1:899318/149‑1 (MQ=255)
   gTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCt                                                                                                                                                    <  1:1670106/150‑1 (MQ=255)
   gTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCt                                                                                                                                                    <  1:21266/150‑1 (MQ=255)
      gTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTCCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGt                                                                                                                                                 <  2:669550/150‑1 (MQ=255)
       ttACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGtt                                                                                                                                                >  1:794415/1‑150 (MQ=255)
       ttACATTTTTACTTTGGCAGTCATACCGGGACTGACACGGACCTGTTCGAGAATGGCGGACCAGTCCTGGCGACCGCGTCCCGCTGCGCGCGCCTGGCTATAATCCTCCCGTGAGGCGGCCCCCAGCTGCATGTGCCAATGCCGCTgggg                                                                                                                                                <  2:778857/150‑3 (MQ=255)
        tACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCTGCATGTGCCAATGCCGCTggggg                                                                                                                                               <  1:1034405/150‑4 (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)
                tACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATc                                                                                                                                       >  2:1463227/1‑150 (MQ=255)
                     tGGCAGTCATCCCGGCACTTACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAggg                                                                                                                                  <  2:385325/150‑1 (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)
                          gTCATCCCGGCACTGACACGGACCTATTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATg                                                                                                                             >  1:381101/1‑150 (MQ=255)
                          ctCATCCCGGCACGGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCTCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGATGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATg                                                                                                                             <  1:1513588/149‑1 (MQ=255)
                            cATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGCCCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGcc                                                                                                                           >  2:429960/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)
                                                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)
                                                           aTGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAg                                                                                            >  2:702347/1‑150 (MQ=255)
                                                                ggACCAGTCCTGGCGAGCGGGACGCGCTGCGCGCGCCTGGCTATAATCCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCggg                                                                                       <  2:135172/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)
                                                                      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)
                                                                                       ccGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACtt                                                                <  2:1463598/150‑1 (MQ=255)
                                                                                        cGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACttt                                                               <  2:381101/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)
                                                                                                gcgcCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGcc                                                       <  1:611981/150‑1 (MQ=255)
                                                                                                   ccTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAgg                                                    <  2:1505532/150‑1 (MQ=255)
                                                                                                   ccTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGAAAGATCGATAATGAAGGGGGTAGAAAGATCTCCGATGAGGAATTTGTTTGGCAAgg                                                    >  2:933205/1‑150 (MQ=255)
                                                                                                      ggCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAg                                                 >  2:1201798/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                                           >  2:651108/1‑150 (MQ=255)
                                                                                                            aaaCCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGa                                           >  1:1260361/1‑150 (MQ=255)
                                                                                                             aaCCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGaa                                          >  2:1281647/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)
                                                                                                                              cccccAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACcggcgg                         >  2:782524/1‑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)
                                                                                                                                         aTTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCtc              <  1:1263812/150‑1 (MQ=255)
                                                                                                                                           tGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCtcat            <  1:237110/150‑1 (MQ=255)
                                                                                                                                            ggCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCtcatc           <  1:1581058/150‑1 (MQ=255)
                                                                                                                                              cacaTGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCAc         >  2:482033/1‑150 (MQ=255)
                                                                                                                                              cacaTGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCAc         >  2:589067/1‑150 (MQ=255)
                                                                                                                                                caTGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACtt       >  1:1087033/1‑150 (MQ=255)
                                                                                                                                                 aTGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACttt      <  2:679181/150‑1 (MQ=255)
                                                                                                                                                  tGCAGCGGGGTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTg     <  2:251386/150‑1 (MQ=255)
                                                                                                                                                    cAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGAc   >  1:356846/1‑150 (MQ=255)
                                                                                                                                                    cAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGAc   >  2:1458054/1‑150 (MQ=255)
                                                                                                                                                    cAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGGTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGAc   >  2:1637263/1‑150 (MQ=255)
                                                                                                                                                    cAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTAATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCAGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATAACTTTGAc   >  2:398086/1‑150 (MQ=255)
                                                                                                                                                     aGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACg  <  2:479119/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTACCGCTCATCACTTTGACG  >  NC_000913/4072668‑4072966

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