Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,483,912 T→A 18.0% V53D (GTC→GAC)  evgA → response regulator in two‑component regulatory system with EvgS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,483,9120TA18.0% 401.7 / 76.0 166V53D (GTC→GAC) evgAresponse regulator in two‑component regulatory system with EvgS
Reads supporting (aligned to +/- strand):  ref base T (72/64);  new base A (15/15);  total (87/79)
Fisher's exact test for biased strand distribution p-value = 8.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 171 total aligned reads displayed.

TGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTT  >  NC_000913/2483813‑2484011
                                                                                                   |                                                                                                   
tGAGCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCTGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGt                                                                                                     <  1:480606/100‑1 (MQ=255)
 gATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTc                                                                                                    >  2:1121331/1‑100 (MQ=255)
 gATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTc                                                                                                    >  2:570329/1‑100 (MQ=255)
  aTCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCg                                                                                                   >  1:283132/1‑100 (MQ=255)
   tCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGa                                                                                                  >  1:653952/1‑100 (MQ=255)
   tCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGa                                                                                                  >  1:1373463/1‑100 (MQ=255)
   tCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGa                                                                                                  <  1:1450513/100‑1 (MQ=255)
    cAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGat                                                                                                 >  2:44491/1‑100 (MQ=255)
     aaaaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGata                                                                                                <  1:1733791/100‑1 (MQ=255)
     aaaaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGata                                                                                                <  2:2777764/100‑1 (MQ=255)
      aaaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGatat                                                                                               >  2:188921/1‑100 (MQ=255)
      aaaaCGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGatat                                                                                               <  2:600341/100‑1 (MQ=39)
         aCGATATTTAAATCTTAGCATATTTGACTGAAGGCGTAAGTTCCGTTCAGCGGTTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATccc                                                                                            <  1:1994173/100‑1 (MQ=255)
          cGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATcccc                                                                                           <  1:1452477/100‑1 (MQ=255)
          cGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcagg                                                                                           <  2:2374495/100‑4 (MQ=38)
          cGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcagg                                                                                           <  1:2072330/100‑4 (MQ=38)
          cGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcag                                                                                            <  2:1466193/99‑3 (MQ=38)
           gATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCg                                                                                          >  1:3020944/1‑100 (MQ=255)
           gATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCg                                                                                          >  2:3923896/1‑100 (MQ=255)
           gATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcaggc                                                                                          <  2:3578902/100‑5 (MQ=38)
           gATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcaggc                                                                                          <  2:3429799/100‑5 (MQ=38)
           gATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcagg                                                                                           <  1:4120332/99‑4 (MQ=38)
            atatTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCgg                                                                                         <  2:2018603/100‑1 (MQ=255)
            atatTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCgg                                                                                         <  2:3344098/100‑1 (MQ=255)
            atatTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCgg                                                                                         <  1:3141465/100‑1 (MQ=255)
             tatTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGa                                                                                        <  1:3249274/100‑1 (MQ=255)
               ttGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGt                                                                                      >  2:2704428/1‑100 (MQ=255)
                  aaaTCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTaa                                                                                   >  1:2251376/1‑100 (MQ=255)
                  aaaTCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTaa                                                                                   >  1:2250052/1‑100 (MQ=255)
                   aaTCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAAc                                                                                  <  1:920296/100‑1 (MQ=255)
                     tCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCCGCGGGTGGAAACCCATAAGCCTGATATCCTCATCATTGCTGTCGATATtccacgagctga                                                                                   >  1:4357082/1‑88 (MQ=255)
                     tCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACgg                                                                                >  1:355048/1‑100 (MQ=255)
                      cTTATCAGAATGAACCGAAGCCTGAAGTGCCGTTTAGCTTTGGCAAACACTTAAGCCTGATATCGTCATTATTGACGTCGCGACCCCTGGAGTTAACgg                                                                                <  1:175868/99‑1 (MQ=255)
                       ttAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTa                                                                              <  1:62047/100‑1 (MQ=255)
                        tAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGCGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTAt                                                                             >  1:4039336/1‑100 (MQ=255)
                          gCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATGGATGTCGAAATCCCCGGAGTTAACGGTATcc                                                                           >  1:2975756/1‑100 (MQ=255)
                              agTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGgt                                                                       >  2:3052302/1‑100 (MQ=255)
                              agTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGgt                                                                       <  1:2157843/100‑1 (MQ=255)
                                  gACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTa                                                                   <  2:3605719/100‑1 (MQ=255)
                                   aCTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAg                                                                  <  1:3923866/100‑1 (MQ=255)
                                   aCTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAg                                                                  <  2:2282730/100‑1 (MQ=255)
                                    cTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGa                                                                 >  2:1980794/1‑100 (MQ=255)
                                     tGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGaa                                                                <  1:1493715/100‑1 (MQ=255)
                                      gAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATACAGGTGTTAGaaa                                                               >  1:1445479/1‑100 (MQ=255)
                                       aaGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAAc                                                              <  1:2285511/100‑1 (MQ=255)
                                        aGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACg                                                             >  2:741869/1‑100 (MQ=255)
                                        aGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACg                                                             >  2:3330080/1‑100 (MQ=255)
                                         ggCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGc                                                            <  2:134253/100‑1 (MQ=255)
                                          gCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGca                                                           <  1:4221566/100‑2 (MQ=255)
                                           cGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGc                                                            >  1:3630314/1‑98 (MQ=255)
                                           cGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTg                                                          <  1:2384172/100‑1 (MQ=255)
                                                 tGCGGTTGAGCGGGTGGAAACACTTAATCCGGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAg                                                        <  2:3997176/96‑1 (MQ=255)
                                                  gCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAgc                                                   >  2:384804/1‑100 (MQ=255)
                                                  gCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAgc                                                   >  1:1707150/1‑100 (MQ=255)
                                                     gTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGcc                                                <  2:2696772/100‑1 (MQ=255)
                                                     gTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGcc                                                <  2:693287/100‑1 (MQ=255)
                                                        cAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAat                                             >  2:606770/1‑100 (MQ=255)
                                                          gCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAatat                                           >  1:3992759/1‑100 (MQ=255)
                                                           cGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAAtata                                          >  2:1075580/1‑100 (MQ=255)
                                                             ngTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCACCAATATAGc                                        >  2:3579285/2‑100 (MQ=255)
                                                               tGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCgg                                      <  2:2286629/100‑1 (MQ=255)
                                                                     cacTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattat                                >  2:3357920/1‑100 (MQ=255)
                                                                       cTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattatta                              <  2:2130994/100‑1 (MQ=255)
                                                                        ttAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattattat                             <  2:4269884/100‑1 (MQ=255)
                                                                         ctccggggATATCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattattatt                            >  1:3741614/8‑100 (MQ=37)
                                                                         ctccggggATATCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAattattatt                            <  2:3741640/93‑1 (MQ=37)
                                                                           ccggggATATCGACATCAATGATGACGATATCTCCGGAGTTAACGGTATCCAGGTGTTAGAAACGATGAGGAAGCGCCAATATAGGTGAATtattattat                          >  1:13924/6‑100 (MQ=25)
                                                                           aGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATtattattat                          >  2:2223639/1‑100 (MQ=255)
                                                                            gCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTCTCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGAGGAattatcattatc                         >  1:3496286/1‑100 (MQ=255)
                                                                            cggggATATCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATc                         >  2:2415965/5‑100 (MQ=38)
                                                                            cggggATATCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATc                         >  2:3535612/5‑100 (MQ=38)
                                                                             ccTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCg                        <  2:2577172/100‑1 (MQ=255)
                                                                             ccTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCg                        <  2:573871/100‑1 (MQ=255)
                                                                             ccTGATATCGTCATCATTGATGTCGATATCAACGGAGTTAACGGTATCCAGGTGTTAGAAACGCAGAGGAAGCGCCAATATAGCGGAATTATTATTATc                         >  1:1991469/1‑99 (MQ=255)
                                                                               tGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGtc                      <  1:2543781/100‑1 (MQ=255)
                                                                               tGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGtc                      >  1:4301731/1‑100 (MQ=255)
                                                                               ggATATCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGtc                      >  1:4107186/2‑100 (MQ=39)
                                                                               ggATATCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGACAATATAGCGGAATTATTATTATCGtc                      >  1:3169344/2‑100 (MQ=38)
                                                                               agATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGAGTTAGAAACGCTGAGGAGGCGACAATATAGTGGAATTATTATTATCg                        >  1:2146136/2‑98 (MQ=255)
                                                                                  tatCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTcc                   <  1:2086068/100‑1 (MQ=255)
                                                                                  tatCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTcc                   <  2:4184874/100‑1 (MQ=255)
                                                                                  tatCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGt                       >  1:106585/1‑96 (MQ=39)
                                                                                   atCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCg                  >  1:2225538/1‑100 (MQ=255)
                                                                                     cGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCGCCGCt                >  2:380213/1‑100 (MQ=255)
                                                                                      gTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTa               >  2:3568458/1‑100 (MQ=255)
                                                                                           catTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAAt          >  2:3027173/1‑100 (MQ=255)
                                                                                           catTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAAt          >  2:2162472/1‑100 (MQ=255)
                                                                                           catTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAAt          >  1:2674968/1‑100 (MQ=255)
                                                                                            atTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATg         <  2:2839445/100‑1 (MQ=255)
                                                                                            atTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATg         >  2:1095404/1‑100 (MQ=255)
                                                                                             ttGATGTCGGTATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGTTGCGCCAATATAGCGGAATTATTATTATCATCTCCGCTAAAAATGa        >  1:2571167/1‑100 (MQ=255)
                                                                                             ttGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGa        <  2:323491/100‑1 (MQ=255)
                                                                                              tGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGAc       <  2:2340881/100‑1 (MQ=255)
                                                                                              tGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGAc       <  1:1276052/100‑1 (MQ=255)
                                                                                              tGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGAc       <  2:984622/100‑1 (MQ=255)
                                                                                               gATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGAcc      <  1:3397650/100‑1 (MQ=255)
                                                                                               gATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGAcc      >  2:2110010/1‑100 (MQ=255)
                                                                                               gATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGAcc      <  2:4073395/100‑1 (MQ=255)
                                                                                                  gTCGATATCCCCGGAGTTACCGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCAtt   <  1:1064802/100‑1 (MQ=255)
                                                                                                   tCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCAttt  >  1:1962883/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTT  >  NC_000913/2483813‑2484011

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: