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,9460CT99.0% 622.9 / ‑2.9 202intergenic (+106/+110)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (94/102);  total (99/103)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

Only 100 of 203 total aligned reads displayed.

TGAAGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTA  >  NC_000913/720847‑721042
                                                                                                   |                                                                                                
tgaaGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTc                                                                                                  <  1:1687598/100‑1 (MQ=255)
   aGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGt                                                                                               >  2:887368/1‑100 (MQ=255)
    ggTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTATTTAGt                                                                                               >  2:3339908/1‑99 (MQ=255)
      tGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAat                                                                                            <  2:2105035/100‑1 (MQ=255)
            tatTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATaaaa                                                                                      >  2:39789/1‑100 (MQ=255)
             atTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAat                                                                                     <  1:4073434/100‑1 (MQ=255)
               tCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGATTTTTATTTAGTAATATaa                                                                                        >  1:1392969/1‑95 (MQ=255)
                 tATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTg                                                                                 <  1:1219590/100‑1 (MQ=255)
                 tATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTg                                                                                 <  1:789416/100‑1 (MQ=255)
                 tATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTg                                                                                 >  1:620283/1‑100 (MQ=255)
                    ccAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAt                                                                              >  2:669201/1‑100 (MQ=255)
                     cAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAAtt                                                                             >  1:2118296/1‑100 (MQ=255)
                         aaaTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTa                                                                         <  1:4757/100‑1 (MQ=255)
                         aaaTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAACTGtt                                                                          >  2:1079027/1‑99 (MQ=255)
                            tACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAttt                                                                      >  2:3034083/1‑100 (MQ=255)
                             aCGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAtttt                                                                     <  2:1197207/100‑1 (MQ=255)
                              cGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAttttt                                                                    <  1:2667708/100‑1 (MQ=255)
                               ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAttttt                                                                    >  2:1737783/1‑99 (MQ=255)
                               ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                   <  2:280275/100‑1 (MQ=255)
                               ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                   <  1:742963/100‑1 (MQ=255)
                                gTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgt                                                                  >  2:1827432/1‑100 (MQ=255)
                                 ttGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtg                                                                 <  1:1963750/100‑1 (MQ=255)
                                 ttGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtg                                                                 <  2:1161084/100‑1 (MQ=255)
                                 tccgatcttgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTGAGTAATATAAAATATTGAATTGTTATTTTtgtg                                                                 <  2:2503275/92‑1 (MQ=255)
                                  tGAACTGTTAAGGTGATAAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtg                                                                 <  2:3200724/99‑1 (MQ=255)
                                  tGAACTGTGAAGGTTATCAATAAAAAATGATCAATTTTATTTTATTTAAGGATGTGCTTTTTACTTGGTAATATAAAATATTGAAGTGGTATTTTtg                                                                   <  2:580511/97‑1 (MQ=255)
                                  tGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtgt                                                                >  1:3128346/1‑100 (MQ=255)
                                   gAACTGTGAAGTTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgt                                                                  <  2:4103151/97‑1 (MQ=255)
                                    aaCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTgtgtgt                                                              >  1:642189/1‑100 (MQ=255)
                                    aaCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTgtgtgt                                                              <  1:18162/100‑1 (MQ=255)
                                     attgtgAAGTTGATCAATAAAAAATGATCAATCATAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTgtgtgt                                                              <  2:740691/97‑1 (MQ=255)
                                     aCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGtgtt                                                             <  1:2066674/100‑1 (MQ=255)
                                      cTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGtgttg                                                            <  1:1998208/100‑1 (MQ=255)
                                      cTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGtgttg                                                            >  1:3011659/1‑100 (MQ=255)
                                      cTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGtgtt                                                             <  1:2781150/99‑1 (MQ=255)
                                      cTGTGAAGGTGATCAATAAAAAATAATCAATTTTAATTTTTTTAATGATGATCTTTTAACTAATTAAAATAAAATATTGATTTTTAATTTTTttgggt                                                              >  2:1017874/1‑98 (MQ=255)
                                       tgtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGtgttgt                                                           >  1:1424375/1‑100 (MQ=255)
                                        gtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAttttttttttttttt                                                          >  2:2597984/1‑90 (MQ=255)
                                        gtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTgttgtt                                                          >  2:2394096/1‑100 (MQ=255)
                                        gtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTgttgtt                                                          >  2:3840944/1‑100 (MQ=255)
                                         tgAAGTTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGttt                                                         <  2:2550031/100‑1 (MQ=255)
                                         tgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtgagatcgga                                                         >  1:2503290/1‑92 (MQ=255)
                                         tgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGttt                                                         >  1:426913/1‑100 (MQ=255)
                                          gAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTa                                                        >  2:942283/1‑100 (MQ=255)
                                          gAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTa                                                        <  1:1153819/100‑1 (MQ=255)
                                          gAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTATGTGTTTTGTTTa                                                        >  2:4136317/1‑100 (MQ=255)
                                            aGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAg                                                      <  1:826548/100‑1 (MQ=255)
                                            aGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAg                                                      >  2:4027003/1‑100 (MQ=255)
                                              gTGATCAATAAAAAATTATCAATCTTAATTTATTTAATGATGAGCTTTTTTCTTAGTAATATAAAATATTGAATTGTCATTTTTGTGTGTTGTTTAAGa                                                     <  2:3945365/99‑1 (MQ=255)
                                               tGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGGGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATa                                                   >  1:15853/1‑100 (MQ=255)
                                               tGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATa                                                   >  1:1049031/1‑100 (MQ=255)
                                               tGATCAATAAAAAATGATCAATCTTAAATTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATa                                                   <  2:1726791/100‑1 (MQ=255)
                                                gATCAATAAAAAATGATCCATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTTAATTGTTATTTTTGTGGGTGGTTTATGATa                                                   >  2:22378/1‑99 (MQ=255)
                                                 aTCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaa                                                 >  1:3953087/1‑100 (MQ=255)
                                                  actaTAAAAAATGTTCAATCTTAATTTATAAAATGATGAACTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTTTTTTGTTTCAGATaaaa                                                >  2:1610317/4‑100 (MQ=255)
                                                   cAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaaaa                                               >  2:1479427/1‑100 (MQ=255)
                                                    aaTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAg                                              >  1:1325624/1‑100 (MQ=255)
                                                     aTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTTTTGTTTAAGATAAAAAGc                                             >  2:3641032/1‑100 (MQ=255)
                                                     aTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaaaa                                               >  2:545234/1‑98 (MQ=255)
                                                     aTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGATTAAGATAAAAAGc                                             >  2:1379379/1‑100 (MQ=255)
                                                         aaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGta                                         >  2:30861/1‑100 (MQ=255)
                                                           aaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatt                                       <  1:1795252/100‑1 (MQ=255)
                                                             tGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtattta                                     <  1:1891065/100‑1 (MQ=255)
                                                               aTCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATGTAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                   <  2:2669709/100‑1 (MQ=255)
                                                               aTCAATCTTAATTTATTTAATGATGAGCTTTTTACGTAGTAATATAAAATATTGAATTGTTTTTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                   >  2:805863/1‑100 (MQ=255)
                                                                tCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttatta                                  <  2:627062/100‑1 (MQ=255)
                                                                   aTCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattatta                               <  2:2975127/100‑1 (MQ=255)
                                                                     cTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACg                             <  1:3714593/100‑1 (MQ=255)
                                                                     cTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACg                             >  1:3693006/1‑100 (MQ=255)
                                                                     cTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACg                             >  2:1518841/1‑100 (MQ=255)
                                                                      ttAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACgg                            <  1:2540568/100‑1 (MQ=255)
                                                                       tAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                           <  1:3543151/100‑1 (MQ=255)
                                                                       tAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTAAAAAAAAAAAAAAAAAGCCGTATTTATTATTACGGc                           <  2:1406600/100‑1 (MQ=255)
                                                                        aaTTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt                          >  2:1770756/1‑100 (MQ=255)
                                                                          tttatttaATGATGGTTTTTTTTCTTAGTATTATAAATGTTTGAATGGTTATTTTTGTGTGTGGTTTGAGATAAAAAGCCGTATTTATTATTACGGcgt                         <  2:2896699/99‑3 (MQ=255)
                                                                           ttatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttt                       <  2:2136883/100‑2 (MQ=38)
                                                                            tatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGAAGTATTTATTATTACGGCTttaa                      >  2:1403992/1‑100 (MQ=255)
                                                                             atttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaat                     >  1:973196/1‑100 (MQ=255)
                                                                             atttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaat                     >  1:3842885/1‑100 (MQ=255)
                                                                             atttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaa                      >  2:2853887/1‑99 (MQ=255)
                                                                             atttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttttta                     <  2:2314547/100‑4 (MQ=37)
                                                                              tttaATGATGGGATTTTTACTTAGTAAGAGAAAATATTGAATTGTTATGTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                    <  2:35552/100‑1 (MQ=255)
                                                                              tttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                    <  2:2096633/100‑1 (MQ=255)
                                                                              tttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                    <  2:524391/100‑1 (MQ=255)
                                                                               ttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatta                   <  1:3451644/100‑1 (MQ=255)
                                                                               ttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAAATAAAAAGCCGTAATAATAAATACGGCtttttatc                   <  2:3685456/100‑6 (MQ=25)
                                                                                taATGATGAGCTTTTTACTTATTAATATAAAATACTGAATTGTTATTTTTATCTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                  <  2:1004764/100‑7 (MQ=21)
                                                                                taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                  <  1:221804/100‑7 (MQ=25)
                                                                                taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                  <  1:3169585/100‑7 (MQ=25)
                                                                                taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                  <  1:3684934/100‑7 (MQ=25)
                                                                                taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                  <  2:3256385/100‑7 (MQ=25)
                                                                                  atgatgAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctta                <  1:2518391/100‑9 (MQ=25)
                                                                                     ntgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaaa             >  2:1231978/2‑100 (MQ=255)
                                                                                     atgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATaaaa             >  1:470264/1‑100 (MQ=255)
                                                                                      tgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAag            >  1:1255730/1‑100 (MQ=255)
                                                                                          cTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag        >  1:2474747/1‑100 (MQ=255)
                                                                                           tttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcagg       <  1:2751100/100‑1 (MQ=255)
                                                                                              ttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTg    <  2:274231/100‑1 (MQ=255)
                                                                                               tACTGAGTAATATAAAATGTTGAATTGGTATTTTGGTGTTTTGTGTAATATAAAAAGGCGTATTTATTGTTATGGGTTTAATTAATAAAAGGCAGGCTGt   >  2:2428256/1‑100 (MQ=255)
                                                                                                aCTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTa  >  2:3875925/1‑100 (MQ=255)
                                                                                                   |                                                                                                
TGAAGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTA  >  NC_000913/720847‑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: