Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,971 C→A 21.2% 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.2% 474.8 / 116.3 175intergenic (‑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/69);  new base A (33/4);  total (102/73)
Fisher's exact test for biased strand distribution p-value = 9.49e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-03

Only 100 of 189 total aligned reads displayed.

GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705873‑3706070
                                                                                                  |                                                                                                           
gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                            >  1:3261483/1‑100 (MQ=255)
gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                            >  1:2513037/1‑100 (MQ=255)
gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                            >  2:62525/1‑100 (MQ=255)
  gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                          >  1:1017289/1‑100 (MQ=255)
  gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                          <  1:884441/100‑1 (MQ=255)
  gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                          <  1:2659379/100‑1 (MQ=255)
     tGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                                       >  2:234636/1‑100 (MQ=255)
      gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                      >  1:3153233/1‑100 (MQ=255)
       aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                     >  1:1015881/1‑100 (MQ=255)
        aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                    >  1:2171396/1‑100 (MQ=255)
         gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                   >  2:470254/1‑100 (MQ=255)
             tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                               >  1:208334/1‑100 (MQ=255)
              gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                              <  2:322772/100‑1 (MQ=255)
                 gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                           >  1:1559613/1‑100 (MQ=255)
                  ccGGATGGCGCTGCACGCAGGGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                          >  2:2432760/1‑100 (MQ=255)
                    ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                        >  1:1482479/1‑100 (MQ=255)
                    ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                        >  2:868382/1‑100 (MQ=255)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATTCGGcc                                                                                       <  2:2913279/100‑1 (MQ=255)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                       <  2:1031116/100‑1 (MQ=255)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                       <  1:2948789/100‑2 (MQ=37)
                        ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                                    <  1:2160077/100‑1 (MQ=255)
                         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                  >  1:2937429/1‑100 (MQ=255)
                         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                  >  2:1012478/1‑100 (MQ=255)
                         gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATaaa                                                                                   <  1:621100/100‑3 (MQ=25)
                          cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                                 >  1:1946221/1‑100 (MQ=255)
                          cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                                 >  1:1886916/1‑100 (MQ=255)
                               cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                            >  1:288168/1‑100 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                          >  1:1545066/1‑100 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                          >  2:660961/1‑100 (MQ=255)
                                   cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                        >  1:1599740/1‑100 (MQ=255)
                                   cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                        <  2:2685832/100‑1 (MQ=255)
                                     gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                                                      <  2:548214/100‑1 (MQ=255)
                                       ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                                    <  1:109213/100‑1 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                  <  2:1923895/100‑1 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                  >  1:420948/1‑100 (MQ=255)
                                           caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                <  2:1983125/100‑1 (MQ=255)
                                           caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                <  2:1272215/100‑1 (MQ=255)
                                           caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                <  2:1048938/100‑1 (MQ=255)
                                             gcgcCACCGAGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                              <  2:2949647/100‑1 (MQ=255)
                                               gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAt                                                            >  1:3075504/1‑100 (MQ=255)
                                                ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                           <  1:192722/100‑1 (MQ=255)
                                                     gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGa                                                     >  1:279554/1‑100 (MQ=255)
                                                           aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAtt‑‑‑tt                                            <  2:523417/99‑1 (MQ=255)
                                                           aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTc                                           >  2:933832/1‑100 (MQ=255)
                                                           aCAGATTGTTGACAAAGTGCGCTTTGTGTCTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTc                                           <  2:3076525/100‑1 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑g                                        <  2:3261465/100‑1 (MQ=255)
                                                                  gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGcc                                  >  2:17837/1‑100 (MQ=255)
                                                                    tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCt                                 <  1:481165/99‑1 (MQ=255)
                                                                    tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTg                                >  1:3003748/1‑100 (MQ=255)
                                                                     gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGa                               <  1:1199259/100‑1 (MQ=255)
                                                                     gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGa                               >  2:118135/1‑100 (MQ=255)
                                                                     gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGa                               >  1:2692331/1‑100 (MQ=255)
                                                                      aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGAt                              <  2:2910192/100‑1 (MQ=255)
                                                                      aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGAt                              <  2:1651688/100‑1 (MQ=255)
                                                                         aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAg                           <  1:1324162/100‑1 (MQ=255)
                                                                          tttGCGTTTTGTTTATGTCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGc                          <  1:83484/98‑1 (MQ=255)
                                                                          aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGc                          <  2:2405360/100‑1 (MQ=255)
                                                                           gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCg                         >  1:2363386/1‑100 (MQ=255)
                                                                            tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                        >  2:3191119/1‑100 (MQ=255)
                                                                             gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGa                       >  1:1653061/1‑100 (MQ=255)
                                                                              cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                      >  2:2161221/1‑100 (MQ=255)
                                                                               gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                     <  1:2141095/100‑1 (MQ=255)
                                                                               gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                     <  1:1384898/100‑1 (MQ=255)
                                                                               gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                     >  2:738201/1‑100 (MQ=255)
                                                                               gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                     <  2:752072/100‑1 (MQ=255)
                                                                                 aatttgtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAAT‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAgcgc                   >  1:2477115/10‑100 (MQ=18)
                                                                                  ttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCa                  >  1:1647782/9‑100 (MQ=17)
                                                                                    gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCAtn                <  2:450400/100‑2 (MQ=255)
                                                                                    gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                <  1:2389975/100‑1 (MQ=255)
                                                                                     tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATT‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCa               >  1:2293212/6‑100 (MQ=17)
                                                                                      ttATGCCGGATGCGGCGTTAACGCCTTATCCGGTCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg              <  2:1249503/100‑1 (MQ=255)
                                                                                      ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg              <  1:93256/100‑1 (MQ=255)
                                                                                      gtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAAT‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATc                >  2:1441463/5‑98 (MQ=17)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAg              >  2:2271386/5‑100 (MQ=18)
                                                                                       taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  2:169374/4‑100 (MQ=17)
                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  2:1429717/4‑100 (MQ=18)
                                                                                       taggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAATTGTGTAAATTCAATAAATTGATCGATC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  1:2082083/4‑100 (MQ=21)
                                                                                       tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg             <  2:3033155/100‑1 (MQ=255)
                                                                                       tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg             <  1:585578/100‑1 (MQ=255)
                                                                                       tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg             <  1:285737/100‑1 (MQ=255)
                                                                                        aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGCAATCTT‑‑‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGc            >  2:1102279/3‑100 (MQ=18)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATT‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:646642/3‑100 (MQ=17)
                                                                                         tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa           <  1:219719/100‑1 (MQ=255)
                                                                                         ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:2725752/2‑100 (MQ=17)
                                                                                         ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:1239303/2‑100 (MQ=17)
                                                                                         ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:1869722/2‑100 (MQ=17)
                                                                                         ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc            >  1:2411478/2‑100 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:940196/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:371943/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATAT‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcga          >  1:2970106/2‑98 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATAT‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcga          >  1:673225/2‑98 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:2185810/2‑100 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGAT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:1148090/2‑100 (MQ=17)
                                                                                          gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt          <  1:500415/100‑1 (MQ=255)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt         >  2:1989500/1‑100 (MQ=18)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  1:1801059/1‑99 (MQ=17)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑AT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  1:81050/1‑99 (MQ=18)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATAT‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcgat         >  1:172112/1‑97 (MQ=21)
                                                                                           ccGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCT‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGctga         >  1:2803652/1‑97 (MQ=18)
                                                                                                tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc    <  2:1069237/100‑1 (MQ=255)
                                                                                                  |                                                                                                           
GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705873‑3706070

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