Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,984 T→G 24.0% 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,9840TG24.0% 587.5 / 198.8 259intergenic (‑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 (89/108);  new base G (59/3);  total (148/111)
Fisher's exact test for biased strand distribution p-value = 9.24e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.20e-02

Only 100 of 263 total aligned reads displayed.

TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083
                                                                                                   |                                                                                                        
ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                                          <  1:750074/100‑1 (MQ=255)
   ctgagggATGGCGCTGCACGCAGAGGCAGCGCGTAAACGGGTGTACAGATTGTTGTCAAAGTGCGATTTGTGTATGGCGTAGCCCGCGTTAACACCttc                                                                                                        <  2:362567/94‑2 (MQ=255)
    tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                                      >  1:1776977/1‑100 (MQ=255)
        ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                  >  1:921884/1‑100 (MQ=255)
         gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                 <  2:3184218/100‑1 (MQ=255)
         gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCGCATCCGGc                                                                                                  <  2:1624355/99‑1 (MQ=255)
           tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                               <  2:3786337/100‑1 (MQ=255)
           tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                               >  1:1736440/1‑100 (MQ=255)
            ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                                              >  2:348834/1‑100 (MQ=255)
             gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                            >  1:1709447/1‑100 (MQ=255)
              cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGGTAACGCCTTATCCGGCCTAC‑aa                                                                                           >  1:105899/1‑100 (MQ=255)
               gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                                          >  1:3664594/1‑100 (MQ=255)
                 tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                                        <  1:3178730/100‑1 (MQ=255)
                  gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                                       <  1:3646792/100‑1 (MQ=255)
                   cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                      <  1:3253302/100‑1 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                    <  1:2944035/100‑1 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGACGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                    >  1:2927163/1‑100 (MQ=255)
                      gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                   <  2:3907159/100‑1 (MQ=255)
                       cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                                  <  1:3704685/100‑1 (MQ=255)
                           ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                              >  1:3424797/1‑100 (MQ=255)
                           agCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                              >  2:1754591/2‑100 (MQ=255)
                               caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                          >  1:4151751/1‑100 (MQ=255)
                               caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                          <  1:144595/100‑1 (MQ=255)
                               caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                          <  2:3087705/100‑1 (MQ=255)
                                aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATa                                                                         <  1:3314644/100‑1 (MQ=255)
                                aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCAGCCTAC‑AAAATCGTGCAAATTCAATa                                                                         <  1:2713647/100‑1 (MQ=255)
                                 gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                                        >  1:2310467/1‑100 (MQ=255)
                                        cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGAATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                                    <  2:2903065/97‑1 (MQ=255)
                                          ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTCC‑AAAATCGTGCAAATTCAATACATTGc                                                                   >  2:2491604/1‑96 (MQ=255)
                                           gTATACAGATTGTTGATAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAt                                                            <  2:2437113/100‑1 (MQ=255)
                                             ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAt‑‑tt                                                        <  1:2488960/100‑1 (MQ=255)
                                              taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGTATTCGGCGTTAACTCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAt‑‑ttt                                                       >  1:2268689/1‑100 (MQ=255)
                                               aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTc                                                      <  1:2098404/100‑1 (MQ=255)
                                                cgatctgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCg                                                     <  2:2288436/95‑1 (MQ=255)
                                                cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCg                                                     <  2:2927146/100‑1 (MQ=255)
                                                cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCg                                                     <  2:2268675/100‑1 (MQ=255)
                                                  gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTa                                                   >  1:2568812/1‑100 (MQ=255)
                                                  gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTa                                                   <  2:4151717/100‑1 (MQ=255)
                                                     tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGc                                                >  2:2981095/1‑100 (MQ=255)
                                                        tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTg                                             <  1:485548/100‑1 (MQ=255)
                                                        tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTg                                             >  2:857927/1‑100 (MQ=255)
                                                        tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTg                                             >  2:110380/1‑100 (MQ=255)
                                                        tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTg                                             >  2:1890898/1‑100 (MQ=255)
                                                        tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTg                                             >  2:2378451/1‑100 (MQ=255)
                                                        cgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTg                                             >  2:281776/2‑100 (MQ=255)
                                                         gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGa                                            >  1:2729033/1‑100 (MQ=255)
                                                          aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGAt                                           >  1:3552480/1‑100 (MQ=255)
                                                           cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATa                                          <  2:2646751/100‑1 (MQ=255)
                                                               gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCg                                      >  1:3743788/1‑100 (MQ=255)
                                                                tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCgg                                     >  2:75980/1‑100 (MQ=255)
                                                                 gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGa                                    >  1:871732/1‑100 (MQ=255)
                                                                  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGa                                    >  1:233793/1‑99 (MQ=255)
                                                                  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAg                                   <  2:3729640/100‑1 (MQ=255)
                                                                    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAgcg                                 <  1:281779/100‑1 (MQ=255)
                                                                     attttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAATCGTGCTAATTCAATATATTGC‑A‑GAAA‑‑CCTTGTAGGCCTGATAAGCGTAgcgc                                >  1:1406379/10‑100 (MQ=17)
                                                                      ttGTTTATGCCGGATGCGCCGTTAACGCCTTATCCGGCCTAC‑AAACTCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCa                               <  1:1157764/100‑1 (MQ=255)
                                                                      atttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑‑GAATC‑CCTTGTAGGCCTGATAAGCGTAGCGCAt                              >  1:2671933/9‑100 (MQ=18)
                                                                       tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCAt                              <  2:1924994/100‑1 (MQ=255)
                                                                        ttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGC‑A‑TGAG‑‑CCACGTAGGCCTGATAAGCGTAGCGCa                               >  1:1294540/7‑98 (MQ=17)
                                                                         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCa                            <  1:3819240/100‑1 (MQ=255)
                                                                         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCa                            <  1:735648/100‑1 (MQ=255)
                                                                         tgtaggCCGGATAAGGCGTTCACTCCGCATCCGGCCTAC‑AAAAGCACGCAAATTCAATATATTGC‑A‑GAGA‑‑TTGCGTAGGCCTGATAAGCGTAGCGCATc                             >  1:3992390/6‑99 (MQ=11)
                                                                          gtaggCCGGATAAGGCGCTTACGCCGCATCCGGCCTAC‑AAGAACGTGCAAATTCAACAAATTGTGATGCAT‑‑CCT‑GTAGGCCTGATAAGCGTAGCGCATCa                            >  1:1023285/5‑100 (MQ=18)
                                                                           taggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAAGCGTGCAAATTCAATACATTGC‑ATGGGC‑‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAgg                          >  1:3892185/4‑100 (MQ=17)
                                                                           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGC‑A‑‑ATA‑‑TGACGTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  2:3794943/4‑100 (MQ=18)
                                                                           tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTCC‑A‑GAAT‑‑TTTCGTAGACCTGATAAGAGGAGCGCATCAgg                          <  2:2567147/100‑1 (MQ=255)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGC‑ATGGGC‑‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:3332077/2‑100 (MQ=17)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTGT‑G‑GGCT‑‑CCACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:1516650/2‑100 (MQ=21)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGC‑A‑GGAG‑‑CTGCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:2057983/2‑100 (MQ=17)
                                                                             ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGc                         >  1:1366761/2‑99 (MQ=18)
                                                                              gtcggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGC‑A‑GGAA‑‑CGATGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  2:1482452/10‑99 (MQ=21)
                                                                              gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt                       >  2:933234/1‑100 (MQ=255)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGC‑A‑GAAA‑‑TCATGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  1:1641078/1‑99 (MQ=18)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGC‑A‑GTTT‑‑ATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaa                       >  1:3919453/1‑99 (MQ=21)
                                                                               ccGGATGCGGAGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGACCTGATAAGCTGAGCGCATAACGCAtt                      >  2:2547629/1‑100 (MQ=18)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑A‑‑ATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCaatt                     >  1:2392307/7‑97 (MQ=21)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGC‑A‑AGAT‑‑TTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCacgct                    >  1:2524187/7‑96 (MQ=17)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGC‑AGGGAC‑‑TAT‑GTAGGCCTGATAAGCATAGCGCATCAGGCagctt                    >  1:4016767/7‑96 (MQ=21)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGC‑A‑GAGA‑‑TTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAattt                    >  1:3613338/7‑100 (MQ=17)
                                                                                 ggacaaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGC‑A‑AGAT‑‑TTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCacgct                    >  2:1197906/7‑96 (MQ=17)
                                                                                  gataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGC‑A‑GGAC‑‑CAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                   >  2:2828621/6‑100 (MQ=18)
                                                                                   aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                  <  2:1699575/100‑1 (MQ=255)
                                                                                    tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                 >  1:4047714/1‑100 (MQ=255)
                                                                                      aggcgTTAACGCCTTTTTTGGCCTAC‑TAAATCATGCAAATTCAATTCATTGC‑A‑GAAT‑‑TTTCGTAGGCCTAATAAGCGGAGCGCATCAGGCATTtttgc                 >  2:3554594/2‑98 (MQ=18)
                                                                                       ngcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATCtttgcctt              >  2:1041620/2‑100 (MQ=255)
                                                                                       ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt              <  2:11452/100‑1 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  2:320447/100‑1 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            <  2:481015/100‑1 (MQ=255)
                                                                                         cgTTAACGCCTTATCCGACCAAC‑‑AAATTATGCAACTCCACAATATTGC‑A‑GAAA‑‑TTAAGTAGGACTTATAAGCATAACGCATAAAACAtttattcgttt             >  2:329578/1‑90 (MQ=255)
                                                                                          gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTTGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc           <  2:883385/100‑1 (MQ=255)
                                                                                          gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc           <  1:3215199/100‑1 (MQ=255)
                                                                                            tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa         >  1:4091024/1‑100 (MQ=255)
                                                                                             aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt        >  1:3365096/1‑100 (MQ=255)
                                                                                              aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       >  1:3396360/1‑100 (MQ=255)
                                                                                               cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAACGCATCAGGCATTTTTGCCTTTGCCATCa      >  2:3967615/1‑100 (MQ=255)
                                                                                                gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCan     <  2:535133/100‑2 (MQ=255)
                                                                                                gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg     >  1:259524/1‑100 (MQ=255)
                                                                                                 ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt    >  1:3514199/1‑100 (MQ=255)
                                                                                                 ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt    >  1:3105986/1‑100 (MQ=255)
                                                                                                 ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt    >  1:3015510/1‑100 (MQ=255)
                                                                                                   |                                                                                                        
TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGC‑A‑GAAT‑‑TTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083

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: ‑