Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 721,009 T→A 6.0% intergenic (+169/+47) ybfK → / ← kdpE uncharacterized protein/response regulator in two‑component regulatory system with KdpD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913721,0090TA 6.0% 356.1 / 8.3 129intergenic (+169/+47)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base T (51/70);  new base A (3/5);  total (54/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

Only 100 of 139 total aligned reads displayed.

TTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCATATAAATA  >  NC_000913/720861‑721155
                                                                                                                                                    |                                                                                                                                                      
ttCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                                                                                                                                        <  2:510242/149‑1 (MQ=255)
ttCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                                                                                                                                        <  1:1234782/149‑1 (MQ=255)
 tCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttatta                                                                                                                                                       <  2:1540364/149‑1 (MQ=255)
     tCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTAc                                                                                                                                                   <  1:1032605/149‑1 (MQ=255)
     tCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTAc                                                                                                                                                   >  2:428833/1‑149 (MQ=255)
       cAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACgg                                                                                                                                                 >  2:363251/1‑149 (MQ=255)
        aGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                                                                                >  2:1546808/1‑149 (MQ=255)
        aGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                                                                                >  1:390106/1‑149 (MQ=255)
        aGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                                                                                >  1:110386/1‑149 (MQ=255)
          gAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt‑‑‑‑t                                                                                                                                          <  2:1406882/149‑1 (MQ=255)
          gAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt‑‑‑‑t                                                                                                                                          >  2:429880/1‑149 (MQ=255)
           aaaTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt‑‑‑‑tt                                                                                                                                         >  1:195232/1‑149 (MQ=255)
           aaaTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCt‑‑‑‑tt                                                                                                                                         <  2:1317706/149‑1 (MQ=255)
             aTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑ttaa                                                                                                                                       >  2:926547/1‑149 (MQ=255)
               aCGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑ttaatt                                                                                                                                     <  2:584851/149‑1 (MQ=255)
               aCGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCt‑‑‑‑ttttat                                                                                                                                     <  2:1136943/149‑5 (MQ=255)
                cGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑ttaatta                                                                                                                                    <  1:3274/149‑1 (MQ=255)
                   ttGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑Ttaattaat                                                                                                                                  <  2:1284893/148‑1 (MQ=255)
                   ttGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑Ttaattaat                                                                                                                                  >  1:1284893/1‑148 (MQ=255)
                      aaCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATaaaa                                                                                                                              >  2:1584737/1‑149 (MQ=255)
                      aaCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATaaaa                                                                                                                              <  1:303396/149‑1 (MQ=255)
                      aaCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATaaaa                                                                                                                              <  2:1209797/149‑1 (MQ=255)
                      aaCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATaaaa                                                                                                                              >  2:1360298/1‑149 (MQ=255)
                        cTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAagg                                                                                                                            >  2:864762/1‑149 (MQ=255)
                        cTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGATTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCTTATTTATTATTACGGAT‑‑‑‑TTAATTAATAAAaga                                                                                                                            >  2:1059660/1‑148 (MQ=255)
                         tgtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAaggc                                                                                                                           >  2:669737/1‑149 (MQ=255)
                          gtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAaggca                                                                                                                          <  2:1576737/149‑1 (MQ=255)
                           tgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAaggcag                                                                                                                         >  1:1525770/1‑149 (MQ=255)
                              aGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCt                                                                                                                      <  1:717195/149‑1 (MQ=255)
                                 tGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTa                                                                                                                   <  1:647078/149‑1 (MQ=255)
                                 tGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTa                                                                                                                   <  2:1279436/149‑1 (MQ=255)
                                  taTCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTAt                                                                                                                  <  1:405879/148‑1 (MQ=255)
                                  gATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTAt                                                                                                                  <  2:1478237/149‑1 (MQ=255)
                                    tCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTa                                                                                                                >  2:474829/1‑149 (MQ=255)
                                      aaTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTaaa                                                                                                              <  1:429880/149‑1 (MQ=255)
                                       aTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTaaaa                                                                                                             <  1:1508729/149‑1 (MQ=255)
                                           aaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATATAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTaa                                                                                                         >  2:168723/1‑149 (MQ=255)
                                           aaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTaa                                                                                                         <  1:827247/149‑1 (MQ=255)
                                             aaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAata                                                                                                       <  2:1525770/149‑1 (MQ=255)
                                               tGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAATCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATAtt                                                                                                     <  2:1017085/149‑1 (MQ=255)
                                               tGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATAtt                                                                                                     <  1:1256792/149‑1 (MQ=255)
                                                 aTCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTTAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAGTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGTCAGGCTGTATTAAAATTAATATTCa                                                                                                   >  2:977795/1‑149 (MQ=255)
                                                   cAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCaaa                                                                                                 <  2:1094675/149‑1 (MQ=255)
                                                    aaTCTTAATTTATTTAATTATGAGCTTTTTACTAAGTAATATAAAATATTTAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                                                                                <  1:271248/149‑1 (MQ=255)
                                                    aaTCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAg                                                                                                <  2:1314238/149‑1 (MQ=255)
                                                          aaTTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATaaa                                                                                          >  2:1474469/1‑149 (MQ=255)
                                                           atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTTAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                                                                         >  2:49311/1‑149 (MQ=255)
                                                           atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                                                                         >  2:219243/1‑149 (MQ=255)
                                                           atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                                                                         >  2:1544126/1‑149 (MQ=255)
                                                              tatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGa                                                                                      <  2:1581379/149‑1 (MQ=255)
                                                              tatttaATGATGAGCTTTTTAATCAGTAATATAAAATATTTAATTGTTATTTTTTTGAGTTTTTTAAGATAAAAAGCCGTATTTATAATTACGGGT‑‑‑‑TTAATTAATAAAATGCATTCTGTATTAAAATTAATACTCAAATCATAAACCGa                                                                                      <  2:107984/149‑1 (MQ=255)
                                                                tttaATTATCAGCTTTATACTCAGTAATATAAAATATTGACTAGTTATTTTTGTGTGTTGTTTAAGATAAAAAGGCGTATTTATTATTACGGCT‑‑‑‑TAAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                    <  2:105502/149‑1 (MQ=255)
                                                                tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                    <  2:1307938/149‑1 (MQ=255)
                                                                tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                    <  1:1082515/149‑1 (MQ=255)
                                                                tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                    <  2:533313/149‑1 (MQ=255)
                                                                tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                    <  1:1278407/149‑1 (MQ=255)
                                                                tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                    <  2:390106/149‑1 (MQ=255)
                                                                tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCAGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                    <  1:343389/149‑1 (MQ=255)
                                                                  taATGATGAGCTTTTTACTCAGTAATATAAAATATTTAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGc                                                                                  <  1:494752/149‑1 (MQ=255)
                                                                   aaTGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                                                                                 <  1:121215/149‑1 (MQ=255)
                                                                     tgatgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCaa                                                                               <  2:110386/149‑1 (MQ=255)
                                                                       atgaGATTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGAAAAAAATCAGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATa                                                                             <  2:57900/149‑1 (MQ=255)
                                                                        tgaGCTTTTTACTCAGTAATATAAAATATTGATTTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGAT‑‑‑‑TTATTTATTAAAAGGCAGGATTTATTAAAATTAATATTCAAAGCTTAAACCGATAGCCAATAc                                                                            >  1:1022942/1‑149 (MQ=255)
                                                                            cTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTATTTAATAAAAGGCAGGCTGTATTTAAATTAATATTCAAAGAATAAACCGATAGCCAATACCTGt                                                                        >  2:283674/1‑149 (MQ=255)
                                                                             tttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                                                                       <  2:775174/149‑1 (MQ=255)
                                                                             tttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                                                                       <  2:1151381/149‑1 (MQ=255)
                                                                              ttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                                                                      <  2:1021221/149‑1 (MQ=255)
                                                                              ttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                                                                      <  2:676580/149‑1 (MQ=255)
                                                                                ttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCa                                                                    <  2:378365/149‑1 (MQ=255)
                                                                                 tACTCAGTAATATAAAATATTTAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                                                                   <  1:49311/149‑1 (MQ=255)
                                                                                 tACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                                                                   <  2:365671/149‑1 (MQ=255)
                                                                                  aCTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGt                                                                  >  2:1454928/1‑149 (MQ=255)
                                                                                   cTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTa                                                                 >  2:558052/1‑149 (MQ=255)
                                                                                    tCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                                                                <  1:212927/149‑1 (MQ=255)
                                                                                      aGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATg                                                              >  1:1334456/1‑149 (MQ=255)
                                                                                              aaaaTATTGAATTGTTATTTTTGTGTGTTTTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACAGATAGCCAATAACGGTTTCAGTAATTAAATTGct                                                      >  2:1546606/1‑148 (MQ=255)
                                                                                               aaaTATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATg                                                         >  2:318437/1‑149 (MQ=38)
                                                                                                 atatTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgg                                                   <  2:452045/149‑1 (MQ=255)
                                                                                                   atTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGAGTgggc                                                 >  1:1083835/1‑149 (MQ=255)
                                                                                                    ttGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgggcg                                                <  1:1360298/149‑1 (MQ=255)
                                                                                                       aaTTGTTATTTTTGTGTGTTTTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgggcgggc                                             <  2:633434/149‑1 (MQ=255)
                                                                                                           gTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTTCGGCT‑‑‑‑TTAATTAATAAAAGGCAGGATGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCATTACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGa                                         >  2:727807/1‑149 (MQ=255)
                                                                                                             tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATc                                       <  1:1559040/149‑1 (MQ=255)
                                                                                                                            tttAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAATTTAATATTCAAAGCATAATCCGAGAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCTGTTTTTg                        >  2:69205/1‑149 (MQ=255)
                                                                                                                            tttAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCATGTTCCAGTTTTTg                        >  1:403362/1‑149 (MQ=255)
                                                                                                                            tttAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAAATAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTg                        >  1:862586/1‑149 (MQ=255)
                                                                                                                               aaGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCg                     >  1:542688/1‑149 (MQ=255)
                                                                                                                                ttaattaaaGCCGTAATAATAAATACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGACAATACCGGTTTCAGTAATGAAATGGCGTGGGCTGGCGGGATCCTGTTCCAGTTTTTGTCga                    >  1:1038589/7‑148 (MQ=25)
                                                                                                                                aGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCAAAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGc                    <  1:463196/149‑1 (MQ=255)
                                                                                                                                 gATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCa                   >  2:957358/1‑149 (MQ=255)
                                                                                                                                 gATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCa                   <  1:865597/149‑1 (MQ=255)
                                                                                                                                 gATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCa                   >  2:187842/1‑149 (MQ=255)
                                                                                                                                    aaaaaGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGAt                <  2:386478/149‑1 (MQ=255)
                                                                                                                                     aaaaGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATg               >  1:1551486/1‑149 (MQ=255)
                                                                                                                                        aGCCGTATTTATTATTACGGAT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAACGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTcc            <  1:428833/149‑1 (MQ=255)
                                                                                                                                          ccGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCa          >  1:755489/1‑149 (MQ=255)
                                                                                                                                           cGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCACAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCat         <  1:1584737/149‑1 (MQ=255)
                                                                                                                                             tatttattATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATTAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCatat       <  1:1466264/149‑1 (MQ=255)
                                                                                                                                                tttttattaCGGCT‑‑‑‑TTAATTTATAAACGGCAGTCTTTTTTATAATTATTATTCACAGCATAAACGGATAGCCAATACCGGTTTCTGTAATGAAATGGCGTTGGCGGTCTGGATCCTGTTCCAGATTTTGTCGCAGATGTCCCATATaaa    <  1:255840/146‑1 (MQ=255)
                                                                                                                                                  attattaCGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCATATAAATa  <  1:558052/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCATATAAATA  >  NC_000913/720861‑721155

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