Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 720,946 C→T 100% intergenic (+106/+110) 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_000913720,9460CT98.7% 484.5 / ‑3.1 150intergenic (+106/+110)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (79/68);  total (81/69)
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.54e-01

Only 100 of 150 total aligned reads displayed.

AGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTA  >  NC_000913/720850‑721042
                                                                                                |                                                                                                
aGGTGAAGATATTCTATCCAGGAAATACGGTTGTACTGTGAAGGTGATCAATAAAAAATGATCCATCTTCATTTATTTAATGATGAGCTTTTTACTTAgc                                                                                               <  1:1340438/100‑2 (MQ=255)
   tGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTaa                                                                                             >  1:2881229/1‑99 (MQ=255)
   tGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAat                                                                                            >  2:599132/1‑100 (MQ=255)
    gAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGCTGAGCTTTTTACTTAGTAata                                                                                           >  2:3049621/1‑100 (MQ=255)
      aGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAAtata                                                                                         >  2:128673/1‑100 (MQ=255)
      aGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAAtata                                                                                         >  1:2747988/1‑100 (MQ=255)
      aGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAAtata                                                                                         >  2:1794571/1‑100 (MQ=255)
      aGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAAtata                                                                                         >  1:2947392/1‑100 (MQ=255)
      aGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAAtata                                                                                         >  2:728753/1‑100 (MQ=255)
          atTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAat                                                                                     >  2:3389275/1‑100 (MQ=255)
          atTCTATCCAGAAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAat                                                                                     <  1:2210362/100‑1 (MQ=255)
           ttCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAata                                                                                    >  1:4100/1‑100 (MQ=255)
           ttCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAata                                                                                    >  2:2908146/1‑100 (MQ=255)
              tATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTg                                                                                 >  2:627992/1‑100 (MQ=255)
              tATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTg                                                                                 <  1:2260350/100‑1 (MQ=255)
              tATCCAGGAAATACGGTTGAACTGTGAAGGTGATAAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTg                                                                                 <  2:846722/100‑1 (MQ=255)
               aTCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGa                                                                                >  1:3094698/1‑100 (MQ=255)
                  cAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAtt                                                                             >  1:397184/1‑100 (MQ=255)
                         tACAGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAttt                                                                      <  1:596422/100‑1 (MQ=255)
                          aCGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAtttt                                                                     >  2:2102408/1‑100 (MQ=255)
                           cGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAttttt                                                                    <  2:1125313/100‑1 (MQ=255)
                           cGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAttttt                                                                    >  2:1886854/1‑100 (MQ=255)
                           cGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTAttttt                                                                    >  1:862071/1‑100 (MQ=255)
                            ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                   >  1:1639740/1‑100 (MQ=255)
                            ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                   >  1:2702915/1‑100 (MQ=255)
                             gTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgt                                                                  >  2:2093935/1‑100 (MQ=255)
                              ttGAACTGTGTAGGTGATCAATAAACTATGCTCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtg                                                                 <  1:280354/100‑1 (MQ=255)
                              ttGAACTGTGAAGGTGTTCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtg                                                                 <  1:1976375/100‑1 (MQ=255)
                                gAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTGATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtgtg                                                               >  2:1169755/1‑100 (MQ=255)
                                gAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtgtg                                                               >  1:2514480/1‑100 (MQ=255)
                                 aaCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTgtgtgt                                                              <  2:2263036/100‑1 (MQ=255)
                                  aCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATATAATGATGAGCTTTTTACTTAGTCATATAAATTATTGAATTGTTATTATTATGtgtt                                                             >  1:1852519/1‑100 (MQ=255)
                                   cTGTTAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGtgttg                                                            <  1:1499299/100‑1 (MQ=255)
                                    tgtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGtgttgt                                                           <  2:1258742/100‑1 (MQ=255)
                                     gtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTGAATGATGAGCTAAATACTTAGTAATATAAAATAATGAATTGTTAATTTtgat                                                                 >  1:3031576/1‑91 (MQ=255)
                                      tgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGttt                                                         >  2:1519383/1‑100 (MQ=255)
                                       gAAGATAATAAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTa                                                        <  1:285474/100‑1 (MQ=255)
                                         aGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAg                                                      <  1:440361/100‑1 (MQ=255)
                                         aGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAg                                                      <  2:926141/100‑1 (MQ=255)
                                          ggTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGa                                                     >  2:1805157/1‑100 (MQ=255)
                                            ttattaaTAAAAAATGATCAATCTTAATTTATGTAATGATTTGCTTTTTACTTAGTAATATAAAATATTGAATTTTTATTTTTGTGTGTTGTTTAAg                                                      <  1:989708/92‑1 (MQ=255)
                                             gATTAATAAAAAATTATCAATCTTAATTGATTTAATGATGAGCTTTTTATTTAGTAATATAAAATATTGAATAGTCATTTTTGTGTGTTGTTTTAGAt                                                    <  1:2507444/98‑1 (MQ=255)
                                             gATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAGATATTCAATTGTTTATTTTGTATGTTGTTTAAtata                                                   >  1:2283154/1‑99 (MQ=255)
                                             gATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaa                                                  >  1:258664/1‑100 (MQ=255)
                                             gATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaa                                                  >  1:512588/1‑100 (MQ=255)
                                              aTCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTTTTTAATAATATAAAATATTGAATTTTTATTTTTGTGTGTTATTTAAAATa                                                   >  1:218427/1‑98 (MQ=255)
                                              aTCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaa                                                 >  2:2741230/1‑100 (MQ=255)
                                               tCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaaa                                                <  2:316170/100‑1 (MQ=255)
                                               tCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaaa                                                <  2:79060/100‑1 (MQ=255)
                                               tCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaaa                                                <  1:1140695/100‑1 (MQ=255)
                                                cAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTATTTTTTTAATATATTTAATCGTTATTTATGTGTGTTTTTTGAGATaata                                                >  1:1331691/1‑97 (MQ=255)
                                                cAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaaaa                                               >  2:1588214/1‑100 (MQ=255)
                                                 aaTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAg                                              >  2:329567/1‑100 (MQ=255)
                                                  aTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGc                                             <  2:1685316/100‑1 (MQ=255)
                                                  aTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGc                                             <  2:2368995/100‑1 (MQ=255)
                                                   tAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGcc                                            <  2:26044/100‑1 (MQ=255)
                                                    aaaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTTTGTGTTGTTTAAGATataaattcg                                           >  1:561330/1‑96 (MQ=255)
                                                    aaaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCg                                           <  1:1442700/100‑1 (MQ=255)
                                                      aaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGta                                         <  1:1606261/100‑1 (MQ=255)
                                                      aaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGta                                         >  2:2764979/1‑100 (MQ=255)
                                                      aaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGta                                         <  1:2368401/100‑1 (MQ=255)
                                                       aaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtat                                        >  2:1925539/1‑100 (MQ=255)
                                                        aaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatt                                       <  1:2925669/100‑1 (MQ=255)
                                                            aTCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttat                                    <  2:864561/99‑1 (MQ=255)
                                                            aTCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTGAAGATAAAAAGCCGtatttatt                                   <  1:1961088/100‑1 (MQ=255)
                                                             tCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttatta                                  <  1:1184461/100‑1 (MQ=255)
                                                             tCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttatta                                  >  1:209284/1‑100 (MQ=255)
                                                             tCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttatta                                  <  2:1933626/100‑1 (MQ=255)
                                                               aaTCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTGTTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattatt                                <  2:395235/100‑1 (MQ=255)
                                                                aTCTTAATTTATTTAATGATGAGCTTGTTACTTACTAATATAAAATATTGAATTGTTATTTTTGTGTGCTGTTTAAGATAAAAAGCCGTATttatta                                  >  1:1790551/1‑97 (MQ=255)
                                                                 tCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTATTTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTAc                              >  1:461317/1‑100 (MQ=255)
                                                                   ttAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACgg                            <  1:3167754/100‑1 (MQ=255)
                                                                    tAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                           >  2:3102416/1‑100 (MQ=255)
                                                                     aaTTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTTTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt                          <  1:1851076/100‑1 (MQ=255)
                                                                      ntttatttaATGATGAGCTTTTTTTTTTTTTTTATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                         >  1:1214440/2‑100 (MQ=255)
                                                                      atttatttAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                         >  2:2608229/1‑100 (MQ=255)
                                                                      atttatttAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtt                         <  1:1098897/100‑1 (MQ=255)
                                                                      atttatttAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtt                         <  1:1098805/100‑1 (MQ=255)
                                                                       tttatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCttt                        >  2:1758135/1‑100 (MQ=255)
                                                                        ttatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTtta                       <  2:3286217/100‑1 (MQ=255)
                                                                        ttatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTtta                       <  1:713616/100‑1 (MQ=255)
                                                                           tttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                    <  1:858293/100‑1 (MQ=255)
                                                                           tttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                    <  1:995976/100‑1 (MQ=255)
                                                                           tttaATGATGAGCTTTTTACTTAGTAATATAAAAAATTGAATTGTTATTTTTGTGTGTTGTAAAAAATAAAAAGCCGTATTTATTATTACGGCTttaatt                    <  1:1142527/100‑1 (MQ=255)
                                                                            ttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatc                   <  2:73621/100‑6 (MQ=25)
                                                                             taATGATGAGCTTTTTACTTAGTAATATAAAATATTGGATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                  >  1:1810529/1‑100 (MQ=255)
                                                                             taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                  <  1:1662272/100‑1 (MQ=255)
                                                                              aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCTGTATTTATTATTACGGCTTtaattaat                 >  1:76310/1‑100 (MQ=255)
                                                                              aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                 <  1:1574012/100‑1 (MQ=255)
                                                                               atgatgAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                >  1:3323538/1‑100 (MQ=255)
                                                                                 gatgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaa              >  2:1848503/1‑100 (MQ=255)
                                                                                   tgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAag            <  2:2965014/100‑1 (MQ=255)
                                                                                        tttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcagg       <  2:476287/100‑1 (MQ=255)
                                                                                         ttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc      <  1:24547/100‑1 (MQ=255)
                                                                                         ttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc      >  1:520526/1‑100 (MQ=255)
                                                                                         ttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc      <  1:3377980/100‑1 (MQ=255)
                                                                                         ttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc      >  2:2741750/1‑100 (MQ=255)
                                                                                            tACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGt   >  2:336767/1‑100 (MQ=255)
                                                                                            tACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGt   >  1:2653436/1‑100 (MQ=255)
                                                                                             aCTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTgca  >  1:2048116/1‑98 (MQ=255)
                                                                                                |                                                                                                
AGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTA  >  NC_000913/720850‑721042

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