Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 720,946 C→T 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,9460CT93.7% 282.8 / 6.0 95intergenic (+106/+110)ybfK/kdpEuncharacterized protein/response regulator in two‑component regulatory system with KdpD
Reads supporting (aligned to +/- strand):  ref base C (2/4);  new base T (46/43);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

AAGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTA  >  NC_000913/720849‑721042
                                                                                                 |                                                                                                
aaGGTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAg                                                                                                <  1:1833589/100‑1 (MQ=255)
 aagTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAg                                                                                                <  1:180013/97‑1 (MQ=255)
  ggTGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTa                                                                                              >  1:1398404/1‑100 (MQ=255)
    tGAAGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAat                                                                                            >  2:1100572/1‑100 (MQ=255)
      aaGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAatat                                                                                          >  1:1887290/1‑100 (MQ=255)
       aGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCGATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAataaa                                                                                         <  2:1318097/100‑3 (MQ=255)
       aGATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAAtata                                                                                         >  2:1762954/1‑100 (MQ=255)
       aGATATTCTATCCAGGAAATACGGTTAAACTGTTAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAatat                                                                                          <  1:1317688/99‑1 (MQ=255)
        gATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATaa                                                                                        >  2:913230/1‑100 (MQ=255)
        gATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATaa                                                                                        <  1:737775/100‑1 (MQ=255)
        gATATTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATaa                                                                                        <  1:1405555/100‑1 (MQ=255)
           atTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAat                                                                                     >  1:1624423/1‑100 (MQ=255)
           atTCTATCCAGGAAATACGGTTGAACTGTGAAGGTGACTAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAat                                                                                     >  1:505612/1‑100 (MQ=255)
            ttCTATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAata                                                                                    >  2:748710/1‑100 (MQ=255)
               tATCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTg                                                                                 >  1:981901/1‑100 (MQ=255)
                aTCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGa                                                                                >  1:1294730/1‑100 (MQ=255)
                aTCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGa                                                                                >  1:1447977/1‑100 (MQ=255)
                aTCCAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGa                                                                                >  1:1538481/1‑100 (MQ=255)
                  ccAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAt                                                                              >  1:949835/1‑100 (MQ=255)
                   cAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAtt                                                                             >  1:1839621/1‑100 (MQ=255)
                   cAGGAAATACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAtt                                                                             >  2:1230357/1‑100 (MQ=255)
                       aaaTACGGTTGAACTGTGAAGGTGATCAATAAAAAATGAGCAAACTTAATTTATTTAATGATGAGCTTTTTACGTAGTAATATGAAAGATTGAATTg                                                                            <  1:1972/97‑1 (MQ=255)
                          tACGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAttt                                                                      >  1:1418304/1‑100 (MQ=255)
                           aCGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAtttt                                                                     <  2:608895/100‑1 (MQ=255)
                            cGGTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTAttttt                                                                    >  2:625712/1‑100 (MQ=255)
                             ggTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                   <  1:45575/100‑1 (MQ=255)
                              gTTGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtg                                                                   <  1:1701260/99‑1 (MQ=255)
                               ttGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtg                                                                 <  1:423832/100‑1 (MQ=255)
                               ttGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtg                                                                 <  2:661463/100‑1 (MQ=255)
                               ttGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTAAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAACTATTGAATTGTTATTTTtgtg                                                                 <  1:1551302/100‑1 (MQ=255)
                                tGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTtgtgt                                                                >  2:971046/1‑100 (MQ=255)
                                tGAACTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTtgtgt                                                                <  2:875109/100‑1 (MQ=255)
                                  aaCTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGATTTGTTATTTTTgtgtgt                                                              >  1:401553/1‑100 (MQ=255)
                                    cTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGtgttg                                                            >  2:947540/1‑100 (MQ=255)
                                    cTGTGAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGtgttg                                                            >  2:1706665/1‑100 (MQ=255)
                                     tgtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGtgttgt                                                           >  1:257355/1‑100 (MQ=255)
                                      gtgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTgttgtt                                                          >  1:1933949/1‑100 (MQ=255)
                                       tgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGttt                                                         >  2:174202/1‑100 (MQ=255)
                                       tgAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGttt                                                         >  2:1915504/1‑100 (MQ=255)
                                        gAAGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTa                                                        >  1:1350085/1‑100 (MQ=255)
                                         aaGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTaa                                                       >  1:796043/1‑100 (MQ=255)
                                         aaGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTaa                                                       >  2:241686/1‑100 (MQ=255)
                                          aGGTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAg                                                      <  2:1187515/100‑1 (MQ=255)
                                          aGGTGATAAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTaa                                                       <  1:508492/99‑1 (MQ=255)
                                            gTGATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGAt                                                    >  1:52961/1‑100 (MQ=255)
                                              gATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaa                                                  >  2:249251/1‑100 (MQ=255)
                                              gATCAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAAGTGTTATTTTTGTGTGTTGTTTAAGATaa                                                  >  2:1462001/1‑100 (MQ=255)
                                                 cAATAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATaaaaa                                               >  1:1118455/1‑100 (MQ=255)
                                                  aaTAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGGATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAg                                              >  2:7426/1‑100 (MQ=255)
                                                    tAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGcc                                            <  1:601101/100‑1 (MQ=255)
                                                    tAAAAAATGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGcc                                            <  2:1118456/100‑1 (MQ=255)
                                                     naaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCg                                           >  2:827454/2‑100 (MQ=255)
                                                     aaaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCg                                           <  1:748709/100‑1 (MQ=255)
                                                     aaaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCg                                           <  1:1230353/100‑1 (MQ=255)
                                                      aaaaaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGt                                          <  2:981902/100‑1 (MQ=255)
                                                         aaTGATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatt                                       <  1:1635270/100‑1 (MQ=255)
                                                            gATCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttat                                    <  1:209391/100‑1 (MQ=255)
                                                             aTCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGtatttatt                                   >  2:1259405/1‑100 (MQ=255)
                                                             aTCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGATGTTTAAGATAAAAAGCCGtatttat                                    >  1:1743597/1‑99 (MQ=255)
                                                              tCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttatta                                  <  2:69122/100‑1 (MQ=255)
                                                              tCAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttatta                                  <  1:1358242/100‑1 (MQ=255)
                                                              tCAATCTTAATTTATTTAATGATGAGCTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttatta                                  <  2:694960/100‑1 (MQ=255)
                                                               cAATCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattat                                 >  2:297128/1‑100 (MQ=255)
                                                                aaTCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATttattatt                                >  2:85010/1‑100 (MQ=255)
                                                                  tCTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTAc                              <  2:233711/100‑1 (MQ=255)
                                                                   cTTAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACg                             <  1:1100571/100‑1 (MQ=255)
                                                                    ttAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTTTTGTTTAAAATAAAAAGCCGTATttattatta                               >  1:1812006/1‑97 (MQ=255)
                                                                    ttAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACgg                            <  1:1257456/100‑1 (MQ=255)
                                                                     tAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                           >  2:1822758/1‑100 (MQ=255)
                                                                     tAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                           <  1:1539832/100‑1 (MQ=255)
                                                                     tAATTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGc                           <  2:1263855/100‑1 (MQ=255)
                                                                      aaTTTATTTAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCt                          >  2:1414315/1‑100 (MQ=255)
                                                                       atttatttAATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCtt                         <  1:1455830/100‑1 (MQ=255)
                                                                          tatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaa                      <  1:525258/100‑1 (MQ=255)
                                                                          tatttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaa                      >  1:773405/1‑100 (MQ=255)
                                                                           atttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaat                     >  1:1666346/1‑100 (MQ=255)
                                                                             ttaATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTttaatta                   <  1:1360350/100‑1 (MQ=255)
                                                                              taATGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTAATAATAAATACGGCtttttatct                  <  2:1866639/100‑7 (MQ=25)
                                                                               aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                 <  2:1624433/100‑1 (MQ=255)
                                                                               aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                 >  2:1591569/1‑100 (MQ=255)
                                                                               aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                 <  1:1152350/100‑1 (MQ=255)
                                                                               aaTGATGAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTtaattaat                 <  2:819347/100‑1 (MQ=255)
                                                                                atgatgAGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATtaata                <  1:1414307/100‑1 (MQ=255)
                                                                                    tgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAag            <  1:1093356/100‑1 (MQ=255)
                                                                                    tgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAag            <  1:1380823/100‑1 (MQ=255)
                                                                                    tgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAag            <  1:77218/100‑1 (MQ=255)
                                                                                    tgaGCTTTTTACTTAGTAATATAAAATATTGAATTGTGATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTANTTAATAAAag            <  2:668852/100‑1 (MQ=255)
                                                                                       gCTTTTTACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca         >  2:129686/1‑100 (MQ=255)
                                                                                        cTTTTTACTCAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcag        <  1:1915491/100‑1 (MQ=255)
                                                                                         tttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcagg       >  2:1802024/1‑100 (MQ=255)
                                                                                         tttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggca         <  2:1685074/98‑1 (MQ=255)
                                                                                          ttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc      >  1:1466968/1‑100 (MQ=255)
                                                                                          ttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc      <  2:503146/100‑1 (MQ=255)
                                                                                          ttttACTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAaggcaggc      <  2:1900421/100‑1 (MQ=255)
                                                                                              aCTTAGTAATATAAAATATTGAATTGTTATTTTTGTGTGTTGTTTAAGATAAAAAGCCGTATTTATTATTACGGCTTTAATTAATAAAAGGCAGGCTGTa  >  1:397802/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 < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: