Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,263,518 T→A 7.3% N48I (AAT→ATT)  tdcD ← propionate kinase/acetate kinase C, anaerobic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,263,5180TA 7.3% 919.2 / 5.0 300N48I (AAT→ATT) tdcDpropionate kinase/acetate kinase C, anaerobic
Reads supporting (aligned to +/- strand):  ref base T (88/168);  new base A (28/2);  total (130/170)
Fisher's exact test for biased strand distribution p-value = 2.33e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 307 total aligned reads displayed.

CGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTT  >  NC_000913/3263369‑3263667
                                                                                                                                                     |                                                                                                                                                     
     tCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGGGCCAGCGGGGCTGGGTCTCCCCCCATTACg                                                                                                                                                  >  2:2713175/1‑150 (MQ=255)
     tCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACg                                                                                                                                                  <  2:3511893/150‑1 (MQ=255)
       ggTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCGCGGTGCTGGCTCGACCCCAGTTACGGa                                                                                                                                                >  1:3066799/1‑150 (MQ=255)
          aaaaaTACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCCTCGTAGCTGTGGTGATCCACCCGGGCTGGCTCTCCCCCAATTACGGaaaa                                                                                                                                             >  2:2266207/1‑150 (MQ=255)
          aaaaaTACTGCCGCCGTGAGCGATGCGGGGGCCAATACAGGCCACACTGACATTTAACTTCCGTTTTTCCAGTTCAACTGCAATGGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATaa                                                                                                                                             <  2:1373687/150‑1 (MQ=255)
              aTACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAAt                                                                                                                                         <  1:274967/150‑1 (MQ=255)
                  tGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGGGCCACCGGGGCGGGCGCGCCCCCCGATACGGATAAGAATGcca                                                                                                                                     >  1:614515/1‑148 (MQ=255)
                  tGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCAt                                                                                                                                     <  1:1726050/150‑1 (MQ=255)
                  tGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCAt                                                                                                                                     <  1:3782883/150‑1 (MQ=255)
                    ccgccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGTGCCAGCGGGGCTGGCTCTCCCCCATTTAAGGATAAGAATGCAatt                                                                                                                                   >  1:3041901/1‑150 (MQ=255)
                    ccgccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCAGTTACGGATAAGAATGCAttt                                                                                                                                   >  1:5867190/1‑150 (MQ=255)
                    ccgccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCGGGCTCTCCCCCATTTACGGATAAGAATGCAttt                                                                                                                                   >  2:3861187/1‑150 (MQ=255)
                      gccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTc                                                                                                                                 <  1:5207307/150‑1 (MQ=255)
                       ccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTcc                                                                                                                                <  1:5332623/150‑1 (MQ=255)
                        cgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAAGTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCg                                                                                                                               <  1:4999979/150‑1 (MQ=255)
                         gTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGa                                                                                                                              <  1:4148579/150‑1 (MQ=255)
                         gTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGa                                                                                                                              <  1:5925550/150‑1 (MQ=255)
                         gTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGa                                                                                                                              <  1:1142460/150‑1 (MQ=255)
                         gTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGa                                                                                                                              <  1:1320087/150‑1 (MQ=255)
                           gAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCACTGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGt                                                                                                                            <  2:213871/150‑1 (MQ=255)
                           gAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGt                                                                                                                            <  1:2881632/150‑1 (MQ=255)
                             gCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTa                                                                                                                          >  2:3475935/1‑150 (MQ=255)
                              cGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTaa                                                                                                                         <  1:2010099/150‑1 (MQ=255)
                                 tGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATAc                                                                                                                      <  2:5044974/150‑1 (MQ=255)
                                  gCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATAcc                                                                                                                     <  2:618027/150‑1 (MQ=255)
                                    ggTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGt                                                                                                                   <  1:694637/150‑1 (MQ=255)
                                     gTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTc                                                                                                                  <  1:5402112/150‑1 (MQ=255)
                                       ggCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCgg                                                                                                                <  2:4089205/150‑1 (MQ=255)
                                          cAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCaa                                                                                                             <  2:4417422/150‑1 (MQ=255)
                                          cAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCACCGGGGCTGGCTCACCCCCAGATACGGAAAAGAAATCAATTTCCGAGTTAATACCggcggcca                                                                                                             >  1:5863408/1‑148 (MQ=255)
                                                 ggCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGa                                                                                                      <  2:1874331/150‑1 (MQ=255)
                                                 ggCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGa                                                                                                      <  1:4321472/150‑1 (MQ=255)
                                                  gCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCGCCGTGCTGGCTCTCCCCCAATTACGGCTAAGAATGCAATTTACGAGTTAAAAACCTCGGCCAAACCAGaa                                                                                                     >  2:2816228/1‑149 (MQ=255)
                                                  gCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGAc                                                                                                     >  2:5454364/1‑150 (MQ=255)
                                                    cacaCTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACAt                                                                                                   <  1:1577354/150‑1 (MQ=255)
                                                     acacTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGTGCCAGCGGTGCTGGGTCTCCCCCATTTAAGGTTAATAATGCATTTTACGAGTTAATTACGTCGGGAATAACTGACAtt                                                                                                  >  2:4134327/1‑150 (MQ=255)
                                                       acTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGGTAAGAATGCATTTTACGAGTTAATACCGTCGGCAATACCTGACATTaa                                                                                                >  2:3238663/1‑150 (MQ=255)
                                                        cTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAAt                                                                                               <  2:5928759/150‑1 (MQ=255)
                                                        cTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAAt                                                                                               <  1:823470/150‑1 (MQ=255)
                                                         tGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATa                                                                                              <  2:5170531/150‑1 (MQ=255)
                                                            cATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACtt                                                                                           <  1:5088141/150‑1 (MQ=255)
                                                              tttAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGTGCCAGCGGTGCTGGCTCTCCCCCATTTACGGGTAAGAATGCAATTTCCCAGGTAATACCGTCGGCAAAACCTGACAATAAAACATca                                                                                         >  2:2335008/1‑150 (MQ=255)
                                                              tttAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTca                                                                                         <  1:3810272/150‑1 (MQ=255)
                                                              tttAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTca                                                                                         <  1:1775679/150‑1 (MQ=255)
                                                              tttAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTca                                                                                         <  2:5984960/150‑1 (MQ=255)
                                                              tttAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTca                                                                                         <  1:4464587/150‑1 (MQ=255)
                                                               ttAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCAATTTCCGAGTTAATACCGTCGGCAATAACTGACATTAATACTTaac                                                                                        >  2:4333917/1‑150 (MQ=255)
                                                                tAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTcaca                                                                                       <  2:2616402/150‑1 (MQ=255)
                                                                 aaaTTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGTGCCAGCGGTGCTGGCTCGCCCCCGTTTACGGGTGAGACTTCATTTTCCGAGTTAAAACCGTCGGCAAAACCTGACCTTAATAcatgaccg                                                                                      >  2:4190433/1‑143 (MQ=255)
                                                                  aaTTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTAAGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATAACTGACATTAATACATcacact                                                                                     >  2:2868445/1‑148 (MQ=255)
                                                                  aaTTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCACCGGTGCTGGCTCTCCCCCAATTACGGATAAGAATTCAATTTACGAGTTAATAACCGCGGCAAAAACAGAAATTAAAAATTcacata                                                                                     >  1:3706068/1‑148 (MQ=255)
                                                                     tCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCACCGGTGCAGGCGCACCCCCAAATAAGGAGAACAAGCCAATTTAAGAGTTAATAAACGCCGCAAAAACAGACAATTATAAATCACAGTCGc                                                                                  >  2:1027869/1‑150 (MQ=255)
                                                                      ccGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCt                                                                                 <  2:5531685/150‑1 (MQ=255)
                                                                      ccGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCt                                                                                 <  1:1153749/150‑1 (MQ=255)
                                                                       cGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCAATTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTg                                                                                <  1:1437272/150‑1 (MQ=255)
                                                                            ttCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATc                                                                           <  1:4269150/150‑1 (MQ=255)
                                                                               cAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAg                                                                        <  2:4580089/150‑1 (MQ=255)
                                                                                  ttCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTGCCGAGTTAATACCGGCGGCCAGACCGGACAGTAATACATCACCGGCCCGGGCCTCGAGccc                                                                     >  2:4493387/1‑148 (MQ=255)
                                                                                      aaTGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGaa                                                                 <  1:3214260/150‑1 (MQ=255)
                                                                                      aaTGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGaa                                                                 <  1:2240080/150‑1 (MQ=255)
                                                                                      aaTGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGaa                                                                 <  1:5790986/150‑1 (MQ=255)
                                                                                          cAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGTGCCAGCGGTGCTGGCTCTCCCCCATTTACCGATAAGAAGGCATTTTCCGAGTTAATAACGTCGGCAATAACTGACATTAATAATTCACAGTCGCTGGCCACGAGCACGGAAAAca                                                             >  2:3636738/1‑149 (MQ=255)
                                                                                            aTTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTa                                                           <  1:3083438/150‑1 (MQ=255)
                                                                                            aTTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTa                                                           <  1:1285335/150‑1 (MQ=255)
                                                                                             ttGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCACCGGTGCTGGCTCTCCCCCAGTTAAGGATAAGAATGCAATTTACGAGTTAATACCCTCGGCAAAACCAGACATTAATACATCACAGTAGCTGGCAACGAGCCCGGAAAACATaa                                                          >  2:3111555/1‑150 (MQ=255)
                                                                                               gCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATc                                                        <  1:2298753/150‑1 (MQ=255)
                                                                                               gCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATc                                                        <  2:3044506/150‑1 (MQ=255)
                                                                                                  ttCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCCATTACGGATAAGTATGCAATTTCCGAGTTAATACCGTCGGCAATAACTGACATTAATAATTCACAGTAGCTGGCATAGTGCACGGGAAACTTAAtagga                                                     >  2:5592619/1‑146 (MQ=255)
                                                                                                  ttCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGaa                                                     <  2:2893910/150‑1 (MQ=255)
                                                                                                  ttCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGaa                                                     <  1:4810901/150‑1 (MQ=255)
                                                                                                  ttCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATAACTGACATTAATAATTCACAGTCGCTGGCCTCGAGCACGGAAAACTTAATCCaa                                                     >  2:4795428/1‑150 (MQ=255)
                                                                                                  ttCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACCGATAAGAATTCATTTTACGAGTTAAAAACCTCCGCCAAACCAGAAAATAAAAAATAACACACGCTGGCCACGCGCACCGAAAACTTaaaagaa                                                     >  2:5669030/1‑145 (MQ=255)
                                                                                                        gCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCAATTACGGATAAGAATGCATTTTACGAGTTAATACCGTCGGCAATACCTGACATTAATACATCACAGTCGCTGGCATCGAGCACGGGAAACTTAATCGAAGACGaa                                               >  2:3738916/1‑150 (MQ=255)
                                                                                                                 tAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCTTTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTta                                      <  1:575886/150‑1 (MQ=255)
                                                                                                                 tAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTta                                      <  2:1256167/150‑1 (MQ=255)
                                                                                                                   gCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTtaat                                    <  1:4240112/150‑1 (MQ=255)
                                                                                                                       tggtgAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAcc                                >  2:4953491/1‑150 (MQ=255)
                                                                                                                       tggtgAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAcc                                >  2:377311/1‑150 (MQ=255)
                                                                                                                        ggagAGCCAGCACCGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCa                               >  2:5298866/4‑150 (MQ=255)
                                                                                                                           gAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaaaa                            <  1:4875744/150‑1 (MQ=255)
                                                                                                                           gAGCCAGCACCGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaaaa                            >  1:4392127/1‑150 (MQ=255)
                                                                                                                              ccAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAaacaa                         >  1:4008903/1‑150 (MQ=255)
                                                                                                                              ccAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAaacaa                         >  1:1960458/1‑150 (MQ=255)
                                                                                                                              ccAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAaacaa                         >  2:4736858/1‑150 (MQ=255)
                                                                                                                              ccAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAaacaa                         >  2:1700278/1‑150 (MQ=255)
                                                                                                                               cAGCGGTGCTGGCTCTCCCCCATTTACGTATAAGATTGCATTTTCCGAGTTAATACCGTCGGCAATACCTGAAATTAATACTTCACAGTCGCTGGCATCGAGCACGTAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAaacaac                        <  1:939355/150‑1 (MQ=255)
                                                                                                                                aGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAAcc                       >  1:3463144/1‑150 (MQ=255)
                                                                                                                                 gCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCg                      <  2:1487776/150‑1 (MQ=255)
                                                                                                                                  cGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCgg                     >  1:2570028/1‑150 (MQ=255)
                                                                                                                                    gTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATAACGTCGGCAATACCTGACAATAATACTTCACAGTCGCTTGCATCGAGCACGGAAAACTTAAAAGAAGACGAACCACAGTTAATAACCAAAAAAACCGGaa                   >  2:1233675/1‑150 (MQ=255)
                                                                                                                                        tGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAattc               >  2:3400253/1‑150 (MQ=255)
                                                                                                                                        tGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAattc               >  2:775678/1‑150 (MQ=255)
                                                                                                                                              tcCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATtcattcat         <  1:1033112/150‑1 (MQ=255)
                                                                                                                                                ccccATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCAtct       >  1:3398990/1‑150 (MQ=255)
                                                                                                                                                 cccATTTACGGATAAGAGTGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCAtctc      >  1:3909197/1‑150 (MQ=255)
                                                                                                                                                 cccATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCAtctc      >  1:385266/1‑150 (MQ=255)
                                                                                                                                                  ccATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACCGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACTAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATctct     <  1:4531531/150‑1 (MQ=255)
                                                                                                                                                  acATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATctct     >  2:3178162/2‑150 (MQ=255)
                                                                                                                                                   cATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCtt    <  2:1207814/150‑1 (MQ=255)
                                                                                                                                                     tttACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCtttt  <  1:2138030/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCCATTTACGGATAAGAATGCATTTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAATTCATTCATCTCTTTT  >  NC_000913/3263369‑3263667

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