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% 440.9 / 124.0 182intergenic (‑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 (69/71);  new base A (33/9);  total (102/80)
Fisher's exact test for biased strand distribution p-value = 7.48e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.14e-01

Only 100 of 191 total aligned reads displayed.

GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705873‑3706070
                                                                                                  |                                                                                                          
gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                           >  1:3034722/1‑100 (MQ=255)
     tGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                                      <  2:71224/100‑1 (MQ=255)
     tGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                                      <  1:3045970/100‑1 (MQ=255)
      gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                     >  1:1229326/1‑100 (MQ=255)
       aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                     >  1:3186369/1‑99 (MQ=255)
         gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                   <  2:1888034/99‑1 (MQ=255)
         gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                  >  1:3673142/1‑100 (MQ=255)
           aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                <  1:3240886/100‑1 (MQ=255)
              gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                             >  1:721500/1‑100 (MQ=255)
              gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                             >  1:2170085/1‑100 (MQ=255)
               ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                            >  1:1952097/1‑100 (MQ=255)
               ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                            >  1:1994346/1‑100 (MQ=255)
                 gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                          >  2:590897/1‑100 (MQ=255)
                 gacGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                          >  2:2579282/3‑100 (MQ=255)
                  ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                         <  1:3371274/100‑1 (MQ=255)
                   cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                        <  2:2024779/100‑1 (MQ=255)
                   cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCgg                                                                                        <  1:2364407/100‑1 (MQ=37)
                   cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCgg                                                                                        <  2:2199478/100‑1 (MQ=37)
                    ggATGGCGCTGCACGCAGAGGCACAGCGCTACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                       <  1:367745/100‑1 (MQ=25)
                    ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                       <  2:654093/100‑1 (MQ=37)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                      >  2:2642354/1‑100 (MQ=255)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                      <  1:1439088/100‑2 (MQ=37)
                     gATGGCGCGGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTGTGGTTATGCCGGATGCGGCGTTATCGCCTTATCCGGcc                                                                                      <  1:1279210/100‑1 (MQ=25)
                       tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATa                                                                                    <  1:1804698/100‑1 (MQ=37)
                       tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATa                                                                                    <  2:3423502/100‑1 (MQ=37)
                         gcgcTGCACGCAGAGGCACAGCGTCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                 <  1:511058/100‑1 (MQ=255)
                         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATCCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                 <  1:2011271/100‑1 (MQ=255)
                         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCAGCCTAC‑a                                                                                 >  1:363104/1‑100 (MQ=255)
                           gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCTGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                               >  1:935600/1‑100 (MQ=255)
                              gCACGCAGAGGCACAGCGCCACCGTATATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                             >  1:1924241/1‑99 (MQ=255)
                               cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                           <  2:573253/100‑1 (MQ=255)
                                aCGCAGAGGCACAGCGCCACCGGGTATACAGACTGTTGAGAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑agaa                                                                              >  1:3529737/1‑96 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                         >  1:2084578/1‑100 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTg                                                                         >  2:3473030/1‑100 (MQ=21)
                                  gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                        <  1:2676683/100‑1 (MQ=255)
                                  gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                        <  2:1254493/100‑1 (MQ=255)
                                  gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                        <  2:749898/100‑1 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCa                                                                <  2:3103257/100‑1 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCa                                                                <  2:857784/100‑1 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCa                                                                >  1:2202194/1‑100 (MQ=255)
                                           caacgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTGTGCCGGATGCGGCGTTAACACCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCa                                                                >  1:2379309/4‑98 (MQ=255)
                                            aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATa                                                             <  2:1920368/100‑1 (MQ=255)
                                             gagcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATAc                                                            <  2:2470253/98‑1 (MQ=255)
                                               gcCACTGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACAt                                                          >  2:1213499/1‑100 (MQ=255)
                                                ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACAtt                                                         <  2:2182920/100‑1 (MQ=255)
                                                ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGACGGATGCGGCGTTAACGCCTTATTCGGCCTAC‑AAAATCGTGC‑AAATTCAATACAtt                                                         >  1:1894181/1‑100 (MQ=255)
                                                  aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGc                                                       <  2:1647400/100‑1 (MQ=255)
                                                  aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGc                                                       <  2:247163/100‑1 (MQ=255)
                                                      ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGaa                                                 <  2:1264040/100‑1 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑g                                        <  2:492490/100‑1 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑g                                        <  2:1989021/100‑1 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑g                                        >  2:3218815/1‑100 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑g                                        <  1:2800359/100‑1 (MQ=255)
                                                                ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAgg                                    <  1:1917468/100‑1 (MQ=255)
                                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGcc                                  <  2:2951956/100‑1 (MQ=255)
                                                                   ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCt                                 >  2:1823407/1‑100 (MQ=255)
                                                                    tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTg                                >  1:2712378/1‑100 (MQ=255)
                                                                    tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTg                                >  1:3438451/1‑100 (MQ=255)
                                                                       cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATa                             >  1:3196531/1‑100 (MQ=255)
                                                                          aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGc                          <  1:3751427/100‑1 (MQ=255)
                                                                            tGCGCGTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCgg                        <  1:1307701/100‑1 (MQ=255)
                                                                               gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGTCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                     >  1:2611423/1‑100 (MQ=255)
                                                                               gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGTTAAGCg                         >  2:2423031/1‑97 (MQ=255)
                                                                               gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGa                       <  1:2325865/98‑1 (MQ=255)
                                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTGC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAg                      >  1:857020/1‑98 (MQ=255)
                                                                                cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                    <  1:1419714/100‑1 (MQ=255)
                                                                                 cttttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAAGCATGC‑AAATTCAATATATTGC‑‑AGAGT‑TCAC‑‑GTAGGCCTGATAAGCGTAgcgc                   >  2:2007855/10‑100 (MQ=17)
                                                                                  ttttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAATCGTGC‑TAATTCAATATATTGC‑‑AGAAA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCa                  >  2:1505621/9‑100 (MQ=18)
                                                                                  ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                  >  2:3356334/1‑100 (MQ=255)
                                                                                  ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                  >  2:301389/1‑100 (MQ=255)
                                                                                  ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                  <  1:838261/100‑1 (MQ=255)
                                                                                  ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                  <  2:1240099/100‑1 (MQ=255)
                                                                                   tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGC‑CAATTCAATATATTGC‑‑AGGAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCAt                 >  1:2726716/8‑100 (MQ=18)
                                                                                   tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTAT‑AAATTCAATGTATTGC‑‑AGTTT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCAt                 >  2:2698219/8‑100 (MQ=21)
                                                                                   tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                 >  1:1198358/1‑100 (MQ=255)
                                                                                    ttgtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCC‑AAATTCAATATATTGC‑‑AGAAA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATc                >  2:2307935/7‑100 (MQ=18)
                                                                                    gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                >  1:2451771/1‑100 (MQ=255)
                                                                                      gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAATTC‑TGC‑AAATTCAACAAATTGC‑‑‑GAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  1:1548395/5‑100 (MQ=18)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTT‑AAATTCAATATATTGC‑‑ACCAT‑CCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  2:1563847/5‑100 (MQ=21)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGA‑AAATTCAATATATTGC‑‑AGGAG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  2:2838790/5‑100 (MQ=18)
                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTAT‑AAATTCAATGTATTGC‑‑AGTTT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg             >  2:2889443/4‑100 (MQ=25)
                                                                                       tATGCCGGATGCGGCGTTAACGCCTTATCCGGCTTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg             <  1:3496329/100‑1 (MQ=255)
                                                                                        nggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑ATACCGCGC‑AAATTCAATATATTGCGGAGAAA‑AT‑‑‑‑GTAGGCCTGATAAGCGTAGCGCGTCAGGct           >  1:2043341/3‑99 (MQ=18)
                                                                                        aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGC‑AAATTCAACATATTGC‑‑AATTC‑TCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  2:1063249/3‑99 (MQ=17)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGC‑AAATTCAATATATTGC‑‑AGGAA‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGc            >  2:172828/3‑100 (MQ=18)
                                                                                         tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa           <  2:406035/100‑1 (MQ=255)
                                                                                         tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa           >  1:1379054/1‑100 (MQ=255)
                                                                                         ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGC‑AAATTCAATATATTGC‑‑AGAGT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:386876/2‑100 (MQ=14)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGC‑CAATTCAATATATTGC‑‑AGGAC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:1389603/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGC‑AAATTCAATATATTGC‑‑‑‑AAT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:2116612/2‑100 (MQ=18)
                                                                                          gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt          >  2:1460892/1‑100 (MQ=255)
                                                                                          gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt          >  1:2885851/1‑100 (MQ=255)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCC‑AAATTCAATATATTGC‑‑AGAAA‑TCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa          >  2:2975010/1‑99 (MQ=18)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGC‑AAATTCAATAAATTGC‑‑A‑AAT‑CAACATGTAGGCCGGATAAGCGCAGCGCATCAGGCa           >  1:2571710/1‑100 (MQ=18)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGC‑TAATTCAATATATTGT‑‑GGGCT‑CCAC‑‑GTAGGCCTGATAAGCGTAGCTCATCAGGCag          >  1:2497197/1‑99 (MQ=18)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGC‑TAATTCAATATATTGC‑‑AGGGA‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCag          >  1:1217899/1‑99 (MQ=18)
                                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  2:2232917/100‑1 (MQ=255)
                                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  1:2341624/100‑1 (MQ=255)
                                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  1:3025232/100‑1 (MQ=255)
                                                                                                  cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  <  1:3742334/100‑1 (MQ=255)
                                                                                                  |                                                                                                          
GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGC‑AAATTCAATACATTGC‑‑AGAAT‑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: ‑