Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,984 T→G 33.1% intergenic (‑194/+114) 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,9840TG33.1% 391.6 / 260.0 218intergenic (‑194/+114)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (78/68);  new base G (69/3);  total (147/71)
Fisher's exact test for biased strand distribution p-value = 9.13e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.99e-03

Only 100 of 224 total aligned reads displayed.

TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705886‑3706083
                                                                                                  |                                                                                                          
 gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                          <  2:322772/100‑1 (MQ=255)
    gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                       >  1:1559613/1‑100 (MQ=255)
     ccGGATGGCGCTGCACGCAGGGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                      >  2:2432760/1‑100 (MQ=255)
     ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCg                                                                                                      <  1:3066418/100‑1 (MQ=255)
       ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                    >  1:2163844/1‑100 (MQ=255)
       ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                    >  1:1685028/1‑100 (MQ=255)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATTCGGcc                                                                                                   <  2:2913279/100‑1 (MQ=255)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                   <  1:2819020/100‑1 (MQ=255)
        gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                                   <  1:2948789/100‑2 (MQ=37)
         aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTACCGCCTTATCCGGCCt                                                                                                  <  2:2775713/100‑1 (MQ=255)
           ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                                                >  2:1675526/1‑100 (MQ=255)
            gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                              >  1:1225803/1‑100 (MQ=255)
            gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATaaa                                                                                               <  1:198795/100‑3 (MQ=25)
             cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                                             >  1:1946221/1‑100 (MQ=255)
              gcTGCACGCAGAAGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                                            >  1:2153862/1‑100 (MQ=255)
               cTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaaa                                                                                           <  2:1032499/100‑1 (MQ=255)
                 gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                                         >  2:1439776/1‑100 (MQ=255)
                 gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATc                                                                                         >  1:1753293/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:1585138/1‑100 (MQ=255)
                        gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                                                                  <  2:548214/100‑1 (MQ=255)
                         agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                                 <  1:773615/100‑1 (MQ=255)
                          ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                                                <  1:109213/100‑1 (MQ=255)
                           gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                               <  1:2685436/100‑1 (MQ=255)
                            cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                              >  1:2234566/1‑100 (MQ=255)
                            cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                              >  2:779388/1‑100 (MQ=255)
                            cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCAGATGCGGCGTTAACGCCTTATACGGCCTACAAAATCGAGCaaa                                                                                  >  2:609378/1‑96 (MQ=255)
                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                            >  1:551611/1‑100 (MQ=255)
                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                            <  2:342601/100‑1 (MQ=255)
                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                            >  2:1117900/1‑100 (MQ=255)
                              caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                            <  2:1983125/100‑1 (MQ=255)
                               aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATa                                                                           >  1:970267/1‑100 (MQ=255)
                                   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                       <  1:192722/100‑1 (MQ=255)
                                   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                       <  2:1892072/100‑1 (MQ=255)
                                        gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGa                                                                 >  1:279554/1‑100 (MQ=255)
                                             taCAGATTGTTGACAAAGTGCGCTCTGTTTATGCCGGATGCGGCGTTCACGCCTTATCCGGCCTACAAAATCGTGCAAAGTCAATACATTGCAGAAtt‑‑‑tt                                                       <  2:750135/100‑1 (MQ=255)
                                              aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTc                                                      <  1:2858016/100‑1 (MQ=255)
                                                  aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAg                                                  <  1:3098765/100‑1 (MQ=255)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGcc                                               >  2:17837/1‑100 (MQ=255)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGcc                                               >  2:3068836/1‑100 (MQ=255)
                                                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGcc                                               >  1:867736/1‑100 (MQ=255)
                                                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCt                                              <  1:481165/99‑1 (MQ=255)
                                                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATa                                          <  2:874867/100‑1 (MQ=255)
                                                            aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAg                                        <  1:1324162/100‑1 (MQ=255)
                                                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGc                                       <  2:2405360/100‑1 (MQ=255)
                                                              gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCg                                      >  1:2363386/1‑100 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAgc                                  <  1:2141095/100‑1 (MQ=255)
                                                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAgc                                  <  2:1924913/100‑1 (MQ=255)
                                                                    aatttgtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAAT‑‑‑CATGTAGGCCTGATAAGCGTAgcgc                                >  1:2477115/10‑100 (MQ=18)
                                                                     ttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑‑TGCGTAGGCCTGATAAGCGTAGCGCa                               >  1:1647782/9‑100 (MQ=17)
                                                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATc                             >  1:811930/1‑100 (MQ=255)
                                                                        tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATT‑‑‑ATGTAGGCCTGATAAGCGTAGCGCATCa                            >  1:2293212/6‑100 (MQ=17)
                                                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGTCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAg                           <  2:1249503/100‑1 (MQ=255)
                                                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTCGTAGGCCTGATAAGCGTAGCGCATCAgg                          >  2:1409460/4‑100 (MQ=21)
                                                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAATTGTGTAAATTCAATAAATTGATCGATC‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAgg                          >  1:2082083/4‑100 (MQ=21)
                                                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAgg                          <  2:3033155/100‑1 (MQ=255)
                                                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAgg                          <  2:1015879/100‑1 (MQ=255)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATAT‑‑‑GACGTAGGCCTGATAAGCGTAGCGCATCAGGcg                        >  2:3097622/3‑99 (MQ=18)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑‑TTGGTAGGCCTGATAAGCGCAGCGCATCAGGc                         >  1:2967283/3‑100 (MQ=25)
                                                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATT‑‑‑ATGTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  1:646642/3‑100 (MQ=17)
                                                                            ggtcggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACG‑‑‑ATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:2614586/11‑100 (MQ=21)
                                                                            ggcaggataaggcgTTCACGCCGCATCCTGCCTACAAAAGCACGCAAATTCAATATATTGCAGAGAT‑‑‑GGCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:1192797/11‑100 (MQ=11)
                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGc                         >  1:2411478/2‑100 (MQ=21)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATC‑‑‑CTCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:2185810/2‑100 (MQ=21)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑‑TGCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:1880016/2‑100 (MQ=18)
                                                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGAT‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:1148090/2‑100 (MQ=17)
                                                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTT‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCag                       >  1:1801059/1‑99 (MQ=17)
                                                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                      >  1:1412405/1‑100 (MQ=21)
                                                                              ccggataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                     >  1:395009/9‑100 (MQ=18)
                                                                              ccggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAGC‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  2:2246570/9‑98 (MQ=17)
                                                                              ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      <  1:592825/100‑1 (MQ=255)
                                                                               cggataaggcgTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGATC‑‑‑CACGTAGGCCTGATAAGCGTAGCGCATCAGGCagtt                     >  1:714282/8‑97 (MQ=21)
                                                                                ggataaggcgTTTACGCCTTATCCGGCCTACAAATTCTTCCAAATTCAATATATTGCAGAAAT‑‑‑CATGTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                    >  1:2925598/7‑100 (MQ=18)
                                                                                ggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAAC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:1294449/7‑100 (MQ=18)
                                                                                ggataaggcgTTTACGCCGCATCCGGCCTACAAATCCTGCAAATTCAACAAATTGCGAATCC‑‑‑CTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgc                   >  1:2999883/7‑98 (MQ=18)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACC‑‑‑ATGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:1083567/7‑100 (MQ=18)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTA‑‑‑TTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                    >  2:525366/7‑96 (MQ=25)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATT‑‑‑TTCGTAGGCCTGATAAGCGTAGCGCATCAGGCacgct                    >  2:460648/7‑96 (MQ=17)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGC‑‑‑TGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:919571/7‑100 (MQ=18)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATT‑‑‑ATGTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                    >  1:93945/7‑100 (MQ=17)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                   >  1:2694792/7‑98 (MQ=21)
                                                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGT‑‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                   >  1:1239947/7‑98 (MQ=21)
                                                                                ggataaggcgTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      >  2:2873490/7‑98 (MQ=18)
                                                                                 gataaggcgTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGAC‑‑‑TATGTAGGCCTGATAAGCATAGCGCATCAGGCagcttt                   >  1:748490/6‑95 (MQ=21)
                                                                                 gataaggcgTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                   >  1:1502061/6‑100 (MQ=18)
                                                                                  aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  2:2692317/100‑1 (MQ=255)
                                                                                   tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                 <  2:1069237/100‑1 (MQ=255)
                                                                                     cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct               >  2:2122041/1‑100 (MQ=255)
                                                                                       gcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgccttt             <  2:2347368/100‑1 (MQ=255)
                                                                                        cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  2:429645/100‑1 (MQ=255)
                                                                                         gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc           >  1:2389130/1‑100 (MQ=255)
                                                                                            aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt        >  1:284058/1‑100 (MQ=255)
                                                                                              cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa      >  1:297365/1‑100 (MQ=255)
                                                                                               gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     >  1:1053770/1‑100 (MQ=255)
                                                                                               gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     <  2:1559609/100‑1 (MQ=255)
                                                                                                  ttATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt  <  2:420945/100‑1 (MQ=255)
                                                                                                  ttATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt  <  1:62525/100‑1 (MQ=255)
                                                                                                  ttATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt  <  1:1339287/100‑1 (MQ=255)
                                                                                                  |                                                                                                          
TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATT‑‑‑TTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705886‑3706083

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: