Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 721,006 T→A 5.1% intergenic (+166/+50) 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,0060TA 5.1% 379.6 / 4.5 129intergenic (+166/+50)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base T (51/71);  new base A (3/4);  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 = 9.97e-01

Only 100 of 136 total aligned reads displayed.

TTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCATATA  >  NC_000913/720861‑721151
                                                                                                                                                 |                                                                                                                                                     
ttCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                                                                                                                                    <  2:510242/149‑1 (MQ=255)
ttCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                                                                                                                                    <  1:1234782/149‑1 (MQ=255)
 tCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttatta                                                                                                                                                   <  2:1540364/149‑1 (MQ=255)
    aTCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATTATGAGATTTTTACTCAGTAAAATAAAATATTGAATTGTTATTTTTTTGTTTTTTTTAAGATAAAAAGCCGAATttattatta                                                                                                                                                >  1:633434/1‑149 (MQ=255)
     tCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTAc                                                                                                                                               <  1:1032605/149‑1 (MQ=255)
     tCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTAc                                                                                                                                               >  2:428833/1‑149 (MQ=255)
        aGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                                                                            >  1:110386/1‑149 (MQ=255)
        aGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                                                                            >  1:390106/1‑149 (MQ=255)
        aGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                                                                                                                                            >  2:1546808/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)
               aCGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCt‑‑‑‑ttttat                                                                                                                                 <  2:1136943/149‑5 (MQ=255)
                cGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑ttaatta                                                                                                                                <  1:3274/149‑1 (MQ=255)
                cGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑ttaatta                                                                                                                                >  2:1559040/1‑149 (MQ=255)
                 ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑ttaattaa                                                                                                                               <  1:253544/149‑1 (MQ=255)
                   ttGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑Ttaattaat                                                                                                                              <  2:1284893/148‑1 (MQ=255)
                   ttGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑Ttaattaat                                                                                                                              >  1:1284893/1‑148 (MQ=255)
                   ttGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATtaata                                                                                                                             >  1:1413704/1‑149 (MQ=255)
                      aaCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATaaaa                                                                                                                          >  2:343389/1‑149 (MQ=255)
                      aaCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATaaaa                                                                                                                          >  2:1584737/1‑149 (MQ=255)
                      aaCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATaaaa                                                                                                                          >  2:1360298/1‑149 (MQ=255)
                        cTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAagg                                                                                                                        >  2:463196/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:1421904/149‑1 (MQ=255)
                           tgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAaggcag                                                                                                                     >  1:1525770/1‑149 (MQ=255)
                                gTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGATCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGt                                                                                                                <  1:415315/149‑1 (MQ=255)
                                 tGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTa                                                                                                               <  1:647078/149‑1 (MQ=255)
                                 tGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTa                                                                                                               <  2:1279436/149‑1 (MQ=255)
                                  gATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTAt                                                                                                              <  2:1478237/149‑1 (MQ=255)
                                  gATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTAt                                                                                                              >  2:681318/1‑149 (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)
                                           aaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTaa                                                                                                     <  1:1241703/149‑1 (MQ=255)
                                                taTCAATCTTAATTTTTTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTc                                                                                                <  1:310043/148‑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)
                                                           atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                                                                     >  2:1544126/1‑149 (MQ=255)
                                                           atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                                                                     >  2:219243/1‑149 (MQ=255)
                                                           atttatttAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAAc                                                                                     <  1:1546808/149‑1 (MQ=255)
                                                              tatttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGa                                                                                  <  2:1581379/149‑1 (MQ=255)
                                                                tttaATTATCAGCTTTATACTCAGTAATATAAAATATTGACTAGTTATTTTTGTGTGTTGTTTAAGATAAAAAGGCGTATTTATTATTACGGCT‑‑‑‑TAAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                <  2:105502/149‑1 (MQ=255)
                                                                tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                <  2:533313/149‑1 (MQ=255)
                                                                tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                <  2:1307938/149‑1 (MQ=255)
                                                                tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                <  1:1278407/149‑1 (MQ=255)
                                                                tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                <  2:390106/149‑1 (MQ=255)
                                                                tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                <  1:380474/149‑1 (MQ=255)
                                                                tttaATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATa                                                                                <  1:1082515/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)
                                                                   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‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGt                                                                    <  2:195232/149‑1 (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)
                                                                                ttACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCa                                                                <  2:378365/149‑1 (MQ=255)
                                                                                 tACTCAGTAATATAAAATATTTAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                                                               <  1:49311/149‑1 (MQ=255)
                                                                                 tACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAg                                                               <  2:365671/149‑1 (MQ=255)
                                                                                    tCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                                                            <  1:212927/149‑1 (MQ=255)
                                                                                    tCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTaa                                                            >  2:407061/1‑149 (MQ=255)
                                                                                      aGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATg                                                          <  2:1082017/149‑1 (MQ=255)
                                                                                             taAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGc                                                   >  2:575713/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)
                                                                                                                   tgtgtgTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTc                             <  1:926547/149‑1 (MQ=255)
                                                                                                                   tgtgtgTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTc                             <  2:367291/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‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGc                >  2:166506/1‑149 (MQ=255)
                                                                                                                                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)
                                                                                                                                 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)
                                                                                                                                           cGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCACAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCat     <  1:1584737/149‑1 (MQ=255)
                                                                                                                                              atttattATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTTTATTAAAATTAATATTCAACGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCAtata  <  2:599789/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCT‑‑‑‑TTAATTAATAAAAGGCAGGCTGTATTAAAATTAATATTCAAAGCATAAACCGATAGCCAATACCGGTTTCAGTAATGAAATGGCGTGGGCGGGCGGGATCCTGTTCCAGTTTTTGTCGCAGATGTCCCATATA  >  NC_000913/720861‑721151

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: