Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,248,064 T→A 5.4% K351* (AAA→TAA)  dhaM ← putative dihydroxyacetone‑specific PTS enzymes: HPr, EI components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,248,0640TA 5.4% 256.2 / 4.7 74K351* (AAA→TAA) dhaMputative dihydroxyacetone‑specific PTS enzymes: HPr, EI components
Reads supporting (aligned to +/- strand):  ref base T (33/37);  new base A (2/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

TAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTGCTTTC  >  NC_000913/1247918‑1248212
                                                                                                                                                  |                                                                                                                                                    
tAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTaa                                                                                                                                                    <  1:408131/149‑1 (MQ=255)
tAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTaa                                                                                                                                                    <  3:70945/149‑1 (MQ=255)
  gAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGa                                                                                                                                                  >  2:738258/1‑149 (MQ=255)
  gAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGa                                                                                                                                                  >  1:216754/1‑149 (MQ=255)
  gAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGa                                                                                                                                                  >  2:510481/1‑149 (MQ=255)
   aaTAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGATTGGGTCAGGTGGACCAGGGTTCGATGCAGAAGATCGTCCACACCAGTATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGac                                                                                                                                                 <  2:291100/149‑2 (MQ=255)
    aTAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAAc                                                                                                                                                <  3:693419/149‑1 (MQ=255)
       gTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTtg                                                                                                                                             <  3:382295/149‑1 (MQ=255)
         tGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTtgct                                                                                                                                           <  2:620416/149‑1 (MQ=255)
           gCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTAAAGAACttagcca                                                                                                                                         <  4:225998/149‑6 (MQ=255)
             gAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCa                                                                                                                                       >  2:704949/1‑149 (MQ=255)
             gAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCa                                                                                                                                       >  2:410731/1‑149 (MQ=255)
               gTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAgg                                                                                                                                     >  1:151718/1‑149 (MQ=255)
                ttAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCttt                                                                                                                                                      >  4:499363/1‑131 (MQ=255)
                ttAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCttt                                                                                                                                                      <  3:499363/131‑1 (MQ=255)
                     cTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATAtt                                                                                                                               <  4:295880/149‑1 (MQ=255)
                          ggAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGcc                                                                                                                          <  2:151718/149‑1 (MQ=255)
                                 ttcttTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAAt                                                                                                                   <  1:704949/149‑1 (MQ=255)
                                 ttcttTCGATTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAAt                                                                                                                   <  2:634225/149‑1 (MQ=255)
                                    tttCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGtt                                                                                                                >  3:129628/1‑149 (MQ=255)
                                    tttCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTc                                                                                                               <  4:257808/149‑1 (MQ=255)
                                    tttCGTTTGGGTCAGGTGGACCACGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGtt                                                                                                                <  1:506910/149‑1 (MQ=255)
                                      tCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCa                                                                                                              <  2:143096/149‑1 (MQ=255)
                                      tCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCa                                                                                                              <  2:170278/149‑1 (MQ=255)
                                             ggTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGaa                                                                                                       >  3:421952/1‑149 (MQ=255)
                                              gTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAg                                                                                                      >  4:563243/1‑149 (MQ=255)
                                              gTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAg                                                                                                      <  1:738258/149‑1 (MQ=255)
                                              gTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTGCATGCTGAAg                                                                                                      >  3:224088/1‑149 (MQ=255)
                                                cAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGa                                                                                                    <  2:7813/149‑1 (MQ=255)
                                                          aGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGATATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTgc                                                                                          >  2:725900/1‑149 (MQ=255)
                                                           gggTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTgcc                                                                                         <  3:296791/149‑1 (MQ=255)
                                                                cGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTAAAGAACtta                                                                                                                                             >  1:287459/1‑91 (MQ=255)
                                                                cGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTAAAGAACtta                                                                                                                                             <  2:287459/92‑2 (MQ=255)
                                                                 gATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTgccgccgcc                                                                                   >  1:83614/1‑149 (MQ=255)
                                                                     cAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCcagca                                                                               >  2:66212/1‑149 (MQ=255)
                                                                       gaagaTCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGt                                                                             >  1:141099/1‑149 (MQ=255)
                                                                        aagaTCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGtt                                                                            >  4:44470/1‑149 (MQ=255)
                                                                             cGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGGTCCGGa                                                                       <  4:421952/149‑1 (MQ=255)
                                                                              gTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGat                                                                      >  1:545271/1‑149 (MQ=255)
                                                                               tCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGatc                                                                     <  1:510481/149‑1 (MQ=255)
                                                                               tCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGatc                                                                     >  3:248908/1‑149 (MQ=255)
                                                                                ccACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGatca                                                                    <  2:37506/149‑1 (MQ=255)
                                                                                  acaTCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGatcatc                                                                  <  1:712558/149‑1 (MQ=255)
                                                                                     tCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTaa                                                               <  4:248908/149‑1 (MQ=255)
                                                                                         tataGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAgt                                                           <  4:224088/149‑1 (MQ=255)
                                                                                           taGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAgtgt                                                         <  1:13550/149‑1 (MQ=255)
                                                                                           taGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCCTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAgtgt                                                         >  4:47929/1‑149 (MQ=255)
                                                                                            aGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTa                                                        <  3:180104/149‑1 (MQ=255)
                                                                                            aGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTa                                                        <  2:580223/149‑1 (MQ=255)
                                                                                               gAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATgg                                                     >  4:407316/1‑149 (MQ=255)
                                                                                               gAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATACTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATgg                                                     >  4:564143/1‑149 (MQ=255)
                                                                                                aGCTTGTAGATATTCATCATCCAGTTGCTGATATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGt                                                    <  1:725900/149‑1 (MQ=255)
                                                                                                   ttGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGAc                                                 <  3:411962/149‑1 (MQ=255)
                                                                                                        gATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGaa                                            >  2:420266/1‑149 (MQ=255)
                                                                                                            ttCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACa                                                             <  3:406594/128‑1 (MQ=255)
                                                                                                            ttCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACa                                                             >  4:406594/1‑128 (MQ=255)
                                                                                                                tcatcCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGc                                    >  4:551909/1‑149 (MQ=255)
                                                                                                                tcatcCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGc                                    >  1:234300/1‑149 (MQ=255)
                                                                                                                tcatcCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGc                                    <  3:551909/149‑1 (MQ=255)
                                                                                                                 catcCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCg                                   <  2:141099/149‑1 (MQ=255)
                                                                                                                 catcCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCg                                   <  1:53569/149‑1 (MQ=255)
                                                                                                                    ccAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGAt                                       >  3:728162/1‑142 (MQ=255)
                                                                                                                    ccAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGAt                                       <  4:728162/142‑1 (MQ=255)
                                                                                                                     cAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCaa                               <  2:116349/149‑1 (MQ=255)
                                                                                                                               tATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAg                     >  4:735623/1‑149 (MQ=255)
                                                                                                                                     tGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCt               >  2:763720/1‑149 (MQ=255)
                                                                                                                                      ggCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTg              >  1:662192/1‑149 (MQ=255)
                                                                                                                                         tAAGTTCTTTAGGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCAt                                                                                                             >  3:490132/1‑51 (MQ=255)
                                                                                                                                         tAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCAt                                                                                                             <  4:490132/51‑1 (MQ=255)
                                                                                                                                         tAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGc            <  4:263303/148‑1 (MQ=255)
                                                                                                                                         tAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGc            >  3:263303/1‑148 (MQ=255)
                                                                                                                                            gTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCt        >  1:74139/1‑149 (MQ=255)
                                                                                                                                            gTTCTTAAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCt        >  2:601791/1‑149 (MQ=255)
                                                                                                                                             ttCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATc                           <  2:372604/129‑1 (MQ=255)
                                                                                                                                             ttCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATc                           >  1:372604/1‑129 (MQ=255)
                                                                                                                                             ttCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTg       >  3:71395/1‑149 (MQ=255)
                                                                                                                                                tttAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTGCtt    >  1:111658/1‑149 (MQ=255)
                                                                                                                                                tttAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTGCtt    >  3:191022/1‑149 (MQ=255)
                                                                                                                                                 ttAAGACCTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGTGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTGCttt   <  3:564143/149‑1 (MQ=255)
                                                                                                                                                  tAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTGCTTTc  <  3:564110/149‑1 (MQ=255)
                                                                                                                                                  tAAGAACTTGCTGCCAGGCATATTCTGACGTGCAATGTTCATGCTGAAGAAGTTCGCTTGCCGCCGCCAGCACTTCCGGATCATCTAACAGTGTCTGGTGACCAGAAAAGATTGCGGCAAGATCGTCAAGCCCGCTGGCTTCTGCTTTc  <  2:692323/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TAGAATAGTTGGCGAGTTAAACTGCGGGAGTTCTTCTTTCGTTTGGGTCAGGTGGACCAGGGTGCGATGCAGAAGATCGTCCACATCAATATAGCGAGCTTGTAGATATTCATCATCCAGTTGCTGGTATTGCTGGCTAAGTTCTTTAAGAACTTGCTGCCAGGCATATTCTGCCGTGCAATGTTCATGCTGAAGGAGTTCGCTTGCCGCCGCCAGCAGTTCCGGATCATCTAACAGTGTATGGTGACCAGAAAAGATTGCGGCAATATCGTCAAGCCCGCTGGCTTCTGCTTTC  >  NC_000913/1247918‑1248212

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: