Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,250,962 T→C 5.7% intergenic (‑124/‑104) dhaK ← / → dhaR dihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit/dhaKLM operon transcription activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,250,9620TC 5.7% 304.1 / 6.1 88intergenic (‑124/‑104)dhaK/dhaRdihydroxyacetone kinase, PTS‑dependent, dihydroxyacetone‑binding subunit/dhaKLM operon transcription activator
Reads supporting (aligned to +/- strand):  ref base T (41/42);  new base C (3/2);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

ATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGTGGCGCTTTTAACAACGATGGTCGGGGCATA  >  NC_000913/1250821‑1251101
                                                                                                                                             |                                                                                                                                           
atTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGaa                                                                                                                                       <  2:750196/148‑1 (MQ=255)
 ttGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGaaa                                                                                                                                      >  1:1039525/1‑148 (MQ=255)
    aTCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACat                                                                                                                                   <  2:1039525/148‑1 (MQ=255)
         ttttttCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTg                                                                                                                                         >  2:530870/1‑137 (MQ=255)
         ttttttCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTg                                                                                                                                         <  1:530870/137‑1 (MQ=255)
                 ttattTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTa                                                                                                                           <  3:857329/143‑1 (MQ=255)
                 ttattTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTa                                                                                                                           >  4:857329/1‑143 (MQ=255)
                  tattTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTaatgaat                                                                                                                     >  4:465168/1‑148 (MQ=255)
                    ttttGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGt                                                                                                                   <  3:465168/148‑1 (MQ=255)
                             aGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTcc                                                                                                                <  2:166780/142‑1 (MQ=255)
                             aGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTcc                                                                                                                >  1:166780/1‑142 (MQ=255)
                             aGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCat                                                                                                              <  2:111505/144‑1 (MQ=255)
                             aGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCat                                                                                                              >  1:111505/1‑144 (MQ=255)
                                aaTTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCAt                                                                                                                                             >  2:632578/1‑110 (MQ=255)
                                aaTTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCAt                                                                                                                                             <  1:632578/110‑1 (MQ=255)
                                   tACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGaa                                                                                                                                       >  1:1027272/1‑113 (MQ=255)
                                   tACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGaa                                                                                                                                       <  2:1027272/113‑1 (MQ=255)
                                           ggCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACgtgt                                                                                            >  1:163492/1‑148 (MQ=255)
                                             cATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCatat                                                                                                            <  4:838740/130‑1 (MQ=255)
                                             cATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCatat                                                                                                            >  3:838740/1‑130 (MQ=255)
                                                 gatgatGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAAt                                                                                      <  2:163492/148‑1 (MQ=255)
                                                  atgatgTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACg                                                                                               >  3:607634/1‑138 (MQ=255)
                                                  atgatgTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACg                                                                                               <  4:607634/138‑1 (MQ=255)
                                                       gTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTAAATGAATGTTCCATATTGAAACTTTTACg                                                                                               >  2:796521/1‑133 (MQ=255)
                                                       gTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTAAATGAATGTTCCATATTGAAACTTTTACg                                                                                               <  1:796521/133‑1 (MQ=255)
                                                            aCGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAATCATATGCTTTAATGAAGGTTCCATATTGAATCTTTTACGTGTa                                                                                           <  4:126612/132‑1 (MQ=255)
                                                            aCGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACata                                                                                                                                  <  2:288862/93‑1 (MQ=255)
                                                            aCGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACata                                                                                                                                  >  1:288862/1‑93 (MQ=255)
                                                            aCGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTa                                                                                           >  3:126612/1‑132 (MQ=255)
                                                            aCGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCACTATGAAACat                                                                                                                                   >  2:1118420/1‑92 (MQ=255)
                                                            aCGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCACTATGAAACat                                                                                                                                   <  1:1118420/92‑1 (MQ=255)
                                                             cGACACGGCGACCACTGAGTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAGAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGtt                                                                                                                  >  4:786335/1‑108 (MQ=255)
                                                             cGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCAt                                                                                                                                             <  3:935981/81‑1 (MQ=255)
                                                             cGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGtt                                                                                                                  <  3:786335/108‑1 (MQ=255)
                                                             cGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATGATGTTTCAt                                                                                                                                             >  4:935981/1‑81 (MQ=255)
                                                                             gactgCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGaa                                                                                                                                       >  4:987045/1‑71 (MQ=255)
                                                                             gactgCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTGTCATAGTGaa                                                                                                                                       <  3:987045/71‑1 (MQ=255)
                                                                                  ccGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTa                                                                                                 >  2:972758/1‑104 (MQ=255)
                                                                                  ccGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTa                                                                                                 <  1:972758/104‑1 (MQ=255)
                                                                                    gATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACtt                                                                                  >  1:657487/1‑117 (MQ=255)
                                                                                    gATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACtt                                                                                  <  2:657487/117‑1 (MQ=255)
                                                                                      tGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTa                                                                                           <  2:450567/106‑1 (MQ=255)
                                                                                      tGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTa                                                                                           >  1:450567/1‑106 (MQ=255)
                                                                                          tCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGa                                                                  >  4:972713/1‑127 (MQ=255)
                                                                                          tCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGa                                                                  <  3:972713/127‑1 (MQ=255)
                                                                                             aTTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTc                                                                                                                 >  4:274496/1‑77 (MQ=255)
                                                                                             aTTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTc                                                                                                                 <  3:274496/77‑1 (MQ=255)
                                                                                              ttGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCat                                                                                                              <  3:779233/79‑1 (MQ=255)
                                                                                              ttGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCat                                                                                                              >  4:779233/1‑79 (MQ=255)
                                                                                              ttGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCg                                                                    >  2:873339/1‑121 (MQ=255)
                                                                                              ttGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCg                                                                    <  1:873339/121‑1 (MQ=255)
                                                                                              ttGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCa                                                     >  2:705112/1‑136 (MQ=255)
                                                                                              ttGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCa                                                     <  1:705112/136‑1 (MQ=255)
                                                                                                         gCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCg                                                                         <  1:112991/105‑1 (MQ=255)
                                                                                                         gCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCg                                                                         >  2:112991/1‑105 (MQ=255)
                                                                                                          cAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTa                                                                                 <  2:592449/96‑1 (MQ=255)
                                                                                                          cAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTa                                                                                 >  1:592449/1‑96 (MQ=255)
                                                                                                               gaaaagaaaaTTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACg                                                                                               >  4:456213/1‑77 (MQ=255)
                                                                                                               gaaaagaaaaTTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACg                                                                                               <  3:456213/77‑1 (MQ=255)
                                                                                                               gaaaagaaaaTTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGTGGCGCttt                        >  3:364773/1‑148 (MQ=255)
                                                                                                                    gaaaaTTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGa                                                                  >  2:1068663/1‑101 (MQ=255)
                                                                                                                    gaaaaTTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGa                                                                  <  1:1068663/101‑1 (MQ=255)
                                                                                                                    gaaaaTTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCggg                                        <  1:331932/127‑1 (MQ=255)
                                                                                                                    gaaaaTTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCggg                                        >  2:331932/1‑127 (MQ=255)
                                                                                                                    gaaaaTTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGTGGCGCTTTTaaca                   <  4:364773/148‑1 (MQ=255)
                                                                                                                     aaaaTTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATa                                                                                     <  1:63350/81‑1 (MQ=255)
                                                                                                                     aaaaTTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATa                                                                                     >  2:63350/1‑81 (MQ=255)
                                                                                                                     aaaaTTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATg                                   <  2:661997/131‑1 (MQ=255)
                                                                                                                     aaaaTTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATg                                   >  1:661997/1‑131 (MQ=255)
                                                                                                                               cAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTaaa                                                                               >  2:409648/1‑77 (MQ=255)
                                                                                                                               cAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTaaa                                                                               <  1:409648/77‑1 (MQ=255)
                                                                                                                               cAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGa                                                                  <  1:58019/90‑1 (MQ=255)
                                                                                                                               cAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGa                                                                  >  2:58019/1‑90 (MQ=255)
                                                                                                                               cAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCggg                                        <  2:306590/116‑1 (MQ=255)
                                                                                                                               cAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCggg                                        >  1:306590/1‑116 (MQ=255)
                                                                                                                               cAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGat                                      <  2:313488/118‑1 (MQ=255)
                                                                                                                               cAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGat                                      >  1:313488/1‑118 (MQ=255)
                                                                                                                                aaaTTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTaaa                                                                               >  3:604347/1‑76 (MQ=255)
                                                                                                                                aaaTTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAAt                                                                             >  3:144867/1‑78 (MQ=255)
                                                                                                                                aaaTTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAAt                                                                             <  4:144867/78‑1 (MQ=255)
                                                                                                                                aaaTTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGTGGCGCTTTTa                      >  3:968556/1‑133 (MQ=255)
                                                                                                                                aaaTTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGTGGCGCTTTTa                      <  4:968556/133‑1 (MQ=255)
                                                                                                                                aaaTTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGTGGCGCTTTTAACAACGa               >  3:877528/1‑140 (MQ=255)
                                                                                                                                aaaTTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGTGGCGCTTTTAACAACGa               <  4:877528/140‑1 (MQ=255)
                                                                                                                                aaaTCATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTaaa                                                                               <  4:604347/76‑1 (MQ=255)
                                                                                                                                  aTTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTa                                                                                                 <  2:52740/56‑1 (MQ=255)
                                                                                                                                  aTTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTa                                                                                                 >  1:52740/1‑56 (MQ=255)
                                                                                                                                    tATGTTTCACTATGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATa                                                   >  1:445435/1‑100 (MQ=255)
                                                                                                                                    tATGTTTCACTATGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATa                                                   <  2:445435/100‑1 (MQ=255)
                                                                                                                                     aTGTTTCACTATGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGTGGCGCTTTTAACAACGATGGTCGGGGCata  >  1:457036/1‑148 (MQ=255)
                                                                                                                                        tttCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAg                                                 >  4:59249/1‑98 (MQ=255)
                                                                                                                                        tttCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAg                                                 <  3:59249/98‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                           
ATTGATCAATTTTTTCATTATTTTGCTCCAGCAATTACGGTAGGGCATGGATGATGTTCAACGACACGGCGACCACTGACTGCCGATGAATCCATTGTGCATCAGGCAAGGGAAAAGAAAATTAGCCCAAATTATGTTTCATAGTGAAACATATGCTTTAATGAATGTTCCATATTGAAACTTTTACGTGTATTAATACTTAAAATTGCGAGCCGGAACACCTTTTGTCATAAGGGATGCGGGATATGAGTGGCGCTTTTAACAACGATGGTCGGGGCATA  >  NC_000913/1250821‑1251101

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