Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,971 C→A 23.1% intergenic (‑181/+127) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9710CA23.1% 401.7 / 106.1 157intergenic (‑181/+127)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (57/63);  new base A (30/6);  total (87/70)
Fisher's exact test for biased strand distribution p-value = 1.17e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01

Only 100 of 170 total aligned reads displayed.

GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705873‑3706070
                                                                                                  |                                                                                                        
gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                         >  2:2129017/1‑100 (MQ=255)
gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                         >  1:2871723/1‑100 (MQ=255)
  gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                       <  1:756992/100‑1 (MQ=255)
  gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGAACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                       <  1:1226476/100‑1 (MQ=255)
   aaTGAAGTATTGCTGCCGGATGTCGCTGCACGCAGAGGCACAGCGCCCCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgccgcg                                                                                                      >  1:2057767/1‑100 (MQ=255)
   aaTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                      <  2:1710239/100‑1 (MQ=255)
   aaTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                      <  2:1545933/100‑1 (MQ=255)
    aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                     <  2:964478/100‑1 (MQ=255)
         ctATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                <  1:536707/99‑1 (MQ=255)
           aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCTACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                              <  1:2643111/100‑1 (MQ=255)
             tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                            >  2:727229/1‑100 (MQ=255)
                tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGccagatc                                                                                         >  2:571589/1‑95 (MQ=255)
                 gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                        >  2:1736738/1‑100 (MQ=255)
                 gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                        <  1:2782064/100‑1 (MQ=255)
                   cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                      <  1:984871/100‑1 (MQ=255)
                   cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                      <  2:382320/100‑1 (MQ=255)
                    ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTATGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                     <  1:1018245/100‑1 (MQ=37)
                     gacggCGCTGCACGCAGCGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTCATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                    <  1:3324301/97‑1 (MQ=255)
                      aTGGCGCTGCACGCAGAGGCACAGCGGCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                   <  2:612960/100‑1 (MQ=255)
                      aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                   <  2:602555/100‑3 (MQ=37)
                         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                               <  2:2890647/100‑1 (MQ=255)
                           gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                             >  2:1212625/1‑100 (MQ=255)
                            cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaaa                                                                            >  2:2763076/1‑100 (MQ=255)
                              gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAACGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                          >  1:3248448/1‑100 (MQ=255)
                              gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                          >  1:3209056/1‑100 (MQ=255)
                               cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                         >  2:1939794/1‑100 (MQ=255)
                               cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                         <  1:67872/100‑1 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                       <  2:3148290/100‑1 (MQ=255)
                                  gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                      <  1:2910478/100‑1 (MQ=255)
                                  gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                      <  1:2990643/100‑1 (MQ=255)
                                  gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                      >  1:3065622/1‑100 (MQ=255)
                                   cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                     >  2:3006712/1‑100 (MQ=255)
                                    agagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaa                                                                    <  2:2808960/100‑1 (MQ=255)
                                     gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaa                                                                    <  2:2148281/99‑1 (MQ=255)
                                     gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCAGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaaa                                                                   >  2:512624/1‑100 (MQ=255)
                                       ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                 <  2:2183261/100‑1 (MQ=255)
                                        gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                >  2:3310428/1‑100 (MQ=255)
                                        gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                <  1:2417294/100‑1 (MQ=255)
                                        gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                <  2:148912/100‑1 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                >  2:1081574/1‑99 (MQ=255)
                                           caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                             <  1:1981949/100‑1 (MQ=255)
                                             gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                           <  1:406756/100‑1 (MQ=255)
                                             gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                           <  1:3554152/100‑1 (MQ=255)
                                                 cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                       >  2:2532193/1‑100 (MQ=255)
                                                 cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                       >  2:44263/1‑100 (MQ=255)
                                                        tataCAGATTTTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                        <  1:512622/92‑1 (MQ=255)
                                                        tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑At‑t                                              >  1:1918168/1‑100 (MQ=255)
                                                        tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑At‑t                                              >  2:2711680/1‑100 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑g                                        >  2:2307703/1‑100 (MQ=255)
                                                              gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTa                                      >  1:2973093/1‑100 (MQ=255)
                                                                ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAgg                                    >  1:494795/1‑100 (MQ=255)
                                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGcc                                  >  2:1767517/1‑100 (MQ=255)
                                                                   ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCt                                 <  1:2013829/100‑1 (MQ=255)
                                                                     gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGa                               >  1:3534522/1‑100 (MQ=255)
                                                                      aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGAt                              >  2:1639038/1‑100 (MQ=255)
                                                                       cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCTTTCTCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATa                             <  1:2492857/100‑1 (MQ=255)
                                                                       cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATa                             <  2:2271243/100‑1 (MQ=255)
                                                                          aGTGCGCTTTGGTTATGCCGTATGCGGCGTTAACGCCTTAACCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAg                           <  1:2134437/99‑1 (MQ=255)
                                                                           gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCg                         >  2:3193447/1‑100 (MQ=255)
                                                                            tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCgg                        <  2:2713378/100‑1 (MQ=255)
                                                                              cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAg                      <  2:1789689/100‑1 (MQ=255)
                                                                              cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAg                      >  1:599879/1‑100 (MQ=255)
                                                                               gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                     <  1:2422677/100‑1 (MQ=255)
                                                                               gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAAACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                     >  2:2501397/1‑100 (MQ=255)
                                                                                tgatttgtaggCCGGATAACGCATTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCAGG‑AG‑TTAT‑‑GTAGTCCTGATAAGCGTAgcgc                   >  1:1612746/11‑100 (MQ=18)
                                                                                    ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAAAACATGCAAATTCAATAGATTGCAGAGAT‑TAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATc                >  2:3109539/7‑100 (MQ=17)
                                                                                    gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                <  1:281037/100‑1 (MQ=255)
                                                                                     tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa               <  2:3458914/100‑1 (MQ=255)
                                                                                      ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg              <  2:2811572/100‑1 (MQ=255)
                                                                                      ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg              <  1:1949751/100‑1 (MQ=255)
                                                                                       taggCCGGATAAGGCGTTTACGTCGCATCCGGCCTAC‑AAAAGAGTGCAAATTCAATACATTGCATG‑GG‑CCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  1:1260320/4‑100 (MQ=17)
                                                                                       taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGC‑‑A‑AT‑TCTCTTGTAGGCCTGATAAGCGTAGCGCATCAgg             >  2:2022583/4‑100 (MQ=18)
                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  2:665379/4‑100 (MQ=18)
                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCAAG‑AT‑TTTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  1:3007797/4‑100 (MQ=17)
                                                                                       tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg             <  1:788065/100‑1 (MQ=255)
                                                                                        aggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGCAGA‑AA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:2744517/3‑100 (MQ=18)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCACC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:943090/3‑100 (MQ=21)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAAATTGTGTAAATTCAATAAATTGATCG‑AT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:1950722/3‑100 (MQ=21)
                                                                                        aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc            <  1:2766662/100‑1 (MQ=255)
                                                                                         tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa           >  2:866177/1‑100 (MQ=255)
                                                                                         tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa           >  1:3083606/1‑100 (MQ=255)
                                                                                         tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATCCAATACATCGCAGA‑AT‑TTTC‑‑GTAGGCCTGACAAGCGGAGCGCATCAGGCa           >  1:694511/1‑100 (MQ=255)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCACC‑AT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:3343404/2‑100 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCAGG‑AG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:2520086/2‑100 (MQ=17)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:1269527/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:2990578/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCAGG‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:2343883/2‑99 (MQ=21)
                                                                                          gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt          >  1:781198/1‑100 (MQ=255)
                                                                                          gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt          >  1:1659209/1‑100 (MQ=255)
                                                                                          gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt          >  1:1510454/1‑100 (MQ=255)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt         >  2:2264169/1‑100 (MQ=18)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt         >  1:1808540/1‑100 (MQ=18)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCAGG‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt         >  2:1424177/1‑100 (MQ=21)
                                                                                           ccGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCAGA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  2:3378991/1‑98 (MQ=17)
                                                                                             ggATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt       >  2:1332835/1‑100 (MQ=255)
                                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATTCATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  2:1380178/100‑1 (MQ=255)
                                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  1:1101177/100‑1 (MQ=255)
                                                                                                  cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  <  2:2557941/100‑1 (MQ=255)
                                                                                                  cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  >  2:1382781/1‑100 (MQ=255)
                                                                                                  cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  <  1:1201001/100‑1 (MQ=255)
                                                                                                  |                                                                                                        
GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705873‑3706070

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: ‑