Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,971 C→A 21.3% 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,9710CA21.3% 620.1 / 147.9 231intergenic (‑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 (85/97);  new base A (44/5);  total (129/102)
Fisher's exact test for biased strand distribution p-value = 1.57e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-02

Only 100 of 243 total aligned reads displayed.

GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705873‑3706070
                                                                                                  |                                                                                                         
gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                          >  2:2404183/1‑100 (MQ=255)
gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                          >  1:738494/1‑100 (MQ=255)
  gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                        >  1:3275436/1‑100 (MQ=255)
  gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                        <  1:2987211/100‑1 (MQ=255)
  gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                        <  1:1965877/100‑1 (MQ=255)
    aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                      >  2:2939286/1‑100 (MQ=255)
      gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                    <  1:2097709/100‑1 (MQ=255)
      gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                    <  1:86004/100‑1 (MQ=255)
       aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                   <  1:3577241/100‑1 (MQ=255)
       aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                   >  1:1002174/1‑100 (MQ=255)
        aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                   >  2:3462030/1‑99 (MQ=255)
        aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                  <  1:252699/100‑1 (MQ=255)
         gTATTGCCGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                 >  2:1940974/1‑100 (MQ=255)
          tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                 <  2:2532002/99‑1 (MQ=255)
          tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                                <  1:3422502/100‑1 (MQ=255)
           aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCAGCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                                <  1:3209110/99‑1 (MQ=255)
            ttGCTGCCGGATGTCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                              <  1:602182/100‑1 (MQ=255)
            ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                              <  2:1642827/100‑1 (MQ=255)
            ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                              >  1:2651115/1‑100 (MQ=255)
            ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                              <  1:2744593/100‑1 (MQ=255)
            ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                              <  2:208250/100‑1 (MQ=255)
            ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                              >  1:1525855/1‑100 (MQ=255)
             tgcggcCGGATTGCGCTGCACGGAGAGGCACAGCGCCACCGGGTATAAAGATTGTTGACAAAGTTCGATTTATTTATGCCAGATGCGGAGTTAACggatt                                                                                             >  1:2296347/1‑96 (MQ=25)
                tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                          <  1:2080456/100‑1 (MQ=255)
                 gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                         <  2:2095945/100‑1 (MQ=255)
                   cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                         >  1:3058352/1‑98 (MQ=255)
                    ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                      <  2:677890/100‑1 (MQ=255)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                     <  2:1620752/100‑1 (MQ=255)
                          cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaac                                                                                <  2:2210724/100‑7 (MQ=25)
                             tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAt                                                                            <  2:2583514/100‑1 (MQ=255)
                              gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                           <  2:441404/100‑1 (MQ=255)
                               cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                          >  1:1604397/1‑100 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                        <  2:2634915/100‑1 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                        >  2:2055967/1‑100 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                        <  2:2370937/100‑1 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                        >  2:260996/1‑100 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                        <  1:3120527/100‑1 (MQ=255)
                                  gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                       <  1:3608290/100‑1 (MQ=255)
                                   cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                      >  1:3561829/1‑100 (MQ=255)
                                   cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                      >  2:2339947/1‑100 (MQ=255)
                                      agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                   <  2:1633978/100‑1 (MQ=255)
                                       ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                                  <  2:3226075/100‑1 (MQ=255)
                                        gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                 <  1:1387367/100‑1 (MQ=255)
                                        gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                 <  1:2930836/100‑1 (MQ=255)
                                        gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                 <  1:903891/100‑1 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                >  2:3608372/1‑100 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                >  2:1215504/1‑100 (MQ=255)
                                          acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                               >  2:1905814/1‑100 (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                                                       <  1:316794/100‑1 (MQ=255)
                                                     gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGCCCTACAAAATCGTGCAAATTCAATACATTGCAgc                                                    >  1:2570598/1‑99 (MQ=255)
                                                        tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAt‑‑t                                              >  1:3164643/1‑100 (MQ=255)
                                                        tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAt‑‑t                                              >  2:3545010/1‑100 (MQ=255)
                                                           aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTc                                           <  1:2951912/100‑1 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑g                                        <  2:27927/100‑1 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑g                                        <  1:660342/100‑1 (MQ=255)
                                                             aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑Gt                                       <  1:338780/100‑1 (MQ=255)
                                                             aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑Gt                                       >  2:2750626/1‑100 (MQ=255)
                                                             aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑Gt                                       <  1:98090/100‑1 (MQ=255)
                                                             aGATTGTTGAATAAGTGCGCTTTGTTTATGCCGGATTCGGCGTTAACGCCTTATCCGGACTACAAAATCGTGCAAATTCAATACCTTGCAGAAT‑‑TTTC‑‑Gt                                       <  1:1663792/100‑1 (MQ=17)
                                                               aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAg                                     <  2:1525862/100‑1 (MQ=255)
                                                                   ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCt                                 <  1:1838717/100‑1 (MQ=255)
                                                                   ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCt                                 >  2:3422694/1‑100 (MQ=255)
                                                                       cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATa                             <  1:2688472/100‑1 (MQ=255)
                                                                       cAAAGTGCGCTTGGTTTATGCCGGATGCGGCGTTCACGTCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATa                             <  2:74374/100‑1 (MQ=255)
                                                                           gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCg                         <  2:3276946/100‑1 (MQ=255)
                                                                              tccgatctttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCATAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAg                      <  2:2254921/93‑1 (MQ=255)
                                                                               gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                     >  2:2311046/1‑100 (MQ=255)
                                                                                  atttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCAt                 >  2:1236679/9‑100 (MQ=18)
                                                                                     tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa               >  2:1459496/1‑100 (MQ=255)
                                                                                      ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg              >  2:3065013/1‑100 (MQ=255)
                                                                                      gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCa               >  2:873352/5‑100 (MQ=18)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGAT‑‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  2:115973/5‑100 (MQ=21)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAG‑‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  2:2676870/5‑100 (MQ=21)
                                                                                       taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGC‑‑AAT‑‑TCTCTTGTAGGCCTGATAAGCGTAGCGCATCAgg             >  2:3142398/4‑100 (MQ=18)
                                                                                       tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg              <  1:1409824/99‑1 (MQ=255)
                                                                                        aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAA‑‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:2006942/3‑100 (MQ=18)
                                                                                        aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGT‑‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:1558833/3‑100 (MQ=17)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGAT‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:1026973/3‑100 (MQ=17)
                                                                                         ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt          >  1:1062888/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGAC‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:1375208/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTT‑‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa           >  2:1118406/2‑100 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGAT‑‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:1284612/2‑100 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCT‑‑CCGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:3254941/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGA‑‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa           >  1:1164148/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGA‑‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:2442943/2‑100 (MQ=17)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt         >  2:365000/1‑100 (MQ=18)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGC‑‑AAT‑‑TCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg          >  1:2134425/1‑98 (MQ=18)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGG‑‑CCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa          >  2:913008/1‑99 (MQ=17)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAG‑‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  2:1328885/1‑99 (MQ=17)
                                                                                           ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt         >  2:3079328/1‑100 (MQ=255)
                                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  1:2233374/100‑1 (MQ=255)
                                                                                                tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc    <  1:2659298/100‑1 (MQ=255)
                                                                                                 gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc   <  1:970168/100‑1 (MQ=255)
                                                                                                  cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  >  2:1755694/1‑100 (MQ=255)
                                                                                                  cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  >  1:3470182/1‑100 (MQ=255)
                                                                                                  |                                                                                                         
GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705873‑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: