Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,985 T→C 30.4% intergenic (‑195/+113) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9850TC30.4% 478.2 / 249.3 299intergenic (‑195/+113)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (91/117);  new base C (82/9);  total (173/126)
Fisher's exact test for biased strand distribution p-value = 3.84e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.67e-03

Only 100 of 303 total aligned reads displayed.

GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705887‑3706084
                                                                                                  |                                                                                                       
gctgcCGGATGGCGTTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                        <  2:3423004/100‑1 (MQ=255)
 ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                       <  2:1636553/100‑1 (MQ=255)
  tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                      >  1:2368879/1‑100 (MQ=255)
    ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCCTATcc                                                                                                     <  2:3943944/99‑1 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCGTTATCCgggc                                                                                                 >  1:3479435/1‑98 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                 <  2:4061381/100‑1 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                 <  1:958875/100‑1 (MQ=255)
       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                 <  2:2935128/100‑1 (MQ=255)
        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                                <  1:3624837/100‑3 (MQ=37)
         tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                               <  1:3333944/100‑1 (MQ=255)
         tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATa                                                                                               <  2:307625/100‑1 (MQ=37)
          ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                                              <  1:1576143/100‑1 (MQ=255)
           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                             >  1:2450178/1‑100 (MQ=255)
           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                             >  2:635904/1‑100 (MQ=255)
            cgcTGCACGCAGAGGCACAGCGCCACCGGGTATTCAGATTGTGGACAAAGTGCGCGATGTGTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATaaac                                                                                            <  2:1334577/100‑4 (MQ=18)
             gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCTCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                                           <  1:3772266/100‑1 (MQ=255)
             gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                                           >  1:2827324/1‑100 (MQ=255)
                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGt                                                                                      >  1:389396/1‑100 (MQ=255)
                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGCTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                   >  2:451639/1‑100 (MQ=255)
                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                                   >  1:293470/1‑100 (MQ=255)
                        agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                                <  2:924853/100‑1 (MQ=255)
                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                              <  1:372316/100‑1 (MQ=255)
                          gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                              <  1:3825801/100‑1 (MQ=255)
                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                             >  1:769454/1‑100 (MQ=255)
                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                                           >  1:2672826/1‑100 (MQ=255)
                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATa                                                                          <  1:464300/100‑1 (MQ=255)
                                  ccACCGGGTCTACACATTGTTGCCAACGTGCGCTTTGTTTATGCCAGATGCGGCGTTAACGCCTTATCCGGCCTACAAACTCGTCCAAATTCAATACAtt                                                                      >  2:3822937/1‑100 (MQ=17)
                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                                      <  2:3299980/100‑1 (MQ=255)
                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                                    >  1:2314389/1‑100 (MQ=255)
                                            taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAt‑‑ttt                                                        <  1:3043425/100‑1 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCg                                                      <  1:2779013/100‑1 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCg                                                      >  1:2535200/1‑100 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCg                                                      <  1:4030972/100‑1 (MQ=255)
                                              cAGATTGTTGACAAAGTGCGCTTTGTGTATGCCGGATGCTGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTc                                                       <  2:1720530/99‑1 (MQ=255)
                                               aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGt                                                     >  1:2651557/1‑100 (MQ=255)
                                               aGATTGGTGACAAAGCGCGCGTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGt                                                     <  1:1371869/100‑1 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTa                                                    <  1:222236/100‑1 (MQ=255)
                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTa                                                    <  2:699352/100‑1 (MQ=255)
                                                  ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAgg                                                  <  2:1489793/100‑1 (MQ=255)
                                                   tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGc                                                 <  2:2627466/100‑1 (MQ=255)
                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGcc                                                >  1:800368/1‑100 (MQ=255)
                                                     ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCAGGCCTACAAAATCGTGCAAATTCAAGACATTGAAGAAT‑‑TGTCGtgtg                                                  >  2:3701934/1‑94 (MQ=14)
                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGcc                                                >  2:2669298/1‑98 (MQ=255)
                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTg                                              <  1:1206307/100‑1 (MQ=255)
                                                        aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGAt                                            <  1:1830529/100‑1 (MQ=255)
                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATa                                           <  2:1919889/100‑1 (MQ=255)
                                                           acgTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCATAGTCCTGATaa                                          >  2:2796824/3‑99 (MQ=255)
                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGc                                        >  1:3009585/1‑100 (MQ=255)
                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCg                                       <  2:955020/100‑1 (MQ=255)
                                                               gcgcTTTGTTTATGCCGGATTCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGa                                     >  2:2524683/1‑100 (MQ=255)
                                                               gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGa                                     <  1:542938/100‑1 (MQ=255)
                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAgc                                   >  1:3551034/1‑100 (MQ=255)
                                                                   attttgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAgcgc                                 >  2:2437288/10‑100 (MQ=17)
                                                                    ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCa                                <  1:635910/100‑1 (MQ=255)
                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATc                              <  2:1163147/100‑1 (MQ=255)
                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATc                              >  2:768741/1‑100 (MQ=255)
                                                                       tttATGCCGGATGCGTCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCa                             <  1:4070090/100‑1 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCa                             >  2:2060689/1‑100 (MQ=255)
                                                                       tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCa                             <  2:3555612/100‑1 (MQ=255)
                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAg                            <  1:1457480/100‑1 (MQ=255)
                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAGTCGTATAAATTCAATGTATTGCAGTTT‑‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAg                            >  1:2081561/5‑100 (MQ=21)
                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGAT‑‑CATGTAGGCCTGATAAGCGTAGCGCATCAgg                           >  2:585792/4‑100 (MQ=17)
                                                                         tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAgg                           <  2:2747110/100‑1 (MQ=255)
                                                                         tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAgg                           <  2:215784/100‑1 (MQ=255)
                                                                         tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAgg                           >  2:305095/1‑100 (MQ=255)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAA‑‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:3602104/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGT‑‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  2:1803595/3‑100 (MQ=14)
                                                                          aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGT‑‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:1492222/3‑100 (MQ=17)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGA‑‑CTATGTAGGCCTGATAAGCATAGCGCATCAGGc                          >  1:3151707/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGAT‑‑ACGTGTAGGCCTGATAAGCGTAGCGCATCAGGc                          >  1:3217352/3‑100 (MQ=18)
                                                                          aggCCGGATAAGGCATTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGAC‑‑TATGTAGGCCTGATAAGCATAGCGCATCAGGc                          >  1:1875500/3‑100 (MQ=18)
                                                                           ggtcggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAAC‑‑GATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:792806/11‑100 (MQ=21)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                        >  1:1985800/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAA‑‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:2572890/2‑100 (MQ=17)
                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGc                          >  1:1037291/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCAT‑‑CCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:3176860/2‑100 (MQ=21)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAG‑‑CTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:1112890/2‑100 (MQ=18)
                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGAT‑‑TATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  2:2372723/2‑100 (MQ=17)
                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGT‑‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                         >  1:976259/1‑99 (MQ=14)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTT‑‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaa                        >  2:3351670/1‑99 (MQ=21)
                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAG‑‑CTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCag                        >  1:1226572/1‑99 (MQ=17)
                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                       >  2:1140502/1‑100 (MQ=255)
                                                                              tggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTT‑‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaatc                      <  1:3426685/93‑4 (MQ=25)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAA‑‑CCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  2:2875162/7‑100 (MQ=18)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAA‑‑CCATGTAGGCCTGATCAGCGTAGCGCATCAGGCAgtt                      >  2:82242/7‑99 (MQ=17)
                                                                               ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGT‑‑TCACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:2021854/7‑100 (MQ=17)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAATTGTGTCAATTCAATAAATTGATCGAT‑‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                       >  2:3641914/7‑96 (MQ=18)
                                                                               ggataaggcgTTCACGCCGCATCCGGCCTACAAAATTGTGTAAATTCAATAAATTGATCGAT‑‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                     >  1:667742/7‑100 (MQ=21)
                                                                               ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt                     >  1:878592/1‑100 (MQ=255)
                                                                                gataaggcgTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAA‑‑ACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAgtttt                    >  1:1284444/6‑100 (MQ=18)
                                                                                gataaggcgTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAA‑‑ACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAgtttt                    >  1:3978452/6‑100 (MQ=18)
                                                                                    cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                >  1:39669/1‑100 (MQ=255)
                                                                                    cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                >  1:2987958/1‑100 (MQ=255)
                                                                                    cggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct                >  1:3401507/1‑100 (MQ=255)
                                                                                     ggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt               >  2:1049565/1‑100 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg             <  1:2602455/100‑1 (MQ=255)
                                                                                       cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg             >  1:1647105/1‑100 (MQ=255)
                                                                                           aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt         <  1:2414212/100‑1 (MQ=255)
                                                                                              gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATctg      >  2:2496787/1‑98 (MQ=255)
                                                                                                  tATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCtt  >  1:2278793/1‑100 (MQ=255)
                                                                                                  |                                                                                                       
GCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTT  >  NC_000913/3705887‑3706084

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