Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,971 C→A 26.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,9710CA26.3% 407.2 / 136.2 186intergenic (‑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 (61/76);  new base A (42/7);  total (103/83)
Fisher's exact test for biased strand distribution p-value = 3.20e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

Only 100 of 196 total aligned reads displayed.

AGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705872‑3706070
                                                                                                   |                                                                                                        
 gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                         >  2:1748949/1‑100 (MQ=255)
   gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGGATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                       <  1:147863/100‑1 (MQ=255)
      tGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                                    <  1:1649064/100‑1 (MQ=255)
          gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCTCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                <  1:559728/100‑1 (MQ=255)
          gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                >  1:680080/1‑100 (MQ=255)
            aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                              <  1:2400577/100‑1 (MQ=255)
               gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                             >  1:1340217/1‑98 (MQ=255)
               gcggcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                           <  1:1974711/97‑1 (MQ=255)
                 tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                         <  1:1258918/100‑1 (MQ=255)
                  gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                        >  1:481410/1‑100 (MQ=255)
                  gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATca                                                                                        >  2:1258923/1‑99 (MQ=255)
                  gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCCTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                        >  1:1397792/1‑100 (MQ=255)
                   ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGccgcatccg                                                                                       <  2:1512444/100‑8 (MQ=255)
                    cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                      >  2:438106/1‑100 (MQ=255)
                     ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                     <  2:2319675/100‑1 (MQ=255)
                      gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                    <  2:839367/100‑1 (MQ=255)
                      gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                    <  2:680083/100‑1 (MQ=255)
                      gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                    <  1:2060932/100‑2 (MQ=37)
                        tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                  >  2:223994/1‑100 (MQ=255)
                           cgcTGCACGCAGAGGCACAGCGCCAGCGGGTATACAGATTGTGGACAATGTGCGCTTTGTTTATGCCGGAAGCGGCGTTAACGCCTTATCCGGACTACa                                                                               >  1:2249205/1‑99 (MQ=25)
                            gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                             <  2:2143171/100‑1 (MQ=255)
                            gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaaca                                                                              <  1:1596315/100‑8 (MQ=25)
                             cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaacaa                                                                             <  2:1539336/100‑9 (MQ=25)
                              tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaagt                                                                           >  1:1770645/1‑98 (MQ=255)
                              tGCACGCACAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCTGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                             <  1:1911095/98‑1 (MQ=255)
                               gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                          >  1:1431514/1‑100 (MQ=255)
                               gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                          >  2:1597978/1‑100 (MQ=255)
                                cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                         >  2:11003/1‑100 (MQ=255)
                                cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                         <  2:638664/100‑1 (MQ=255)
                                cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                         <  2:1124983/100‑1 (MQ=255)
                                 aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGt                                                                        <  2:373007/100‑1 (MQ=255)
                                 aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGAGTTAACGCCTTATCCGGCATACAAAATCGt                                                                        >  1:594317/1‑100 (MQ=255)
                                  cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                       >  2:2710107/1‑100 (MQ=255)
                                  cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                       <  2:2109035/100‑1 (MQ=255)
                                   gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                      <  1:1515062/100‑1 (MQ=255)
                                     agagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaa                                                                    <  1:745594/100‑1 (MQ=255)
                                      nagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                                                   >  2:292453/2‑100 (MQ=255)
                                      gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                                                   <  1:1084670/100‑1 (MQ=255)
                                        ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                                 >  2:212578/1‑100 (MQ=255)
                                          cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                <  1:1305625/99‑1 (MQ=255)
                                          cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                               <  2:2103742/100‑1 (MQ=255)
                                            caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                             >  2:2791515/1‑100 (MQ=255)
                                            caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTAGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATGGTGCAAATTCa                                                               >  1:1228583/1‑98 (MQ=255)
                                                 ccACCGGGTATAAAGATTGTTGACAAAGTGCGCTTTGTTTCTGACAGATTCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGTAAATTCCATACAta                                                        <  1:1736197/100‑2 (MQ=17)
                                                      gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑Ag                                                  >  2:1586866/1‑99 (MQ=255)
                                                          ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAt                                               >  1:1897524/1‑98 (MQ=255)
                                                           taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAtttt                                            >  1:470312/1‑100 (MQ=255)
                                                           taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAAtttt                                            >  2:1487690/1‑100 (MQ=255)
                                                            aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTc                                           <  1:2169490/100‑1 (MQ=255)
                                                             cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑g                                        <  2:1868313/100‑1 (MQ=255)
                                                             cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑g                                        <  2:2439752/100‑1 (MQ=255)
                                                             cAGATTGTTGACAAAGTGCGCTTTGGTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑g                                        <  1:2484867/100‑1 (MQ=255)
                                                              aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑Gt                                       >  2:2511059/1‑100 (MQ=255)
                                                               gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTa                                      <  1:920925/100‑1 (MQ=255)
                                                               gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTa                                      <  1:1719291/100‑1 (MQ=255)
                                                               gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTa                                      >  2:814308/1‑100 (MQ=255)
                                                               gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTa                                      <  2:186109/100‑1 (MQ=255)
                                                                ttccgatctaaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAg                                     <  1:1174405/91‑1 (MQ=255)
                                                                 ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAgg                                    <  1:2342527/100‑1 (MQ=255)
                                                                   gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGcc                                  <  1:105792/100‑1 (MQ=255)
                                                                      gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTg                                <  1:2197914/100‑1 (MQ=255)
                                                                        cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATa                             >  2:1687996/1‑100 (MQ=255)
                                                                        cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATa                             <  2:922526/100‑1 (MQ=255)
                                                                         aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAgagatcggaa                            >  2:1174410/1‑91 (MQ=255)
                                                                             tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCgg                        <  1:1544549/100‑1 (MQ=255)
                                                                             tGCGCTTTGTTTATGACGGATGCGGCGTTAACGCCTTATCCTGCCTACAAAATCGTGCAAATTCTATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCg                         <  1:835977/99‑1 (MQ=255)
                                                                               cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGa                       >  1:2699127/1‑99 (MQ=255)
                                                                                  cttttgtaggCCGGATAAGGCGTTCGCGCCGCATCCGGCCTACAAAAGTTTGCAAATTCAATAAATTGC‑‑AGAATTCAT‑‑GTAGGCCTGATAAGCGAAgcgc                   >  1:2650431/10‑100 (MQ=18)
                                                                                    tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                 <  2:1546142/100‑1 (MQ=255)
                                                                                     ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGC‑‑AATATGAC‑‑GTAGGCCTGATAAGCGTAGCGCATCa               <  2:2810056/94‑1 (MQ=18)
                                                                                     ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGC‑‑AAGATTTTC‑‑GTAGGCCTGATAAGCGTAGCGCATc                >  2:1850542/7‑100 (MQ=17)
                                                                                     ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGC‑‑AGGATACGT‑‑GTAGGCCTGATAAGCGTAGCGCATc                >  2:2034384/7‑100 (MQ=18)
                                                                                      tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa               <  2:972948/100‑1 (MQ=255)
                                                                                       ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg              <  2:689519/100‑1 (MQ=255)
                                                                                       ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg              >  1:2879080/1‑100 (MQ=255)
                                                                                       ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg              >  2:290374/1‑100 (MQ=255)
                                                                                        taggCCGGATAAGGCGTTCACTCCGCATCCGGCCTACAAAAGCACGCAAATTCAATATATTGC‑‑AGAGATTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  1:2210726/4‑100 (MQ=11)
                                                                                        taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGC‑‑AGGAGCTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  2:2277299/4‑100 (MQ=17)
                                                                                        tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGACTGATAAGCGGAGCGCATCAgg             >  1:1303748/1‑100 (MQ=255)
                                                                                         aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:2229127/3‑100 (MQ=18)
                                                                                         aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGC‑‑AGAGTTCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:1831243/3‑100 (MQ=17)
                                                                                         aggCCGGATAAGGCGTTCACTCCGCATCCGGCCTACAAAAGCACGCAAATTCAATATATTGC‑‑AGAGATTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:2669880/3‑100 (MQ=11)
                                                                                         aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑AGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:1385275/3‑100 (MQ=18)
                                                                                         aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATAGATTGC‑‑ACCATCCTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:1860045/3‑100 (MQ=18)
                                                                                         aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAACCGCGCAAATTCAATATATTGCGGAGAAAAT‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct           >  2:2657290/3‑99 (MQ=21)
                                                                                         aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc            >  2:1024208/3‑100 (MQ=18)
                                                                                         aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc            >  1:337852/3‑100 (MQ=18)
                                                                                         aTGCGGGATGCGGCGTGAACGGCTTATCCGGCCTACAAAGTCGTATAAATTCTATGTATTGC‑‑AGTTTATTG‑‑GTAGGCCTGATGCGCTGCGCTTATCAGGc            <  1:1069460/100‑1 (MQ=21)
                                                                                         aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc            >  1:28644/1‑100 (MQ=255)
                                                                                          tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa           <  1:290373/100‑1 (MQ=255)
                                                                                          tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCGGATAAGCGGAGCGCATCAGGCa           <  1:1687989/100‑1 (MQ=255)
                                                                                          ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGC‑‑‑‑AATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGct           >  2:1983699/2‑99 (MQ=18)
                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGC‑‑‑‑AATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:941263/2‑99 (MQ=18)
                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGC‑‑AGGAGTTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt          >  2:2019648/2‑100 (MQ=21)
                                                                                           gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt          <  2:1799011/100‑1 (MQ=255)
                                                                                           gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGC‑‑AGAAACCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:2772691/1‑98 (MQ=17)
                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑‑AGGACCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  2:1312645/1‑99 (MQ=18)
                                                                                              ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt       <  2:267808/100‑1 (MQ=255)
                                                                                                 tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  2:1107043/99‑1 (MQ=255)
                                                                                                   cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  >  2:884310/1‑100 (MQ=255)
                                                                                                   |                                                                                                        
AGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGC‑‑AGAATTTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705872‑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: