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.4% 406.4 / ‑3.2 126intergenic (+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 (57/66);  total (59/67)
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.79e-01

Only 100 of 126 total aligned reads displayed.

AAGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTA  >  NC_000913/720849‑721042
                                                                                                 |                                                                                                
aaGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAg                                                                                                <  2:2242260/100‑1 (MQ=255)
  ggTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTa                                                                                              >  2:1901923/1‑100 (MQ=255)
  ggTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTa                                                                                              >  2:2221958/1‑100 (MQ=255)
   gTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTaa                                                                                             <  1:1056065/100‑1 (MQ=255)
       aGATTTTCTATCCAGGAAATAAAGATGAACTATAAAACTGATAAATAAAAACTGATCAAACTTAATTTATTTAATTATGAGCATTTTATTTCGTACtata                                                                                         >  2:2272877/1‑100 (MQ=255)
       aGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAAtata                                                                                         >  1:519512/1‑100 (MQ=255)
       aGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAAtata                                                                                         <  1:380073/100‑1 (MQ=255)
       aGATATTCTATCCAGGAAATACGGATGAACTGTGAAGCTGATCAATAAAAAATGAACAATCTTAATTTATTTAATGATGAGCTTTTTACTTAATAAtata                                                                                         >  2:1284712/1‑100 (MQ=255)
          tatTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATaaaa                                                                                      >  1:1895877/1‑100 (MQ=255)
          tatTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATaaaa                                                                                      <  1:1740865/100‑1 (MQ=255)
          tatTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATaaaa                                                                                      <  1:1740632/100‑1 (MQ=255)
           atTCTATCCAGGAAATACGGTTGAACTGTGAAAAAAAAAAAAAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAat                                                                                     <  2:37026/100‑1 (MQ=255)
           atTCTATCCAGGAAATACGATTTAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGTGCTTTTTACTTAGTAATATAAAat                                                                                     <  2:1362229/100‑1 (MQ=255)
            ttCTATCCAGGAAATACGGTTGAACTGTGAATGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAat                                                                                     <  2:962970/99‑1 (MQ=255)
             tCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAatat                                                                                   >  2:1977121/1‑100 (MQ=255)
             tCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAatat                                                                                   >  2:1091695/1‑100 (MQ=255)
               tATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTg                                                                                 <  1:1678257/100‑1 (MQ=255)
               tATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTg                                                                                 >  2:463134/1‑100 (MQ=255)
                  cgaGGAAATACGGTTAAACTGTGAAAATGATCAATAAAAAATGATCAATCTTAATGTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAt                                                                              <  2:1675759/98‑1 (MQ=255)
                  ccAGGAAATACGGTTTAACTGTGAAGGTGATCAAGACAAAATGCCCAATCTTAATTTATTTAATGATGAGCTCTTTCCTTAGTACTATAAAATATCGAAt                                                                              >  2:1738012/1‑100 (MQ=255)
                  ccAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAt                                                                              >  1:1068660/1‑100 (MQ=255)
                  ccAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTGAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAt                                                                              >  1:1992503/1‑100 (MQ=255)
                   cAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAtt                                                                             >  2:1595885/1‑100 (MQ=255)
                   cAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAtt                                                                             >  2:298048/1‑100 (MQ=255)
                   cAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAtt                                                                             >  2:55171/1‑100 (MQ=255)
                   cAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAtt                                                                             >  2:1483853/1‑100 (MQ=255)
                    aGGAAATACGTTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTGCTTAGTAATATAAAATATTGAATTg                                                                            <  2:2213375/100‑1 (MQ=255)
                    aGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGa                                                                                <  2:1444793/96‑1 (MQ=255)
                    aGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTg                                                                            <  2:1485574/100‑1 (MQ=255)
                    aGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTg                                                                            >  1:773883/1‑100 (MQ=255)
                     ggAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAt                                                                              <  1:1459743/97‑1 (MQ=255)
                        aaTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTCACTTAGTAATATAAAATATTGAATTGTTAt                                                                        >  1:1025572/1‑100 (MQ=255)
                             ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                   >  2:1864169/1‑100 (MQ=255)
                             ggTTGAACTGTGAAGGTAATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTTCTTAGTAATATAAAATGTTGAATTGTTATTTTtg                                                                   <  2:330363/100‑1 (MQ=255)
                             ggTTAAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                   <  2:330233/100‑1 (MQ=255)
                               tttaaCTGTGAAAATGAAAAAAAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAATATATTGAATTGTTATTTTtgtg                                                                 <  2:1087060/97‑1 (MQ=255)
                               ttGAACTGTGAAAGTGAAAAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtg                                                                 <  2:2071129/100‑1 (MQ=255)
                                tGAACTTTGTAAGTTATCAATAAATAATTTTCAATTTTAATTTATTTAATGATGAGTTTTTTACTTTGTAATATAAAATATTGAAt                                                                              <  2:1013611/86‑1 (MQ=255)
                                 gAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtgt                                                                >  1:943807/1‑99 (MQ=255)
                                  aaCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTgtgtgt                                                              <  2:773874/100‑1 (MQ=255)
                                    ntgtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGtgttg                                                            >  2:923553/2‑100 (MQ=255)
                                     tgttAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGtgtt                                                             <  2:474265/98‑1 (MQ=255)
                                     tgtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGtgttgt                                                           <  2:805506/100‑1 (MQ=255)
                                      gtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTgttgtt                                                          >  2:1914765/1‑100 (MQ=255)
                                       tgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGttt                                                         <  1:1901945/100‑1 (MQ=255)
                                        gAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTa                                                        <  1:2026713/100‑1 (MQ=255)
                                        gAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTa                                                        >  1:234332/1‑100 (MQ=255)
                                        gAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTa                                                        >  1:1700860/1‑100 (MQ=255)
                                            gTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGa                                                     <  2:2248246/99‑1 (MQ=255)
                                             tGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATa                                                   >  2:2156144/1‑100 (MQ=255)
                                              gATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaa                                                  >  1:105602/1‑100 (MQ=255)
                                                 cAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaaaa                                               >  1:798693/1‑100 (MQ=255)
                                                 cAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTAATTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTtaaaataaaa                                                >  2:195563/1‑99 (MQ=255)
                                                 aaaTAAAAAACGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaaaa                                               <  2:1568124/99‑1 (MQ=255)
                                                  aaTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTATTAATATAAAATATTGAATTGTTATTTTTCTTTGTTGTTTAAGATaaaa                                                >  2:743563/1‑98 (MQ=255)
                                                  aaTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATATAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAg                                              >  2:1688248/1‑100 (MQ=255)
                                                  aaTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATACAATATTGAATAGTTATTTTTGTGTGTTGTTTAAGATaaaaa                                               >  2:187274/1‑99 (MQ=255)
                                                   aTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGc                                             <  1:1790727/100‑1 (MQ=255)
                                                     aaaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCg                                           <  2:1892723/100‑1 (MQ=255)
                                                     aaaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCg                                           >  1:1125557/1‑100 (MQ=255)
                                                         aaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatt                                       >  1:545168/1‑100 (MQ=255)
                                                           tGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTTTGTGTTTTTTAAGATAAAAAGCAGtattta                                     >  2:2389400/1‑100 (MQ=255)
                                                            gATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttat                                    <  2:138123/100‑1 (MQ=255)
                                                             aTCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                   <  2:1244450/100‑1 (MQ=255)
                                                              tCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAGAGCCGTATttatta                                  <  2:2075826/100‑1 (MQ=255)
                                                               cAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                 <  1:1645180/100‑1 (MQ=255)
                                                               cAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                 <  1:2014549/100‑1 (MQ=255)
                                                               cAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                 <  2:775650/100‑1 (MQ=255)
                                                               cAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTATGTTGTTTAAGATAAAAAGCCGTATttattat                                 >  2:1879997/1‑100 (MQ=255)
                                                                 aTCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattatta                               <  1:2215075/100‑1 (MQ=255)
                                                                 aTCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattatta                               <  1:49096/100‑1 (MQ=255)
                                                                  tCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTAc                              <  2:2382654/100‑1 (MQ=255)
                                                                  gcTTAATTTATTTCATGATGAGCTTTTTACTAAGTAATATAAATTAATGACTTGTTATTTTTTATTGTTGTTTAATATAAAAATCCGTATttattatttc                              >  2:1691902/2‑98 (MQ=255)
                                                                   cTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACg                             <  2:1912252/100‑1 (MQ=255)
                                                                   cTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACg                             <  1:561372/100‑1 (MQ=255)
                                                                     tAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACgg                            >  2:497685/1‑99 (MQ=255)
                                                                     tAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTGCTTTAAGATAAAAAGCCGTATTTATTATTACg                             >  2:981898/1‑98 (MQ=255)
                                                                      aaTTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt                          <  2:24446/100‑1 (MQ=255)
                                                                         ttatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTtta                       <  2:486893/100‑1 (MQ=255)
                                                                         ttatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttt                       <  1:887744/100‑2 (MQ=38)
                                                                          tatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTtta                       <  2:1325488/99‑1 (MQ=255)
                                                                           atttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAATCCGTATTTATTATTACGGCTttaat                     >  2:1375213/1‑100 (MQ=255)
                                                                            tttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTTTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                    <  2:624532/100‑1 (MQ=255)
                                                                            tttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                    <  1:38964/100‑1 (MQ=255)
                                                                             ttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatc                   <  2:594529/100‑6 (MQ=25)
                                                                              taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                  <  1:2017706/100‑1 (MQ=255)
                                                                              taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                  >  1:1569338/1‑100 (MQ=255)
                                                                              taATGATGAGCTTTTTACTTAGTAATATAAAAAATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                  <  2:113270/100‑7 (MQ=25)
                                                                               aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                 <  1:2280474/100‑1 (MQ=255)
                                                                               aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                 <  1:1911701/100‑8 (MQ=25)
                                                                                 tgatgcGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                 <  1:1914787/98‑1 (MQ=255)
                                                                                 tgatgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTTTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaataa               >  2:2211364/1‑100 (MQ=255)
                                                                                    tgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGGTGTTTAAGATAAAAAGCCGTATTTATTATGACGGCTTTAATTAATaaaa             >  2:2372292/1‑99 (MQ=255)
                                                                                      aGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggc          >  1:1193050/1‑100 (MQ=255)
                                                                                          ttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc      >  1:1468815/1‑100 (MQ=255)
                                                                                             tACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGt   >  2:181757/1‑100 (MQ=255)
                                                                                             tACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGt   >  1:2195272/1‑100 (MQ=255)
                                                                                              aCTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTa  >  2:649087/1‑100 (MQ=255)
                                                                                              aCTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTa  >  1:1507072/1‑100 (MQ=255)
                                                                                              aCTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTa  >  1:379965/1‑100 (MQ=255)
                                                                                                 |                                                                                                
AAGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTA  >  NC_000913/720849‑721042

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