Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,984 T→G 30.2% 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,9840TG30.2% 317.5 / 182.5 179intergenic (‑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 (58/67);  new base G (47/7);  total (105/74)
Fisher's exact test for biased strand distribution p-value = 1.81e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-02

Only 100 of 179 total aligned reads displayed.

GCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCA  >  NC_000913/3705890‑3706079
                                                                                              |                                                                                                   
   ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcagat                                                                                             >  1:321559/1‑96 (MQ=255)
   ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                             >  1:660360/1‑100 (MQ=255)
   ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                             <  2:1714430/100‑1 (MQ=255)
    gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                            <  1:2414257/100‑1 (MQ=255)
    gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                            <  1:2031458/100‑1 (MQ=255)
    gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                            <  1:1642073/100‑2 (MQ=37)
     aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                           <  1:1561210/100‑3 (MQ=37)
        gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTata                                                                                        <  1:1866528/100‑3 (MQ=255)
        gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                                       >  1:1003293/1‑100 (MQ=255)
          gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaaa                                                                                     >  2:2176277/1‑100 (MQ=255)
            tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAAt                                                                                   <  2:6324/100‑1 (MQ=255)
                cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                               >  2:2527939/1‑100 (MQ=255)
                 gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                              <  2:250606/100‑1 (MQ=255)
                  cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaa                                                                            >  1:1979869/1‑100 (MQ=255)
                     agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                          <  2:2110829/100‑1 (MQ=255)
                       gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                        <  2:1756539/100‑1 (MQ=255)
                       gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                        <  1:327576/100‑1 (MQ=255)
                       gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                                        <  1:2665575/100‑1 (MQ=255)
                        cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                       >  2:666547/1‑100 (MQ=255)
                        cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                                       >  1:1308065/1‑100 (MQ=255)
                        cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGACTTATCCGGCCTACAAAATCGTGCAACTTc                                                                        >  2:200813/1‑99 (MQ=255)
                            gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACTCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                   <  1:1941499/100‑1 (MQ=255)
                              gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAt                                                                 >  1:354841/1‑100 (MQ=255)
                                cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTg                                                               <  1:1694706/100‑1 (MQ=255)
                                cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTg                                                               >  2:2631388/1‑100 (MQ=255)
                                 aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                              <  1:1133869/100‑1 (MQ=255)
                                   cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAg                                                           <  2:768951/100‑1 (MQ=255)
                                         taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAtttt                                                   >  2:206616/1‑100 (MQ=255)
                                         taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAtttt                                                   >  2:241697/1‑100 (MQ=255)
                                           cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCg                                                 <  1:729138/100‑1 (MQ=255)
                                           cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCg                                                 <  1:1605600/100‑1 (MQ=255)
                                             gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTa                                               >  2:46913/1‑100 (MQ=255)
                                               ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAgg                                             <  2:1345977/100‑1 (MQ=255)
                                               gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCGGGCCTCCAAAGTCGTGCAAGTTCAATACATTGTAGAAgttt                                                   >  2:2412171/7‑94 (MQ=14)
                                                tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGc                                            <  1:699052/100‑1 (MQ=255)
                                                 gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGcc                                           <  2:710041/100‑1 (MQ=255)
                                                 gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGcc                                           <  2:2429716/100‑1 (MQ=255)
                                                 gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGc                                            >  1:1061909/1‑99 (MQ=255)
                                                  ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCt                                          <  1:554829/100‑1 (MQ=255)
                                                   tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTg                                         <  1:300445/100‑1 (MQ=255)
                                                   tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTg                                         >  1:1220658/1‑100 (MQ=255)
                                                    gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGa                                        >  1:407172/1‑100 (MQ=255)
                                                    gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGa                                        >  2:942258/1‑100 (MQ=255)
                                                     aCAAAGTGTTCTTTGTTGATGCCGGAGGCGGCGTTAACGCCTTATCCGGCCTACATAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGAt                                       <  2:606249/100‑1 (MQ=18)
                                                     aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGAt                                       <  1:2648097/100‑1 (MQ=255)
                                                         aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGc                                   <  1:2450090/100‑1 (MQ=255)
                                                          gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCg                                  <  1:2591667/100‑1 (MQ=255)
                                                           tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCgg                                 <  2:1711901/100‑1 (MQ=255)
                                                           tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCgg                                 <  1:501753/100‑1 (MQ=255)
                                                             cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAg                               >  2:2447668/1‑100 (MQ=255)
                                                               cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATaagcgaagcg                             >  2:308633/1‑100 (MQ=255)
                                                                 ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAgtg                             <  1:1842417/98‑3 (MQ=255)
                                                                   gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATc                         >  1:247865/1‑100 (MQ=255)
                                                                      taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTCACGTAGGCCTGATAAGCGTAGCGCATCAgg                      >  1:534579/4‑100 (MQ=17)
                                                                      tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAgg                      >  1:2365747/1‑100 (MQ=255)
                                                                       aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                    >  1:1392348/3‑100 (MQ=18)
                                                                       aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGc                     >  2:101518/3‑100 (MQ=21)
                                                                       aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGTTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGc                     >  2:6433/3‑100 (MQ=18)
                                                                       aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGTTATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                    >  1:2052756/3‑100 (MQ=21)
                                                                       aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGc                     >  1:1507151/1‑100 (MQ=255)
                                                                       aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGc                     >  1:1507962/1‑100 (MQ=255)
                                                                        ggctggataaggcgTTCACGCCGCATCCGGCCTACAAAATTGTGCAAAATCAATGGATTGCACAGCTAGCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                    >  2:840190/11‑100 (MQ=21)
                                                                        ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                   >  1:483198/2‑100 (MQ=18)
                                                                        ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc                     >  1:287713/2‑99 (MQ=17)
                                                                        ggCCGGATAAGGCGTTCACTCCGCATCCGGCCTACAAAAGCACGCAAATTCACTATATTGCAGAGATTGCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                    >  2:1641264/2‑100 (MQ=11)
                                                                        ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCATGTAGGCCTGATAAGCGTAGCGCATCAGGCa                    >  1:2158119/2‑100 (MQ=18)
                                                                        ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCa                    >  1:2152905/2‑100 (MQ=21)
                                                                        ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCa                    >  1:2281063/2‑100 (MQ=18)
                                                                        agCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                   >  1:225696/2‑100 (MQ=18)
                                                                         gtcggataaggcgTTCACGCCGCATCCGACCCACAAAGTCTTGCAAAATCAATAAATTGCAGGAACGATGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                   >  1:1538163/10‑99 (MQ=21)
                                                                         gCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt                   >  1:594792/1‑100 (MQ=255)
                                                                         gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCag                   >  1:1011896/1‑99 (MQ=18)
                                                                         gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTCACGTAGGCCTGTTAAGCGTAGCGCATCAGGCag                   >  2:537208/1‑99 (MQ=17)
                                                                          ccggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGATCATGTAGGCCTGATAAGCGTAGCGCATAAGGCaat                  >  2:2315923/9‑98 (MQ=14)
                                                                          ccggataaggcgTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaat                  >  2:87525/9‑98 (MQ=18)
                                                                          ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                  >  2:551645/1‑100 (MQ=255)
                                                                          ccGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCagc                  >  2:2251430/1‑98 (MQ=18)
                                                                           cggataaggcgTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                 >  2:1117463/8‑97 (MQ=18)
                                                                           cggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGTTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAttta                >  1:192492/8‑99 (MQ=21)
                                                                            ggataaggcgTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                >  1:1768415/7‑100 (MQ=18)
                                                                            ggataaggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTCACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                >  1:85771/7‑100 (MQ=17)
                                                                            ggataaggcgTTTACGCCGCATCCGGCCTACAAAACAGTGCAAATTCAATATATTGCAATCTTCTGGTAGGCCTGATAAGCATCGCGCATCAGGCaagct                >  2:2446629/7‑96 (MQ=21)
                                                                            ggataaggcgTTCACGCCGCATCCGGCCTACAAGATCGTGCAAATTCAATATATTGCATGAGCTACGTAGGCCTGATAAGCGTAGCGCATCCGGCAgttt                >  2:1073053/7‑100 (MQ=17)
                                                                            ggataaggcgTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                >  2:2731441/7‑96 (MQ=25)
                                                                            ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCATGAAAATTCAATATATTGCAGGAGTTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                >  1:2690385/7‑100 (MQ=17)
                                                                            ggataaggcgTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGATCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                >  1:1709520/7‑100 (MQ=17)
                                                                              aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg              <  2:144561/100‑1 (MQ=255)
                                                                               taaggcgTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGATCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAatttagc             >  2:2112853/4‑97 (MQ=17)
                                                                               tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt               >  2:1544002/1‑98 (MQ=255)
                                                                               tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc             <  1:1879496/100‑1 (MQ=255)
                                                                                gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc            <  1:1086667/100‑1 (MQ=255)
                                                                                  ggcgTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTCACGTAGGCCTGATAAGCGTAGCGCATCAGGCAGTtttgcgtt          <  1:1107760/100‑1 (MQ=17)
                                                                                  ggcgTTAACGCCTTATCCGGCCTACAAAATCGTACAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt          >  2:1162677/1‑100 (MQ=255)
                                                                                    cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg        >  1:312692/1‑100 (MQ=255)
                                                                                     gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc       <  2:1507955/100‑1 (MQ=255)
                                                                                      ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc      <  1:2527957/100‑1 (MQ=255)
                                                                                      ttAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc      <  2:2221099/100‑1 (MQ=255)
                                                                                       tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa     >  2:1454594/1‑100 (MQ=255)
                                                                                       tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa     <  2:748485/100‑1 (MQ=255)
                                                                                         aCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc   <  2:1199266/100‑1 (MQ=255)
                                                                                              |                                                                                                   
GCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCA  >  NC_000913/3705890‑3706079

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: