Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,971 C→A 20.6% 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,9710CA20.6% 388.5 / 81.3 146intergenic (‑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 (52/64);  new base A (24/6);  total (76/70)
Fisher's exact test for biased strand distribution p-value = 8.32e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.33e-01

Only 100 of 155 total aligned reads displayed.

AGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705872‑3706070
                                                                                                   |                                                                                                       
 gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                        <  2:955674/100‑1 (MQ=255)
 gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                        >  1:775867/1‑100 (MQ=255)
  ggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                       >  2:1879488/1‑100 (MQ=255)
    aaTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                     <  1:1408871/100‑1 (MQ=255)
    aaTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                     <  1:2182734/100‑1 (MQ=255)
      tGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                                   <  2:2684959/100‑1 (MQ=255)
       gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                  <  1:1089955/100‑1 (MQ=255)
           tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                               <  2:595133/99‑1 (MQ=255)
            aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                             <  1:2117774/100‑1 (MQ=255)
            aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                             <  2:2313047/100‑1 (MQ=255)
            aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                             >  2:1300874/1‑100 (MQ=255)
                 atctggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                        <  2:321559/96‑1 (MQ=255)
                  gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                       <  1:347443/100‑1 (MQ=255)
                   ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                      <  1:2476364/100‑1 (MQ=255)
                    cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                      >  1:2545208/1‑99 (MQ=255)
                     ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                    <  2:1714430/100‑1 (MQ=255)
                     ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                    <  2:54136/100‑1 (MQ=37)
                      gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                   <  2:748352/100‑1 (MQ=255)
                      gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                   <  1:2031458/100‑1 (MQ=255)
                       aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                  <  1:1561210/100‑3 (MQ=37)
                         ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                                <  1:1707507/100‑1 (MQ=255)
                          gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTata                                                                               <  1:1866528/100‑3 (MQ=255)
                          gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                              >  1:1003293/1‑100 (MQ=255)
                           cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaac                                                                              <  2:913562/100‑7 (MQ=25)
                            gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                            >  2:2176277/1‑100 (MQ=255)
                            gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                            <  1:1148944/100‑1 (MQ=255)
                             cTGCACGCAGAGGCACAGCGCCACCGTGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaac                                                                              <  2:1234556/98‑7 (MQ=21)
                             cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACACCTTATCCGGCCTAC‑aa                                                                             >  2:214648/1‑98 (MQ=255)
                              tGCACGCAGAGGCACAGCTCCAGCGGGTATACAGATTGTTGACATAGTGCGTCTTGCTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                          <  2:1619637/100‑1 (MQ=21)
                              tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                          >  2:2032110/1‑100 (MQ=255)
                              tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                          >  1:571058/1‑100 (MQ=255)
                                  cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                      >  2:2527939/1‑100 (MQ=255)
                                  cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                      >  2:2691398/1‑100 (MQ=255)
                                    cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                    >  2:652063/1‑100 (MQ=255)
                                        ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                >  2:6353/1‑100 (MQ=255)
                                         gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                               <  1:327576/100‑1 (MQ=255)
                                         gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                               <  1:2665575/100‑1 (MQ=255)
                                          cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                              >  1:1308065/1‑100 (MQ=255)
                                          cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                              >  2:666547/1‑100 (MQ=255)
                                          cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                              <  1:1962945/100‑1 (MQ=255)
                                          cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGACTTATCCGGCCTAC‑AAAATCGTGCAACTTc                                                               >  2:200813/1‑99 (MQ=255)
                                            caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                            <  1:652065/100‑1 (MQ=255)
                                              gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACTCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                          <  1:1941499/100‑1 (MQ=255)
                                                  cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                      <  1:1547195/100‑1 (MQ=255)
                                                  cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                      <  1:1694706/100‑1 (MQ=255)
                                                     cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAg                                                   <  2:768951/100‑1 (MQ=255)
                                                       ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAACTTCAATACATTGCAGaa                                                 >  2:1046636/1‑100 (MQ=255)
                                                           taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAt‑ttt                                            >  2:206616/1‑100 (MQ=255)
                                                           taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAt‑ttt                                            >  2:241697/1‑100 (MQ=255)
                                                            aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTc                                           >  1:1065956/1‑100 (MQ=255)
                                                             cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑g                                        <  1:666549/100‑1 (MQ=255)
                                                             cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑g                                        <  1:1605600/100‑1 (MQ=255)
                                                             cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑g                                        >  1:266256/1‑100 (MQ=255)
                                                              aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑g                                        <  2:2545190/99‑1 (MQ=255)
                                                               gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTa                                      >  2:46913/1‑100 (MQ=255)
                                                                  tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGc                                   <  1:699052/100‑1 (MQ=255)
                                                                   gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGcc                                  <  2:2429716/100‑1 (MQ=255)
                                                                   gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGc                                   >  1:1061909/1‑99 (MQ=255)
                                                                    ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCt                                 <  1:554829/100‑1 (MQ=255)
                                                                     tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTg                                <  1:300445/100‑1 (MQ=255)
                                                                     tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTg                                >  1:1220658/1‑100 (MQ=255)
                                                                      gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGa                               >  2:2443537/1‑100 (MQ=255)
                                                                      gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGa                               >  1:407172/1‑100 (MQ=255)
                                                                       aCAAAGTGTTCTTTGTTGATGCCGGAGGCGGCGTTAACGCCTTATCCGGCCTAC‑ATAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGAt                              <  2:606249/100‑1 (MQ=18)
                                                                        cacaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATa                             >  2:1821333/4‑100 (MQ=255)
                                                                        cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATa                             <  2:1750035/100‑1 (MQ=255)
                                                                           aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGc                          >  2:2485976/1‑100 (MQ=255)
                                                                           aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGc                          <  1:2450090/100‑1 (MQ=255)
                                                                           aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGc                          <  1:8254/100‑1 (MQ=255)
                                                                            gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCg                         <  1:2591667/100‑1 (MQ=255)
                                                                             tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGGATAT‑AAAATCGTTCAGATTCAATACATTGCAGAAT‑TTAC‑‑GTAGTCCTGAt                              <  2:667154/94‑1 (MQ=17)
                                                                             tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCgg                        <  1:501753/100‑1 (MQ=255)
                                                                                gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                     >  1:300917/1‑100 (MQ=255)
                                                                                   ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgtg                    <  1:1842417/98‑3 (MQ=255)
                                                                                      tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa               <  1:455788/100‑1 (MQ=255)
                                                                                       ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg              <  1:1105509/100‑1 (MQ=255)
                                                                                        taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCAGAGT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  1:534579/4‑100 (MQ=17)
                                                                                        tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg             >  1:2365747/1‑100 (MQ=255)
                                                                                         aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGC‑‑AAT‑TCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:1379932/3‑100 (MQ=18)
                                                                                         aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCAGTTT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc            >  2:101518/3‑100 (MQ=21)
                                                                                         aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGTTAATTCAATATATTGCAGGGA‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc            >  2:6433/3‑100 (MQ=18)
                                                                                         aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCAGGAG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:2052756/3‑100 (MQ=21)
                                                                                         aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc            >  1:1507962/1‑100 (MQ=255)
                                                                                         aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc            >  1:1507151/1‑100 (MQ=255)
                                                                                          tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa           >  2:848606/1‑100 (MQ=255)
                                                                                          ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCAGAAA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:287713/2‑99 (MQ=17)
                                                                                          ggCCGGATAAGGCGTTCACTCCGCATCCGGCCTAC‑AAAAGCACGCAAATTCACTATATTGCAGAGA‑TTGC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa           >  2:1641264/2‑100 (MQ=11)
                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:2158119/2‑100 (MQ=18)
                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:1350640/2‑100 (MQ=18)
                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCAGTTT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa           >  1:2152905/2‑100 (MQ=21)
                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑AAT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:1624189/2‑100 (MQ=18)
                                                                                          agCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt          >  1:225696/2‑100 (MQ=18)
                                                                                           gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCAGAAA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  1:1011896/1‑99 (MQ=18)
                                                                                           gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGC‑‑AAT‑TCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg          >  1:2395303/1‑98 (MQ=18)
                                                                                           gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCAGAGT‑TCAC‑‑GTAGGCCTGTTAAGCGTAGCGCATCAGGCag          >  2:537208/1‑99 (MQ=17)
                                                                                            ccGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCAGGGA‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCagc         >  2:2251430/1‑98 (MQ=18)
                                                                                                aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  2:144561/100‑1 (MQ=255)
                                                                                                 tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc    <  1:1879496/100‑1 (MQ=255)
                                                                                                 tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc    <  2:2415588/100‑1 (MQ=255)
                                                                                                   cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  >  2:1534336/1‑100 (MQ=255)
                                                                                                   |                                                                                                       
AGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705872‑3706070

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑