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.3% 573.5 / ‑1.8 176intergenic (+106/+110)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (80/93);  total (81/95)
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.76e-01

Only 100 of 176 total aligned reads displayed.

AAGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTA  >  NC_000913/720849‑721042
                                                                                                 |                                                                                                
 aGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGt                                                                                               <  1:798546/100‑1 (MQ=255)
 aGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGt                                                                                               <  1:2952252/100‑1 (MQ=255)
  ggTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTa                                                                                              <  2:2458988/100‑1 (MQ=255)
   gTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTaa                                                                                             >  2:3669893/1‑100 (MQ=255)
    tGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAat                                                                                            >  2:1911704/1‑100 (MQ=255)
           atTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAat                                                                                     <  2:3416766/100‑1 (MQ=255)
           atTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAat                                                                                     >  1:1337640/1‑100 (MQ=255)
           atTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAat                                                                                     >  2:1604187/1‑100 (MQ=255)
            ttCTATCCGGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAata                                                                                    >  1:1260293/1‑100 (MQ=255)
               tataCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTg                                                                                 <  1:839343/100‑1 (MQ=255)
               tATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTg                                                                                 >  1:1951528/1‑100 (MQ=255)
                  ccAGGAAATACGGTTGTACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTgtct                                                                              >  1:3028572/1‑97 (MQ=255)
                  ccAGGAAATACGCTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGact                                                                              >  1:3579239/1‑98 (MQ=255)
                   cAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAtt                                                                             >  1:1086535/1‑100 (MQ=255)
                    aGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTg                                                                            <  1:2477110/100‑1 (MQ=255)
                    aGGAAATACAGTTAAAATATAAAAGTAATAAAAAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAat                                                                                     <  1:709086/91‑1 (MQ=255)
                       aaaTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTa                                                                         <  1:1497592/100‑1 (MQ=255)
                        aaTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTa                                                                         >  1:1405123/1‑99 (MQ=255)
                         aTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAtt                                                                       <  2:2664886/100‑1 (MQ=255)
                          tACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAttt                                                                      >  2:681086/1‑100 (MQ=255)
                           aCGGTTGAACTGTGAAGGTGATCAATAATAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAtttt                                                                     <  1:4142329/100‑1 (MQ=255)
                           aCGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAtttt                                                                     >  1:570190/1‑100 (MQ=255)
                            cGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAttttt                                                                    >  1:4220555/1‑100 (MQ=255)
                            cGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAttttt                                                                    >  1:498209/1‑100 (MQ=255)
                            cGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAttttt                                                                    >  1:865261/1‑100 (MQ=255)
                             ggTTGAACTGTGAAGGTGATAAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                   <  1:232152/100‑1 (MQ=255)
                              gTTGAACTGTGAAGGTGATCACTAAAAAATTATCAACCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                   <  1:2875015/99‑1 (MQ=255)
                               ttGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtg                                                                 <  2:632161/100‑1 (MQ=255)
                               ttGAACTGTGAAAAAAAAAAAAAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtg                                                                 <  1:1998858/100‑1 (MQ=255)
                                 gAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtgtg                                                               <  2:319155/100‑1 (MQ=255)
                                    cTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGtgttg                                                            >  2:3023864/1‑100 (MQ=255)
                                      gtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTgtgtgt                                                              >  1:4124608/1‑96 (MQ=255)
                                      gtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTgttgtt                                                          >  2:2077444/1‑100 (MQ=255)
                                        gAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATCTTCAATTGTTATTTTTGTGTGTTGTTTa                                                        >  2:1447251/1‑100 (MQ=255)
                                          aGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAg                                                      >  2:3581167/1‑100 (MQ=255)
                                          aGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAg                                                      <  2:1730813/100‑1 (MQ=255)
                                          aGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAATACTATAAATTATTGAATTGTTATTTTTGTACGTTGTTTAAg                                                      >  1:3003700/1‑100 (MQ=255)
                                           ggTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGa                                                     <  1:3366203/100‑1 (MQ=255)
                                           ggTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGa                                                     <  2:1046787/100‑1 (MQ=255)
                                            gTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGa                                                     <  1:2166569/99‑1 (MQ=255)
                                             tGATCAATAAAGAATGATCAATCTTAATTGATGTAAGGATGAGCTTGTTACTTAGTGATATAAAATATGGTGTGGTGATTATTGTGTTTTGTTTAAGATa                                                   >  1:2546437/1‑100 (MQ=255)
                                              gATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaa                                                  >  2:3645207/1‑100 (MQ=255)
                                               aTCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaa                                                 <  2:2175082/100‑1 (MQ=255)
                                                tCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTTTTTTTTTTTTTAAAATATTGAATTGTTATTTTTTTTTGTTGTTTAAGATaaa                                                 >  1:1635164/1‑99 (MQ=255)
                                                tCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaaa                                                <  1:1604183/100‑1 (MQ=255)
                                                 taaTAAAAAATGATCAATCTTAATTTATTTAATGATTAGATTTTTACTTAATAATATAAAATATTTAATTGTTATTTTTGTTTGTTGTTTAAAATaaaaa                                               >  1:3514531/2‑100 (MQ=255)
                                                 cAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaaaa                                               >  2:3837673/1‑100 (MQ=255)
                                                  aaTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAg                                              >  2:3208553/1‑100 (MQ=255)
                                                  aaTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTAATTAGTAATATAAAATATTGAATTGTTATTTTTGTGTTTTGTTTAAGATaaaa                                                >  1:3933043/1‑98 (MQ=255)
                                                   aTAAAAAATGATCAATCTTAATTTATTTAATGATGAGTTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGc                                             <  1:3532165/100‑1 (MQ=255)
                                                   aTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGc                                             <  1:840916/100‑1 (MQ=255)
                                                   aTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGc                                             >  2:4070852/1‑100 (MQ=255)
                                                      aaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGt                                          >  1:263452/1‑100 (MQ=255)
                                                       aaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTTTTTTTGTGTGTTGTTTAAGATAAAAAGCCGta                                         >  1:4007918/1‑100 (MQ=255)
                                                       aaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGta                                         >  2:232123/1‑100 (MQ=255)
                                                         aaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatt                                       <  2:1128347/100‑1 (MQ=255)
                                                         aaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatt                                       <  1:3649535/100‑1 (MQ=255)
                                                            gATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATATAAAGACGtatttat                                    >  1:3161125/1‑100 (MQ=255)
                                                             aTCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                   >  2:257422/1‑100 (MQ=255)
                                                                 aTCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattatta                               >  2:2173751/1‑100 (MQ=255)
                                                                   cTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGTTAAAAAGCCGTATTTATTATTACg                             <  1:230739/100‑1 (MQ=255)
                                                                    ttAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACgg                            <  2:1288670/100‑1 (MQ=255)
                                                                    ttAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGGGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACgg                            >  2:1474843/1‑100 (MQ=255)
                                                                     tAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                           <  2:1110239/100‑1 (MQ=255)
                                                                     tAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCATATTTATTATTACGGc                           <  2:331673/100‑1 (MQ=255)
                                                                      aaTTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt                          <  2:595317/100‑1 (MQ=255)
                                                                       atttatttAATGATGAGCTTTTTACTTAGTAATATAAAATTTTGAATTGTTATTTTTGTGTGTTGTTAAAGATAAAAAGCCGTATTTATTATTACGGCtt                         <  1:219684/100‑1 (MQ=255)
                                                                       atttatttAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                         >  1:1058535/1‑100 (MQ=255)
                                                                       atttatttAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGGTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                         <  1:2285782/100‑1 (MQ=255)
                                                                          tatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaa                      <  2:2108394/100‑1 (MQ=255)
                                                                          tatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaa                      >  2:405912/1‑100 (MQ=255)
                                                                           atttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaat                     >  1:2184838/1‑100 (MQ=255)
                                                                           atttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaat                     <  2:1758395/100‑1 (MQ=255)
                                                                            tttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaat                     <  2:869264/99‑1 (MQ=255)
                                                                              taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAATATAAAAAGCCGTAATAATAAATACGGCtttttatct                  <  1:438550/100‑7 (MQ=25)
                                                                              taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                  <  2:1989119/100‑1 (MQ=255)
                                                                              taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                  >  2:1997067/1‑100 (MQ=255)
                                                                              taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaattaa                  <  1:1093263/100‑1 (MQ=255)
                                                                              taATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                  <  1:1341998/100‑7 (MQ=37)
                                                                               aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                 <  1:3784806/100‑1 (MQ=255)
                                                                                atgatgAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTTTGTGTAATTTAAAATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                 <  1:181331/99‑1 (MQ=255)
                                                                                atgatgAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                <  2:4141336/100‑1 (MQ=255)
                                                                                atgatgAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                <  2:3599230/100‑1 (MQ=255)
                                                                                  gatgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTTTTTTTTTTTTAAGATAAAAAGCCGTATTTATTATTACGTCTTTAATTAATaaa              >  1:27051/1‑100 (MQ=255)
                                                                                  gatgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatcttaaa              <  2:1257538/100‑11 (MQ=21)
                                                                                   atgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaaa             >  2:399344/1‑100 (MQ=255)
                                                                                    tgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAag            <  2:3249696/100‑1 (MQ=255)
                                                                                     gaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaaa             >  1:2701352/1‑98 (MQ=255)
                                                                                      agttttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggc          <  1:3782348/97‑1 (MQ=255)
                                                                                      aGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggc          <  2:2923379/100‑1 (MQ=255)
                                                                                       gCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca         >  2:1410776/1‑100 (MQ=255)
                                                                                        cTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag        >  2:371702/1‑100 (MQ=255)
                                                                                         tttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcagg       >  1:206201/1‑100 (MQ=255)
                                                                                          ttttATTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc      <  2:2370668/100‑1 (MQ=255)
                                                                                          ttttACTTAGTAATATAAAATATTGAATTGTTATTTTTTTTTTTTTTTTAATATAAAAAGTCGTATTTATTATTACGTTTTTATTtaataa               >  1:683814/1‑91 (MQ=255)
                                                                                          ttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc      <  2:4177921/100‑1 (MQ=255)
                                                                                          ttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc      <  1:2675892/100‑1 (MQ=255)
                                                                                             tACTTAGTAATATAAAATTTTGCATTGTTCTTTTTGTGTGTGGTTTAAGATAAATAGCCGTCTTTATTATTACCGCTTTAATTAATAAAAGGCAGGCTGt   >  1:249382/1‑100 (MQ=255)
                                                                                             tACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGt   <  2:734934/100‑1 (MQ=255)
                                                                                              aCTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTa  >  2:1845712/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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: