Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,252,708 A→T 5.2% E548V (GAG→GTG)  dhaR → dhaKLM operon transcription activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,252,7080AT 5.2% 252.8 / 3.1 77E548V (GAG→GTG) dhaRdhaKLM operon transcription activator
Reads supporting (aligned to +/- strand):  ref base A (35/38);  new base T (2/2);  total (37/40)
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.91e-01

CGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTAT  >  NC_000913/1252560‑1252855
                                                                                                                                                    |                                                                                                                                                   
cgTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCga                                                                                                                                                     <  1:114147/149‑1 (MQ=255)
cgTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCga                                                                                                                                                     <  5:128665/149‑1 (MQ=255)
    gcagcaTTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGaa                                                                                                                                                  <  7:150755/148‑1 (MQ=255)
        caTTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTg                                                                                                                                             <  2:70080/149‑1 (MQ=255)
          ttCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACTCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGct                                                                                                                                           <  7:85340/149‑1 (MQ=255)
            ccGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGctc                                                                                                                                          <  7:274701/148‑1 (MQ=255)
             cGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGctct                                                                                                                                         <  4:312271/148‑1 (MQ=255)
             cGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGctct                                                                                                                                         <  5:10828/148‑1 (MQ=255)
             cGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGctct                                                                                                                                         >  7:139279/1‑148 (MQ=255)
               gcgcTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCACGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGag                                                                                                                                      >  4:181887/1‑149 (MQ=255)
                   tGGTGAATACCAAATTACGCAGGCGTGAAAAACGCTTCTCTACGCTGCTGAAAATGGATGACGATGCCCTCTCTCGCCTGGTTTCTTGGGCATGGCCCGGCAACGATTTTGACCTTTACAGCGTCATCGAGAATCTTGCTCTGagtagt                                                                                                                                  <  5:138386/149‑1 (MQ=255)
                     gTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGagtagt                                                                                                                                  <  5:452546/147‑1 (MQ=255)
                       gAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGa                                                                                                                                       >  6:126068/1‑140 (MQ=255)
                       gAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGa                                                                                                                                       <  5:126068/140‑1 (MQ=255)
                           aaCAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACgg                                                                                                                          >  8:266558/1‑149 (MQ=255)
                             cAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGgc                                                                                                                        <  8:24392/149‑1 (MQ=255)
                             cAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGgc                                                                                                                        >  1:287355/1‑149 (MQ=255)
                                   aCGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGgc                                                                                                                        <  8:190164/143‑1 (MQ=255)
                                   aCGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGgc                                                                                                                        >  7:190164/1‑143 (MQ=255)
                                       aGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGt                                                                                                              <  3:336124/149‑1 (MQ=255)
                                           ttGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCACCGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGa                                                                                                                                       >  4:326015/1‑120 (MQ=255)
                                           ttGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGa                                                                                                                                       <  3:326015/120‑1 (MQ=255)
                                             gAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCTAGAATCTAGCTCTGAGTAGTGATAACGGGCGCATACGCGGCAGTGa                                                                                                        >  6:291884/1‑149 (MQ=255)
                                             gAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGa                                                                                                        >  8:418818/1‑149 (MQ=255)
                                             gAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGGTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGt                                                                                                          <  7:59057/147‑1 (MQ=255)
                                              aaaaaCGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGAt                                                                                                       >  4:181824/1‑149 (MQ=255)
                                              aaaaaCGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGAt                                                                                                       >  6:64081/1‑149 (MQ=255)
                                                aaaCGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGATCTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGAttt                                                                                                     >  3:121554/1‑149 (MQ=255)
                                                  aCGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGc                                                                                                   >  7:444697/1‑149 (MQ=255)
                                                  aCGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGc                                                                                                   <  8:49931/149‑1 (MQ=255)
                                                   cGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGcc                                                                                                  <  4:315749/149‑1 (MQ=255)
                                                    gCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCg                                                                                                 <  8:139279/149‑1 (MQ=255)
                                                      ttCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTg                                                                                                                                 >  7:266290/1‑115 (MQ=255)
                                                      ttCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTg                                                                                                                                 <  8:266290/115‑1 (MQ=255)
                                                      ttCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGACAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCgg                                                                                                >  6:265298/1‑148 (MQ=255)
                                                        ctctACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGaa                                                                                              <  4:200586/148‑1 (MQ=255)
                                                        ctctACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGaa                                                                                              <  2:287355/148‑1 (MQ=255)
                                                                  cTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGATCTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTAc                                                                                   <  4:121554/149‑1 (MQ=255)
                                                                   tGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGTGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTAcc                                                                                  >  7:383062/1‑149 (MQ=255)
                                                                     aaaaTTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGa                                                                                <  4:454278/149‑1 (MQ=255)
                                                                     aaaaTTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGa                                                                                <  8:444697/149‑1 (MQ=255)
                                                                           gatgaCGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAgg                                                                           <  4:46657/148‑1 (MQ=255)
                                                                           gatgaCGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGc                                                                          >  1:449948/1‑149 (MQ=255)
                                                                            atgaCGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGGCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGc                                                                          <  7:224195/148‑1 (MQ=255)
                                                                             tgaCGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGAttt                                                                                                     <  4:250804/120‑1 (MQ=255)
                                                                             tgaCGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGAttt                                                                                                     >  3:250804/1‑120 (MQ=255)
                                                                               aCGATGCCCTCGCTCTCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGAc                                                                       <  5:291884/148‑1 (MQ=255)
                                                                               aCGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGttt                                                                                     <  7:326652/134‑1 (MQ=255)
                                                                               aCGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGttt                                                                                     >  8:326652/1‑134 (MQ=255)
                                                                                    gCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCTCCGATTTTGGAGGTTACAGCGTCATCGAGAATCTTGCTGGGAGTAGTGGTCACGGGCGCATTCGCGGGAGTGATTTGCCGGAACGTCTGTTTGCCGAGCAGGGGACAgaggg                                                                 >  4:243588/1‑146 (MQ=255)
                                                                                    gCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAg                                                                     <  3:243588/145‑1 (MQ=255)
                                                                                      ccTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGCTTTTGAACTTTACAGCGTCTTCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCg                                                                                 >  4:369025/1‑131 (MQ=255)
                                                                                      ccTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGatgatg                                                               >  6:438308/1‑149 (MQ=255)
                                                                                      ccTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGACAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGatgatg                                                               >  8:403813/1‑149 (MQ=255)
                                                                                       ctcgctcgCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGatgatg                                                               >  4:109007/1‑148 (MQ=255)
                                                                                       ctcgctcgCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGt                                                              >  1:278856/1‑149 (MQ=255)
                                                                                              gcCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCa                                                                             >  3:298914/1‑127 (MQ=255)
                                                                                              gcCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCa                                                                             <  4:298914/127‑1 (MQ=255)
                                                                                                 tGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGTGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCAcc                                                    >  5:160205/1‑149 (MQ=255)
                                                                                                   gTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATACCGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGAGACAGATGATGTCAGCGCCACCCg                                                  >  1:334854/1‑149 (MQ=255)
                                                                                                       cTTGTGCATGCCCAGGCAACGATTTTGACCTTTACAGCGTCATCGAGAATCTTGATCTGAGTATTGCTAACGGGCGCATTCTCGGCAGTGATTTGCCTGAACATCTGTTTAGCGAGCAGGCGACAGATGATGTCAGAGCCACACGCCt                                               <  2:77972/148‑1 (MQ=255)
                                                                                                          gtgCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCTCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTc                                            <  8:312357/148‑1 (MQ=255)
                                                                                                            gCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGATCTGAGTAGTGATAACGGGCGCATTCCCGTCAGTGATTTGCCGGAACATCTGTTTAcc                                                                                  >  5:185313/1‑108 (MQ=255)
                                                                                                                  ccAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGGGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCt                                   <  7:491580/149‑1 (MQ=255)
                                                                                                                   cAGGCAACGATTTTGAACTTTACCGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTg                                  >  2:333736/1‑149 (MQ=255)
                                                                                                                   cAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTg                                  >  4:381089/1‑149 (MQ=255)
                                                                                                                   cAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTg                                  <  3:181824/149‑1 (MQ=255)
                                                                                                                    aGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGt                                 >  7:572/1‑149 (MQ=255)
                                                                                                                    aGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGt                                 <  3:181887/149‑1 (MQ=255)
                                                                                                                        aaCGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCAtt                             >  7:498927/1‑149 (MQ=255)
                                                                                                                        aaCGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCAtt                             >  6:268294/1‑149 (MQ=255)
                                                                                                                                tGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGaa                      <  5:64081/148‑1 (MQ=255)
                                                                                                                                     tttACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGa                 >  8:134260/1‑148 (MQ=255)
                                                                                                                                     tttACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGa                 >  3:98252/1‑148 (MQ=255)
                                                                                                                                          aGCGTCATCGTGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAaga           <  8:383062/149‑1 (MQ=255)
                                                                                                                                           gCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAagag          <  5:438308/149‑1 (MQ=255)
                                                                                                                                            cGTCATCGTGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAgg         <  6:160205/149‑1 (MQ=255)
                                                                                                                                            cGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAgg         <  4:98252/149‑1 (MQ=255)
                                                                                                                                                aTCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAat     >  7:486334/1‑149 (MQ=255)
                                                                                                                                                   gagaATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTg                                  <  3:61170/117‑1 (MQ=255)
                                                                                                                                                   gagaATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTg                                  >  4:61170/1‑117 (MQ=255)
                                                                                                                                                   gagaATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAagag          >  4:160888/1‑141 (MQ=255)
                                                                                                                                                   gagaATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAagag          <  3:160888/141‑1 (MQ=255)
                                                                                                                                                   gagaATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAattat  <  8:486334/149‑1 (MQ=255)
                                                                                                                                                   gagaATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAattat  >  6:232429/1‑149 (MQ=255)
                                                                                                                                                   gagaATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAattaa  >  3:126600/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGTGGCAGCATTCCGGCGCTGGTGAATAACAAATTACGCAGTCTTGAAAAACGCTTCTCTACGCGGCTGAAAATTGATGACGATGCCCTCGCTCGCCTGGTTTCTTGTGCATGGCCAGGCAACGATTTTGAACTTTACAGCGTCATCGAGAATCTTGCTCTGAGTAGTGATAACGGGCGCATTCGCGTCAGTGATTTGCCGGAACATCTGTTTACCGAGCAGGCGACAGATGATGTCAGCGCCACCCGCCTTTCCACCAGTCTGTCATTTGCGGAAGTTGAAAAAGAGGCAATTAT  >  NC_000913/1252560‑1252855

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: