Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 28.6% intergenic (‑195/+113) 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,9850TC28.6% 448.7 / 207.1 257intergenic (‑195/+113)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (92/91);  new base C (69/4);  total (161/96)
Fisher's exact test for biased strand distribution p-value = 9.59e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-03

Only 100 of 261 total aligned reads displayed.

TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705886‑3706084
                                                                                                   |                                                                                                             
tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                               >  1:70557/1‑100 (MQ=255)
 gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                              >  2:2924243/1‑100 (MQ=255)
 gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGACGTTAACGCCTTa                                                                                                              >  1:605633/1‑100 (MQ=255)
   tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                            <  1:2080456/100‑1 (MQ=255)
    gcCGGACGGCGCTGCCCGCAGAGGCACAGCGCCACCGGGTATACAGACTGTTGACACAGTGCGCTTCGTTCCTGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                           <  1:2712024/100‑1 (MQ=37)
     ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                          >  2:3311636/1‑100 (MQ=255)
      cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                           >  1:3058352/1‑98 (MQ=255)
      cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                                         >  1:2863805/1‑100 (MQ=255)
       ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                        <  2:677890/100‑1 (MQ=255)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                       <  2:2483789/100‑1 (MQ=255)
          tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                                     >  1:1392216/1‑100 (MQ=255)
           ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATaa                                                                                                    <  2:1912245/100‑2 (MQ=25)
              cctgtaggaagagGCGCAGTGCCACCGGGTATACAGATTGCTGACAAAGTGAGCTCTGATTATGCCCGATGCGTCGTTGACGCCTTATCCGGCCTAc                                                                                                    >  1:1775629/10‑97 (MQ=17)
                tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAt                                                                                              <  2:2583514/100‑1 (MQ=255)
                  cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                                            <  1:3311611/100‑1 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                                          <  2:2370937/100‑1 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                                          <  2:2634915/100‑1 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                                          >  2:2055967/1‑100 (MQ=255)
                    cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                                          <  1:2677849/100‑1 (MQ=255)
                     gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                                         <  1:3608290/100‑1 (MQ=255)
                      cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                        <  2:1528133/100‑1 (MQ=255)
                           gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                                   <  1:2930836/100‑1 (MQ=255)
                           gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                                   <  1:903891/100‑1 (MQ=255)
                            cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                                  >  2:1215504/1‑100 (MQ=255)
                            cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTCGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                                   >  1:3240987/1‑99 (MQ=255)
                              caGCGCCACCGTGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACTCCTTATCCGGCCTACAAAATCCTGCAAATTCAAt                                                                                <  1:808522/100‑1 (MQ=255)
                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                                <  2:1107637/100‑1 (MQ=255)
                               aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCTCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                                <  1:767953/99‑1 (MQ=255)
                                gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                              >  2:504571/1‑100 (MQ=255)
                                  gcCACCGAGTATGCAGATTGTTGACAAAGTGCTCTTTGTGTATGCCCGATGCGGCGTTAACGCCTTAGCAGGCCTACAAGCTCGTGCAAATTCAAt                                                                                <  1:1725959/96‑1 (MQ=17)
                                   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                           >  1:2629182/1‑100 (MQ=255)
                                     aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                                         <  2:1392222/100‑1 (MQ=255)
                                     aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                                         <  1:316794/100‑1 (MQ=255)
                                       cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAg                                                                      <  1:3509354/100‑1 (MQ=255)
                                           tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑Att                                                               >  1:1992360/1‑100 (MQ=255)
                                           tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑Att                                                               >  1:3164643/1‑100 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑g                                                      <  1:2975994/100‑1 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑g                                                      >  2:3229543/1‑100 (MQ=255)
                                               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑g                                                      <  2:410958/100‑1 (MQ=255)
                                                aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑Gt                                                     <  2:3275461/100‑1 (MQ=255)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGcc                                                >  2:1344721/1‑100 (MQ=255)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGcc                                                <  1:1002677/100‑1 (MQ=255)
                                                      ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCt                                               <  1:2213313/100‑1 (MQ=255)
                                                       tgACAAAGTGAACTTTGTTTATGCCGGATGCGGCATTAACGCCTTATCCGGCCTACAAAATCGTGAAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAgcc                                                 >  1:1433280/1‑95 (MQ=17)
                                                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATa                                           >  1:1070408/1‑100 (MQ=255)
                                                          cAAAGTGCGCTTGGTTTATGCCGGATGCGGCGTTCACGTCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATa                                           <  2:74374/100‑1 (MQ=255)
                                                               tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                      <  1:2655214/100‑1 (MQ=255)
                                                                gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                                     <  2:2290456/100‑1 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                   >  2:2311046/1‑100 (MQ=255)
                                                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                  <  2:2237953/100‑1 (MQ=255)
                                                                     atttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGA‑‑ATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCAt                               >  2:1236679/9‑100 (MQ=18)
                                                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg                            >  2:3065013/1‑100 (MQ=255)
                                                                         gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCG‑‑ATC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                            >  2:115973/5‑100 (MQ=21)
                                                                          taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCA‑‑ATT‑‑‑CTCTTGTAGGCCTGATAAGCGTAGCGCATCAgg                           >  2:3142398/4‑100 (MQ=18)
                                                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGG‑‑ACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                           >  2:1631766/4‑100 (MQ=18)
                                                                          ctggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCA‑‑ATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAg                            <  1:1322762/97‑1 (MQ=21)
                                                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGA‑‑AAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:2006942/3‑100 (MQ=18)
                                                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCATACAAAATCGTGCTAATTCAATATATTGCAGA‑‑AAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:1627716/3‑100 (MQ=18)
                                                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                         >  2:3435657/1‑100 (MQ=255)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAAC‑‑‑TT‑‑GTAGGCCTGATAAGCGCAGCGCATCAgag                          >  2:1322768/2‑98 (MQ=21)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGT‑‑TTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                         >  2:1118406/2‑100 (MQ=21)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAAT‑‑CTC‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:1641664/2‑99 (MQ=18)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGG‑‑AAA‑‑‑CAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCa                         >  1:3572448/2‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGG‑‑GAC‑‑‑TAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa                         >  1:1164148/2‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATT‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:2084962/2‑100 (MQ=17)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGA‑‑GAT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:2442943/2‑100 (MQ=17)
                                                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑GTT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                        >  2:2602012/1‑100 (MQ=255)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGA‑‑ATCC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                       >  2:365000/1‑100 (MQ=18)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCA‑‑ATT‑‑‑CTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg                        >  2:350184/1‑98 (MQ=18)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCA‑‑ATT‑‑‑CTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg                        >  1:2134425/1‑98 (MQ=18)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATG‑‑GGC‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                        >  2:913008/1‑99 (MQ=17)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCA‑‑ATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                        >  2:2106601/1‑99 (MQ=18)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGG‑‑AGT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                       >  2:2698085/1‑100 (MQ=21)
                                                                              ccggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGG‑‑GCC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                       >  2:2789724/9‑98 (MQ=17)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACC‑‑ATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                       >  1:2334547/9‑98 (MQ=21)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGG‑‑AGT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                      >  2:1014758/9‑100 (MQ=21)
                                                                              ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                       >  2:3079328/1‑100 (MQ=255)
                                                                                ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGA‑‑GTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:1898626/7‑100 (MQ=17)
                                                                                ggataaggcgTTTACGCAGCATCCGGCCTACAAAAGCATGCAAATTCGATATATTGCAGA‑‑GTT‑‑‑CAC‑‑GTAAGCCTGATAAGCGTAGCGCATAAGGCAgtt                      >  1:547477/7‑99 (MQ=17)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGT‑‑TTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                     >  2:1478204/7‑96 (MQ=25)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGT‑‑TTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                     >  1:3622922/7‑96 (MQ=25)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGA‑‑GAT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                     >  2:64034/7‑100 (MQ=17)
                                                                                 gataaggcgTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGG‑‑CTC‑‑‑GAC‑‑GTAGGCCTGATAAGCATAGAGCATCAAGCActttt                    >  1:1583303/6‑100 (MQ=18)
                                                                                 gataaggcgTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGG‑‑AAA‑‑‑CAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCAgtttt                    >  1:1604488/6‑100 (MQ=18)
                                                                                    gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                 <  1:970168/100‑1 (MQ=255)
                                                                                     cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAATTTAAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAATCATGATAAGCGGAGTGAATCAGGCATTcttgcct                >  1:2553076/1‑100 (MQ=18)
                                                                                     cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                >  1:3470182/1‑100 (MQ=255)
                                                                                      ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt               >  2:2744096/1‑100 (MQ=255)
                                                                                      ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt               >  1:3435918/1‑100 (MQ=255)
                                                                                         gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc            <  1:3014201/100‑1 (MQ=255)
                                                                                          ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc           <  2:851893/100‑1 (MQ=255)
                                                                                          ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc           <  2:359905/100‑1 (MQ=255)
                                                                                           tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa          >  1:2484278/1‑100 (MQ=255)
                                                                                             aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        >  1:2959480/1‑100 (MQ=255)
                                                                                             aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        >  1:1243904/1‑100 (MQ=255)
                                                                                              cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa       >  1:3139547/1‑100 (MQ=255)
                                                                                              cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa       >  1:1510079/1‑100 (MQ=255)
                                                                                              cGCCGTATCCGGCCTACAAAATCGTGCAAATTCAAGACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc        <  1:1051884/99‑1 (MQ=255)
                                                                                                ccTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt     >  2:1741140/1‑100 (MQ=255)
                                                                                                   tATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt  >  2:1867324/1‑100 (MQ=255)
                                                                                                   |                                                                                                             
TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGA‑‑ATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705886‑3706084

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: