Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,705,985	0	T	C	25.3%	198.1 / 66.7	96	intergenic (‑195/+113)	dppB/dppA	dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (43/28);  new base C (24/0);  total (67/29)
Fisher's exact test for biased strand distribution p-value = 6.09e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.59e-01

TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTC  >  NC_000913/3705886‑3706082                                                                                                    |                                                                                                      tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                       >  2:364784/1‑100 (MQ=255) tgctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCtt                                                                                                       <  2:584583/100‑1 (MQ=255)  gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                      <  2:1464139/100‑1 (MQ=255)  gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                                      >  2:558212/1‑100 (MQ=255)        ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCggn                                                                                                <  1:917350/100‑2 (MQ=255)        ggATGGCGCTGCACGCAGAGACACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                                >  2:1601028/1‑100 (MQ=255)          aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                              >  2:700414/1‑100 (MQ=255)          aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                              <  1:22424/100‑1 (MQ=255)           tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                                             >  2:1081181/1‑100 (MQ=255)                 tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                                      <  1:60431/100‑1 (MQ=255)                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                  >  2:95241/1‑100 (MQ=255)                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                  >  1:255049/1‑100 (MQ=255)                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                  <  1:435353/100‑1 (MQ=255)                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                  >  1:613044/1‑100 (MQ=255)                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                  >  1:223884/1‑100 (MQ=255)                      gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                                 <  1:1564967/100‑1 (MQ=255)                       cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                                >  2:1403461/1‑100 (MQ=255)                       cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                                >  1:500056/1‑100 (MQ=255)                         gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCaaa                                                                              >  2:142290/1‑100 (MQ=255)                          agGCACAGCGCCACCGGGTATACAGATTGTCGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAt                                                                             <  1:214290/100‑1 (MQ=255)                           ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                            >  2:376047/1‑100 (MQ=255)                            gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                           <  1:356241/100‑1 (MQ=255)                            gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATAGTGCAAAtt                                                                            >  2:332043/1‑99 (MQ=255)                              acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCaa                                                                         <  2:934714/100‑1 (MQ=255)                               caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                                        <  2:157182/100‑1 (MQ=255)                                aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATa                                                                       <  1:1403454/100‑1 (MQ=255)                                aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATa                                                                       >  2:672561/1‑100 (MQ=255)                                   gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                                    >  1:1498831/1‑100 (MQ=255)                                    ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                   <  2:609643/100‑1 (MQ=255)                                    ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                   <  2:609744/100‑1 (MQ=255)                                     cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                                  <  2:1505975/100‑1 (MQ=255)                                     cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                                  >  1:605013/1‑100 (MQ=255)                                      aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                                 <  1:228258/100‑1 (MQ=255)                                       ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCa                                                                >  2:1094653/1‑100 (MQ=255)                                       ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCa                                                                <  2:490392/100‑1 (MQ=255)                                         gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑Ga                                                             >  1:300835/1‑100 (MQ=255)                                               aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑Tc                                                     <  2:1647620/100‑1 (MQ=255)                                               aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑Tc                                                     >  2:482641/1‑100 (MQ=255)                                                  gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTa                                                  <  2:1371717/100‑1 (MQ=255)                                                    ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAgg                                                >  2:1250555/1‑100 (MQ=255)                                                      gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGcc                                              <  2:1249185/100‑1 (MQ=255)                                                      gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGcc                                              >  1:1096353/1‑100 (MQ=255)                                                       ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCt                                             >  1:1627216/1‑100 (MQ=255)                                                       ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCt                                             <  2:1031849/100‑1 (MQ=255)                                                        tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTg                                            >  2:152690/1‑100 (MQ=255)                                                        tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTg                                            >  2:198267/1‑100 (MQ=255)                                                         gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATCCATTGCA‑GAAT‑TT‑TCGTAGGCCTGa                                           >  1:1352443/1‑100 (MQ=255)                                                         gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGa                                           >  1:425821/1‑100 (MQ=255)                                                          aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGAt                                          <  1:1317016/100‑1 (MQ=255)                                                             aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAg                                       >  1:795600/1‑100 (MQ=255)                                                                tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCgg                                    >  1:1238273/1‑100 (MQ=255)                                                                 gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGa                                   >  1:554668/1‑100 (MQ=255)                                                                  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAg                                  >  2:193449/1‑100 (MQ=255)                                                                  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAg                                  >  1:865909/1‑100 (MQ=255)                                                                     gatttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑TT‑ATGTAGGCCTGATAAGCGTAGCGCa                              >  2:1115534/10‑100 (MQ=18)                                                                      ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCa                              >  1:423905/1‑100 (MQ=255)                                                                      atttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCCC‑TTGTAGGCCTGATAAGCGTAGCGCAt                             >  2:904568/9‑100 (MQ=18)                                                                       tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCAt                             >  2:1071092/1‑100 (MQ=255)                                                                        gTTTATGCCGGATGCGGCGTTAACGCCTGATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATc                            <  2:71491/100‑1 (MQ=255)                                                                          ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAg                          <  2:298353/100‑1 (MQ=255)                                                                          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑CA‑T‑GTAGGCCTGATAAGCGTAGCGCATCAg                          >  2:1508229/5‑100 (MQ=18)                                                                          gtaggCCGGATAAGGCGTTAACGCCGCATCCGGCATAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAg                          >  2:244991/5‑100 (MQ=18)                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAATCGTGCAAATTCAACATATTGCA‑ATTC‑TC‑TTGTAGGCCTGATAAGCGTAGCGCATCAgg                         >  1:1699646/3‑99 (MQ=18)                                                                            aggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAATCGTGCAAATTCAACATATTGCA‑‑ATT‑CTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc                        >  1:1506162/3‑100 (MQ=18)                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GGAC‑CA‑T‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                       >  1:767661/2‑100 (MQ=18)                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGCA‑‑ATA‑TG‑ACGTAGGCCTGATAAGCGTAGCGCATCAGGcga                      >  2:953868/2‑98 (MQ=21)                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑AT‑TGGTAGGCCTGATAAGCGCAGCGCATCAGGCa                       >  1:1085379/2‑100 (MQ=21)                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑AT‑TGGTAGGCCTGATAAGCGCAGCGCATCAGGCa                       >  2:131046/2‑100 (MQ=21)                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGAT‑CGAT‑CC‑ACGTAGGCCTGATAAGCGTAGCGCATCAGGCa                       >  2:288038/2‑100 (MQ=18)                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑CC‑TCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                       >  2:351287/2‑100 (MQ=21)                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GGAG‑CT‑GCGTAGGCCTGATAAGCGTAGCGCATCAGGCa                       >  2:1338067/2‑100 (MQ=17)                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GAATCCC‑TTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtt                     >  1:1455779/1‑100 (MQ=18)                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCA‑GGGA‑CT‑ATGTAGGCCTGATAAGCATAGCGCATCAGGCag                      >  2:1235653/1‑99 (MQ=18)                                                                               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCAGGGAC‑TA‑T‑GTAGGCCTGATAAGCATAGCGCATCAGGCagc                     >  2:711970/9‑98 (MQ=18)                                                                               ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                     >  2:1528506/1‑100 (MQ=255)                                                                               ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                     >  2:1198727/1‑100 (MQ=255)                                                                                cggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGCA‑ATCT‑TC‑TGGTAGGCCTGATAAGCATCGCGCATCAGGCaagc                    >  2:413237/8‑97 (MQ=18)                                                                                cggataaggcgTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GGAG‑TT‑ATGTAGGCCTGATAAGCGTAGCGCATCAGGCAttta                   >  1:114986/8‑99 (MQ=21)                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GTTT‑AT‑TGGTAGGCCTGATAAGCGCAGCGCATCAGGCaatct                   >  2:681480/7‑96 (MQ=25)                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑CC‑TCGTAGGCCTGATAAGCGTAGCGCATCAGGCaattc                   >  2:36398/7‑96 (MQ=21)                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GGAA‑AC‑ACGTAGGCCTGATAAGCGAAGCGCATCAGGCAgttt                   >  1:1510346/7‑100 (MQ=18)                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GGAA‑AC‑ACGTAGGCCTGATAAGCGAAGCGCATCAGGCAgttt                   >  2:191830/7‑100 (MQ=18)                                                                                  gataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GAAA‑CC‑TTGTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                  >  2:1353041/6‑100 (MQ=18)                                                                                  gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CCAT‑CC‑TCGTAGGCCTGATAAGCGTAGCGCATCAGGCaattct                  >  2:1032851/6‑95 (MQ=25)                                                                                   aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                 <  2:1041634/100‑1 (MQ=255)                                                                                     gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc               <  2:1130083/100‑1 (MQ=255)                                                                                      cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct              >  1:249202/1‑100 (MQ=255)                                                                                      cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct              >  1:390453/1‑100 (MQ=255)                                                                                           ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCACTACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg           <  1:1081175/98‑1 (MQ=255)                                                                                            tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa        >  2:691652/1‑100 (MQ=255)                                                                                             aaCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCAt       >  1:901263/1‑100 (MQ=255)                                                                                              aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc      <  2:255049/100‑1 (MQ=255)                                                                                              aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc      <  2:223884/100‑1 (MQ=255)                                                                                               cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa     >  2:529809/1‑100 (MQ=255)                                                                                               cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa     >  2:1653603/1‑100 (MQ=255)                                                                                                gCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg    <  1:700412/100‑1 (MQ=255)                                                                                                 ccTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGt   >  2:1008388/1‑100 (MQ=255)                                                                                                  cTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc  >  1:523525/1‑100 (MQ=255)                                                                                                    |                                                                                                      TGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GAAT‑TT‑TCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTC  >  NC_000913/3705886‑3706082

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