Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 721,011 T→A 5.2% intergenic (+171/+45) 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,0110TA 5.2% 402.2 / 7.5 131intergenic (+171/+45)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base T (51/73);  new base A (3/4);  total (54/77)
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.76e-01

Only 100 of 138 total aligned reads displayed.

ATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCATATAAATACGCA  >  NC_000913/720865‑721159
                                                                                                                                                  |                                                                                                                                                        
aTCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATTATGAGATTTTTACTCAGTAAAATAAAATATTGAATTGTTATTTTTTTGTTTTTTTTAAGATAAAAAGCCGAATttattatta                                                                                                                                                        >  1:633434/1‑149 (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:110386/1‑149 (MQ=255)
     tgAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt                                                                                                                                                   <  2:1486861/148‑1 (MQ=255)
      gAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt‑‑‑‑t                                                                                                                                              <  2:1406882/149‑1 (MQ=255)
       aaaTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTAATAAATACGGCt‑‑‑‑tt                                                                                                                                             <  1:1316692/149‑1 (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)
            cGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑ttaatta                                                                                                                                        >  2:1559040/1‑149 (MQ=255)
            cGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCt‑‑‑‑ttttatc                                                                                                                                        <  2:1303642/149‑6 (MQ=255)
             ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑ttaattaa                                                                                                                                       <  1:253544/149‑1 (MQ=255)
               ttGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑Ttaattaat                                                                                                                                      >  1:1284893/1‑148 (MQ=255)
               ttGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑Ttaattaat                                                                                                                                      <  2:1284893/148‑1 (MQ=255)
               ttGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATtaata                                                                                                                                     >  1:1413704/1‑149 (MQ=255)
                  aaCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATaaaa                                                                                                                                  >  2:1584737/1‑149 (MQ=255)
                  aaCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATaaaa                                                                                                                                  <  2:1209797/149‑1 (MQ=255)
                  aaCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATaaaa                                                                                                                                  <  1:303396/149‑1 (MQ=255)
                  aaCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATaaaa                                                                                                                                  >  2:343389/1‑149 (MQ=255)
                    cTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAagg                                                                                                                                >  2:864762/1‑149 (MQ=255)
                    cTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAagg                                                                                                                                >  2:463196/1‑149 (MQ=255)
                    cTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGATTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCTTATTTATTATTACGGAT‑‑‑‑TTAATTAATAAAaga                                                                                                                                >  2:1059660/1‑148 (MQ=255)
                     tgtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAaggc                                                                                                                               >  2:669737/1‑149 (MQ=255)
                      gtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAaggca                                                                                                                              <  2:1421904/149‑1 (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)
                              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)
                                       aaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATATAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTaa                                                                                                             >  2:168723/1‑149 (MQ=255)
                                       aaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTaa                                                                                                             <  1:1241703/149‑1 (MQ=255)
                                           tGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATAtt                                                                                                         <  1:1256792/149‑1 (MQ=255)
                                            taTCAATCTTAATTTTTTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTc                                                                                                        <  1:310043/148‑1 (MQ=255)
                                             aTCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTTAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAGTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGTCAGGCTGTATTAAAATTAATATTCa                                                                                                       >  2:977795/1‑149 (MQ=255)
                                                    ttAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATa                                                                                                >  1:533313/1‑149 (MQ=255)
                                                      aaTTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATaaa                                                                                              >  2:1474469/1‑149 (MQ=255)
                                                       atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                                                                             >  2:1544126/1‑149 (MQ=255)
                                                       atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                                                                             >  2:219243/1‑149 (MQ=255)
                                                       atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                                                                             <  1:1546808/149‑1 (MQ=255)
                                                          tatttaATGATGAGCTTTTTAATCAGTAATATAAAATATTTAATTGTTATTTTTTTGAGTTTTTTAAGATAAAAAGCCGTATTTATAATTACGGGT‑‑‑‑TTAATTAATAAAATGCATTCTGTATTAAAATTAATACTCAAATCATAAACCGa                                                                                          <  2:107984/149‑1 (MQ=255)
                                                            tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                        <  1:1082515/149‑1 (MQ=255)
                                                            tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                        <  1:380474/149‑1 (MQ=255)
                                                            tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                        <  2:390106/149‑1 (MQ=255)
                                                            tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                        <  2:533313/149‑1 (MQ=255)
                                                            tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                        <  2:1307938/149‑1 (MQ=255)
                                                            tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCAGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                        <  1:343389/149‑1 (MQ=255)
                                                             ttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGACGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAg                                                                                       >  1:759917/1‑149 (MQ=255)
                                                              taATGATGAGCTTTTTACTCAGTAATATAAAATATTTAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGc                                                                                      <  1:494752/149‑1 (MQ=255)
                                                               aaTGATGAGCTTTTTACTCAGTAATATAAAATATTTAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGcc                                                                                     <  1:1564699/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)
                                                                    tgaGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATAc                                                                                <  1:443084/149‑1 (MQ=255)
                                                                        cTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                                                                            <  2:195232/149‑1 (MQ=255)
                                                                         tttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                                                                           <  2:1151381/149‑1 (MQ=255)
                                                                         tttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGtt                                                                           <  2:775174/149‑1 (MQ=255)
                                                                          ttttACTCAGTAATATAAAATATTTAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                                                                          <  2:1110282/149‑1 (MQ=255)
                                                                          ttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                                                                          <  2:676580/149‑1 (MQ=255)
                                                                          ttttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGttt                                                                          <  2:1021221/149‑1 (MQ=255)
                                                                            ttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCa                                                                        <  2:378365/149‑1 (MQ=255)
                                                                             tACTCAGTAATATAAAATATTTAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                                                                       <  1:49311/149‑1 (MQ=255)
                                                                                tCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                                                                    <  1:212927/149‑1 (MQ=255)
                                                                                tCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                                                                    >  2:407061/1‑149 (MQ=255)
                                                                                  aGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATg                                                                  >  1:1334456/1‑149 (MQ=255)
                                                                                  aGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATg                                                                  <  2:1082017/149‑1 (MQ=255)
                                                                                         taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGc                                                           >  2:575713/1‑149 (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)
                                                                                                   aaTTGTTATTTTTGTGTGTTTTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTgggcgggc                                                 <  2:633434/149‑1 (MQ=255)
                                                                                                       gTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTTCGGCT‑‑‑‑TTAATTAATAAAAGGCAGGATGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCATTACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGa                                             >  2:727807/1‑149 (MQ=255)
                                                                                                        ttATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGAt                                            >  2:1320753/1‑149 (MQ=255)
                                                                                                         tATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATc                                           <  1:1559040/149‑1 (MQ=255)
                                                                                                               tgtgtgTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTc                                     <  2:367291/149‑1 (MQ=255)
                                                                                                               tgtgtgTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTc                                     <  1:926547/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)
                                                                                                                            ttaattaaaGCCGTAATAATAAATACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGACAATACCGGTTTCAGTAATGAAATGGCGTGGGCTGGCGGGATCCTGTTCCAGTTTTTGTCga                        >  1:1038589/7‑148 (MQ=25)
                                                                                                                            aGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGc                        >  2:166506/1‑149 (MQ=255)
                                                                                                                             gATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCa                       >  2:957358/1‑149 (MQ=255)
                                                                                                                             gATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCa                       <  1:865597/149‑1 (MQ=255)
                                                                                                                             gATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGAAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCTTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCa                       >  2:493384/1‑149 (MQ=255)
                                                                                                                                   aaGCCGTAATAATAAATACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTc                 >  2:1210468/1‑149 (MQ=255)
                                                                                                                                    aGCCGTATTTATTATTACGGAT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAACGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTcc                <  1:428833/149‑1 (MQ=255)
                                                                                                                                      ccGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCa              >  1:755489/1‑149 (MQ=255)
                                                                                                                                      ccGTAATAATAAATACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCa              >  2:1545839/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      >  2:1028789/1‑149 (MQ=255)
                                                                                                                                              attattaCGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCATATAAATa      <  1:558052/149‑1 (MQ=255)
                                                                                                                                                 attaCGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAAAATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCATATAAATACGc   <  1:363251/149‑1 (MQ=255)
                                                                                                                                                  ttaCGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCATATAAATACGCa  >  2:1422128/1‑149 (MQ=255)
                                                                                                                                                  ttaCGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCATGTTCCAGTTTTTGTCGCAGATGTCACATATAAATACGCa  >  1:1422812/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                        
ATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCATATAAATACGCA  >  NC_000913/720865‑721159

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: