Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,072,802 G→T 5.8% P257Q (CCG→CAG)  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,8020GT 5.8% 317.5 / 7.4 112P257Q (CCG→CAG) yihU3‑sulpholactaldehyde (SLA) reductase, NADH‑dependent; gamma‑hydroxybutyrate dehydrogenase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base G (55/50);  new base T (2/5);  total (57/55)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.19e-01

Only 100 of 112 total aligned reads displayed.

ATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTAC  >  NC_000913/4072658‑4072949
                                                                                                                                                |                                                                                                                                                   
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)
   ccTTAATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGc                                                                                                                                             <  2:1162764/150‑1 (MQ=255)
    gtttATCCAGTCGTTTAATTTTTTCTTTGGCCGTCATCCCCGCACTGACACGGACCGGTTCCAGAATGGCGGACCAGTCCTGGGGGCCGCGAACCCAGGCGCGCGCCTGGGATTAAACCGACCGTGCGGCGGTACCCAGCGGGAATGGca                                                                                                                                            >  2:183845/2‑150 (MQ=255)
        aTCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATg                                                                                                                                        >  1:1463598/1‑150 (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)
                  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: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)
                                    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)
                                                          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:64275/150‑1 (MQ=255)
                                                                          ggACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCggg                                                                      <  2:59063/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)
                                                                                    tGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGCTCGATCATGAAGGCGGGAGAAAGAtct                                                            >  1:570231/1‑150 (MQ=255)
                                                                                       cgaccgcgaccCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCg                                                         <  1:919813/150‑1 (MQ=255)
                                                                                         accgcgaccCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATGATGAAGGCGGGAGAAAGATCTCCGcg                                                       >  1:177188/1‑149 (MQ=255)
                                                                                            gcgaccCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAg                                                    >  1:1280642/1‑150 (MQ=255)
                                                                                            gcgaccCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAg                                                    >  1:679181/1‑150 (MQ=255)
                                                                                            gcgaccCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAg                                                    >  1:1704599/1‑150 (MQ=255)
                                                                                               accCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCTGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGTCAAGAGCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGAc                                                 <  1:1249144/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)
                                                                                                                   tataAACCTGCCGTGGGGGGGCCCGCAGCGGCATTGGCACATGCAGCTGGGTGGCGACAGCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGATGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAgtgg                             <  1:217379/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)
                                                                                                                         ccTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACCTGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGATGAGGAGTTTGTTTGTCCAGGAAGTGGTGAAGt                       >  2:821764/1‑150 (MQ=255)
                                                                                                                          cTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTg                      <  2:625902/150‑1 (MQ=255)
                                                                                                                            cccGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGGTCGTTTTGGGCTAGATCGATCATGTGGGCGGGTGAATGTTGTGCGCTGAGGGCTTTGTTTGGCGAGGATGTGGTGAAGTgtc                    >  1:236278/1‑148 (MQ=255)
                                                                                                                             ccGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGAcc                   <  1:1690443/150‑1 (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)
                                                                                                                                        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)
                                                                                                                                                |                                                                                                                                                   
ATTCCTTTATCCAGTCGTTACATTTTTACTTTGGCAGTCATCCCGGCACTGACACGGACCTGTTCCAGAATGGCGGACCAGTCCTGGCGACCGCGACCCGCTGCGCGCGCCTGGCTATAAACCTCCCGTGAGGCGGCCCCCAGCGGCATTGGCACATGCAGCTGGTTGGCGACATCAAGGGCGATGCCAAGATCCTTATGGGCAAGATCGATCATGAAGGCGGGAGAAAGATCTCCGCTGAGGACTTTGTTTGGCCAGGAAGTGGTGAAGTGACCTTTACCGGCGGCGGTAC  >  NC_000913/4072658‑4072949

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