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% 487.7 / ‑3.1 151intergenic (+106/+110)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (64/83);  total (65/86)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 151 total aligned reads displayed.

AAGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTA  >  NC_000913/720849‑721042
                                                                                                 |                                                                                                
aaGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAg                                                                                                <  2:1681733/100‑1 (MQ=255)
 aGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGt                                                                                               <  1:1467243/100‑1 (MQ=255)
 aGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGATAATAAAAAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGt                                                                                               <  2:1365033/100‑1 (MQ=255)
   gTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGCTCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTaa                                                                                             >  2:2666212/1‑100 (MQ=255)
     gAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAata                                                                                           <  2:3759514/100‑1 (MQ=255)
      aaGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAatat                                                                                          <  1:3625956/100‑1 (MQ=255)
        gATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATaa                                                                                        <  2:1324807/100‑1 (MQ=255)
         atatTCTATCCAGGAAATATGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGTTGAGCTTTTTACTTAGTAAtata                                                                                         >  2:2508033/1‑98 (MQ=255)
          tatTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATaaa                                                                                       >  2:1858398/1‑99 (MQ=255)
           actcgaTCCAGAAAATACGGTGGCACTGTGAAGGTGATGAATAAAAAATGATCTATCTTAATTTATTAAACGATGAGTTTTTTACTAAGTAATATaaa                                                                                       >  2:2699205/6‑98 (MQ=255)
            ttCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAata                                                                                    >  1:3156284/1‑100 (MQ=255)
            ttCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAata                                                                                    >  2:2064508/1‑100 (MQ=255)
             tCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAACTTTTTACTTAGTAATATAAAatat                                                                                   >  2:3911582/1‑100 (MQ=255)
              cTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATAtt                                                                                  >  1:3335860/1‑100 (MQ=255)
               tATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTg                                                                                 <  2:630890/100‑1 (MQ=255)
                 tCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGaa                                                                               >  1:1280795/1‑100 (MQ=255)
                  ccAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAt                                                                              >  1:1600335/1‑100 (MQ=255)
                  ccAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAt                                                                              >  1:2051559/1‑100 (MQ=255)
                  ccAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAt                                                                              <  1:2223094/100‑1 (MQ=255)
                  ccAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAt                                                                              >  1:2576114/1‑100 (MQ=255)
                  ccAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATCAGCTTTTTACTTAGTAATATAAAATATTGAAt                                                                              >  2:2148885/1‑100 (MQ=255)
                   cAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAtt                                                                             >  2:1324779/1‑100 (MQ=255)
                   cAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAtt                                                                             >  1:779304/1‑100 (MQ=255)
                    aGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTg                                                                            <  1:788307/100‑1 (MQ=255)
                     ggAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAAAGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGt                                                                           >  1:1004468/1‑100 (MQ=255)
                       aaaTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGtt                                                                          <  2:2922024/99‑1 (MQ=255)
                        aaTACGGTTGAACTGTTAAGTTGATCAATAAAAAATTATCAATGTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAt                                                                        <  2:1405208/100‑1 (MQ=255)
                        aaTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAt                                                                        >  1:1666180/1‑100 (MQ=255)
                             ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAttttt                                                                    >  2:2820640/1‑99 (MQ=255)
                             ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                   <  1:2666228/100‑1 (MQ=255)
                             ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                   <  1:3864777/100‑1 (MQ=255)
                              gTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgt                                                                  <  1:2284939/100‑1 (MQ=255)
                              gTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgt                                                                  >  1:339465/1‑100 (MQ=255)
                              gTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgt                                                                  >  2:2249249/1‑100 (MQ=255)
                              gTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                   <  2:2089511/99‑1 (MQ=255)
                               ttGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtg                                                                 <  2:541785/100‑1 (MQ=255)
                               ttGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtg                                                                 <  1:184981/100‑1 (MQ=255)
                                tGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgttt                                                                >  1:1546642/1‑98 (MQ=255)
                                tGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtgt                                                                <  1:3227607/100‑1 (MQ=255)
                                 gAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtgtg                                                               <  1:1917367/100‑1 (MQ=255)
                                  aaCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTgtgtgt                                                              <  1:1349906/100‑1 (MQ=255)
                                   aCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGtgtt                                                             >  1:4124524/1‑100 (MQ=255)
                                    cTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGtgttg                                                            >  1:1476796/1‑100 (MQ=255)
                                      gtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTgttgtt                                                          <  1:1884285/100‑1 (MQ=255)
                                       ttaaGGTGGTTAATAAATTATGATTAATCTTAATTTATTTAATGATGTGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGtgtggt                                                           <  2:3598450/96‑1 (MQ=255)
                                        gAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTa                                                        <  1:1279471/100‑1 (MQ=255)
                                        gAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTa                                                        >  2:3859794/1‑100 (MQ=255)
                                           ggTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGa                                                     <  2:3069354/100‑1 (MQ=255)
                                           ggTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGa                                                     >  1:3280590/1‑100 (MQ=255)
                                           ggTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGa                                                     >  2:1156979/1‑100 (MQ=255)
                                             tGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATa                                                   >  1:2688371/1‑100 (MQ=255)
                                              gATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaa                                                  >  2:2343837/1‑100 (MQ=255)
                                              gATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaa                                                  >  1:1167610/1‑100 (MQ=255)
                                              gATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGAATTTTACTTAGTAATATAAAATATTGAATTGTTATGTTTGTGTGTTGTTTAAGATa                                                   >  2:3622415/1‑99 (MQ=255)
                                              gATCAATAAAAAATGATCAATCTTAATTTATTTAATGATAAGCTTTTTACTTAGTAATATAAAATATTGAGGTGTTATTTTTGTGTGTTGTTTAAGATaa                                                  >  2:4120701/1‑100 (MQ=255)
                                                 cAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaaaa                                               >  2:3262237/1‑100 (MQ=255)
                                                   aTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAgn                                             <  2:2096888/100‑2 (MQ=255)
                                                   aTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGc                                             >  1:4136703/1‑100 (MQ=255)
                                                   aTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGc                                             <  2:1280787/100‑1 (MQ=255)
                                                   aTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGc                                             <  1:3225157/100‑1 (MQ=255)
                                                    tAAAAAATGATCAATCTTAATTTATTTAATTATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGcc                                            <  2:4068978/100‑1 (MQ=255)
                                                    tAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGcc                                            >  1:3289767/1‑100 (MQ=255)
                                                     aaaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCg                                           <  2:3994256/100‑1 (MQ=255)
                                                     aaaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCg                                           <  1:3699913/100‑1 (MQ=255)
                                                     aaaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCg                                           <  1:3155534/100‑1 (MQ=255)
                                                      aaaaaTGATCAATCTTAATTTATTTAATGATGAGCAAATTAATTAGTAATATAAAATATTAAATTGTTATATTTGTGAGTTGTTTAAAATAAAAAACCGt                                          >  2:661450/1‑100 (MQ=255)
                                                       aaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGta                                         <  2:3143503/100‑1 (MQ=255)
                                                       aaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGta                                         >  1:2944430/1‑100 (MQ=255)
                                                        aaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtat                                        >  1:3699712/1‑100 (MQ=255)
                                                           tGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtattta                                     <  1:3046891/100‑1 (MQ=255)
                                                             aTCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                   <  1:2323442/100‑1 (MQ=255)
                                                                 aTCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattatta                               <  2:1114260/100‑1 (MQ=255)
                                                                     tAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                           >  1:1322392/1‑100 (MQ=255)
                                                                      aaTTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTAAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCt                          <  2:2051546/100‑1 (MQ=255)
                                                                       atttatttAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                         <  1:874371/100‑1 (MQ=255)
                                                                         ttatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttt                       <  1:1786357/100‑2 (MQ=38)
                                                                           atttaATGATGAGCTTTTTACTGAGTAATATAAAATATTGAATTGTTATTGTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCttttta                     <  2:1530453/100‑4 (MQ=37)
                                                                            tttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatt                    <  1:951674/100‑1 (MQ=255)
                                                                            tttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttan                    <  2:1559050/100‑5 (MQ=37)
                                                                             ttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatta                   <  2:3393693/100‑1 (MQ=255)
                                                                              taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                  <  1:2045797/100‑7 (MQ=25)
                                                                              taATCATGTGCTTTTTACATAGTAATATAAAATATTTATTTGTTGTTTTTGTGTGTTGTGTCTGATACAGAGCCCTACTTTTTATTACGGCTttaattga                  <  2:3024111/100‑3 (MQ=255)
                                                                               aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCTTATTTATTATTACGGCTTtaattaat                 <  1:2513465/100‑1 (MQ=255)
                                                                               aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                 <  2:3094922/100‑8 (MQ=25)
                                                                               aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatctt                 <  2:3321099/100‑8 (MQ=25)
                                                                                atgattAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTTTGTGTTGTTTAAGATAAAAAGCCGTATTTATTAGTACGGCTTTAATtagta                <  2:3910556/100‑1 (MQ=255)
                                                                                atgatgAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                >  1:1155716/1‑100 (MQ=255)
                                                                                atgatgAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                <  1:3841850/100‑1 (MQ=255)
                                                                                  gatgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatcttaaa              <  2:1751511/100‑11 (MQ=21)
                                                                                    tgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCAGTATTTATTATTACGGCTTTAATTAATAAAag            >  2:56313/1‑100 (MQ=255)
                                                                                     gaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg           >  1:2474767/1‑100 (MQ=255)
                                                                                      aGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTTTTTTTTGTTTAATATAAAAAGCCGTATTTATTATTAATGCTTTAATTAATAAAag            >  2:2146071/1‑98 (MQ=255)
                                                                                      aGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAagg           <  2:2952240/99‑1 (MQ=255)
                                                                                       gCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca         <  1:3636910/100‑1 (MQ=255)
                                                                                        cTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag        <  2:1645189/100‑1 (MQ=255)
                                                                                        cTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag        >  1:1553393/1‑100 (MQ=255)
                                                                                         tttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcagg       <  2:779298/100‑1 (MQ=255)
                                                                                          ttttACTTAGTAATATAAAATATTGCATCGCTATTTTTGTGTTCTTTTAAAGATAAAAAGCCGTCTTTATTATTACTGCTTTAATTAATAAAaggcagtc      >  2:2573159/1‑98 (MQ=255)
                                                                                             tACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGt   <  1:1031481/100‑1 (MQ=255)
                                                                                              aCTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTa  >  2:3994749/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: