Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,232,484	A→G	18.4%	intergenic (+30/+16)	umuC → / ← dsbB	translesion error‑prone DNA polymerase V subunit; DNA polymerase activity/oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,232,484	0	A	G	18.4%	197.2 / 41.9	98	intergenic (+30/+16)	umuC/dsbB	translesion error‑prone DNA polymerase V subunit; DNA polymerase activity/oxidoreductase that catalyzes reoxidation of DsbA protein disulfide isomerase I
Reads supporting (aligned to +/- strand):  ref base A (52/28);  new base G (15/3);  total (67/31)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-01

GCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAA  >  NC_000913/1232385‑1232582                                                                                                    |                                                                                                   gCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCa                                                                                                    <  1:203815/100‑1 (MQ=255)  cAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGAAACAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAg                                                                                                   >  1:269454/1‑100 (MQ=255)    gATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCatg                                                                                                  <  1:2051037/99‑3 (MQ=255)     aTGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAg                                                                                                <  1:272214/100‑1 (MQ=255)        aaGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAgcg                                                                                              <  1:1046008/99‑1 (MQ=255)         aGCGAGCCATGCTTTCACCACGTTATACAACGCGAATTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGcc                                                                                            <  1:856024/100‑1 (MQ=255)          gCGAGCCATGCTTTCACCACGTTATACGACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCg                                                                                           <  1:1718419/100‑1 (MQ=255)              gCCATGCTTTCACCACGTTATACAACGCGACGTTCTGATTTACTGAGGGTCAAATAAATATAACGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATg                                                                                       >  1:990641/1‑100 (MQ=255)               ccATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCTGCGGGAgcgcttttttc                                                                                      <  1:1215345/100‑8 (MQ=255)               ccATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATgg                                                                                      <  2:589861/100‑1 (MQ=255)                cATGCTTTCACCACGTTATACAACGAGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAACAGAGCTCCCGCAAGAGCGCCCAATgg                                                                                      >  2:702933/1‑99 (MQ=255)                 aTGCTTTCACGACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGCAATAAAGCGCTCCTGCGGGAgcgcttttttcct                                                                                    <  2:1995530/100‑10 (MQ=37)                  ttcTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGAtt                                                                                   <  2:2100407/98‑1 (MQ=255)                  ttcTTTCACCAAGTTATACAACGCGAAGTTTTGTTTTACTTTGGGTCAAATAAATATAGAGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGAt                                                                                    <  1:912294/97‑1 (MQ=255)                   gCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTa                                                                                  <  2:229270/100‑1 (MQ=255)                      ttCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCg                                                                               <  1:119637/100‑1 (MQ=255)                       tCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGa                                                                              >  1:1874660/1‑100 (MQ=255)                        caccacGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGAAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGa                                                                              >  2:562998/1‑99 (MQ=255)                         accacGTTATACAACGCGAAGGTCTGATTTAGTGGGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCTCGCAGGAGCGACGAATGGATTAGCGAc                                                                             >  2:1349269/1‑99 (MQ=255)                             ggTTATACAACGCGAAGTTCTGATATACTGAGGGAGAAATAAATATAGCGGCAGGAAAAAAGGGCTCCAGCAGGAGAGACGAATGGATTAGCGACCGa                                                                          >  2:2201374/2‑98 (MQ=255)                             cGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAAc                                                                        <  2:621037/100‑1 (MQ=255)                              gTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACa                                                                       >  2:2669594/1‑100 (MQ=255)                              gTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACa                                                                       <  2:1006929/100‑1 (MQ=255)                                      aCGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATATCGCCAGGGAAAAATTGCTCCCGGAGGAGCTGCGAATGGATTAGCGACCGAACAGATCACGt                                                               >  2:1649692/1‑100 (MQ=255)                                      aCGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGt                                                               >  2:1427077/1‑100 (MQ=255)                                         cgAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGtttt                                                            <  1:2720734/100‑1 (MQ=255)                                            aGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTc                                                         >  2:676215/1‑100 (MQ=255)                                            aGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTc                                                         >  2:1554156/1‑100 (MQ=255)                                            aGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGATCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTc                                                         >  2:675598/1‑100 (MQ=255)                                             gTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCg                                                        <  1:1461224/100‑1 (MQ=255)                                             gTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGATTTTTCg                                                        >  2:1901967/1‑100 (MQ=255)                                              ttCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGTCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGc                                                       <  2:1874638/100‑1 (MQ=255)                                              ttCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGc                                                       >  1:2625413/1‑100 (MQ=255)                                              ttCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGc                                                       >  1:398033/1‑100 (MQ=255)                                                cTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCtt                                                     <  2:28666/100‑1 (MQ=255)                                                cTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATAGAATAGCGACCGAACAGATCACGTTTTTTCGCtt                                                     >  1:2832201/1‑100 (MQ=255)                                                  gATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAACAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTa                                                   >  2:1620871/1‑100 (MQ=255)                                                  gATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAATCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACAGAACAGATCACGTTTTTTCGCTTTa                                                   >  1:1151471/1‑100 (MQ=255)                                                      tACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACg                                               >  2:487427/1‑100 (MQ=255)                                                        cTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATGACGATTTTTCGCTTTAAACGGc                                             >  2:1195324/1‑100 (MQ=255)                                                         tgcgggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACCGAACACGTTTTTTCGCttt                                                    >  2:2283947/4‑92 (MQ=255)                                                          gAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCt                                            >  2:1879649/1‑99 (MQ=255)                                                          gAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTg                                           <  2:9772/100‑1 (MQ=255)                                                          gAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCACAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCt                                            >  2:1346651/1‑99 (MQ=255)                                                           aGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTgg                                          <  1:487430/100‑1 (MQ=255)                                                             ggTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGAGCGAACAAATAACGTTTTTTCGCTTTAAATGGCttgaa                                        >  1:1206981/1‑96 (MQ=255)                                                              gTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTTGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGaa                                       >  2:1935174/1‑100 (MQ=255)                                                               tCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGaaa                                      <  2:347163/100‑1 (MQ=255)                                                               tCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGaaa                                      <  2:100798/100‑1 (MQ=255)                                                                cAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAAt                                     >  1:645458/1‑100 (MQ=255)                                                                 aaataaatATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATc                                    <  1:1411478/100‑1 (MQ=255)                                                                 aaataaatATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATc                                    <  1:3022296/100‑1 (MQ=255)                                                                  aataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATca                                   >  1:2125128/1‑100 (MQ=255)                                                                   ataaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATcac                                  <  1:2918409/100‑1 (MQ=255)                                                                    taaataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATcacc                                 >  1:1124229/1‑100 (MQ=255)                                                                      aataTAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATcaccac                               <  2:2544092/100‑1 (MQ=255)                                                                       atatAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCaccacc                              >  1:424077/1‑100 (MQ=255)                                                                       atatAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGAATTTTCGCTTTAAACGGCTGGGAAATcaccac                               >  2:1439108/1‑99 (MQ=255)                                                                        tataGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCAccacca                             >  2:2862329/1‑100 (MQ=255)                                                                            gCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCAc                         >  1:2064800/1‑100 (MQ=255)                                                                              ggCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTg                       >  1:733872/1‑100 (MQ=255)                                                                              ggCAGGAAAAAAGCGCTCCCACAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTGTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTg                       <  1:643915/100‑1 (MQ=255)                                                                               gCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGc                      >  1:1704111/1‑100 (MQ=255)                                                                                cAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCg                     >  2:1264325/1‑100 (MQ=255)                                                                                  tccattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCTCTTTATACTGCTGGGAAATCACCACCAACACTGCGAc                   >  2:1333540/9‑100 (MQ=37)                                                                                  ggAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGAc                   >  2:1114100/1‑100 (MQ=255)                                                                                  ggAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGAc                   >  1:1206522/1‑100 (MQ=255)                                                                                  ggAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGAc                   <  1:2266184/100‑1 (MQ=255)                                                                                   gAAAAAAGCGCTCCCGCAGGAGCGCCGACTGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACa                  >  1:445889/1‑100 (MQ=255)                                                                                   gAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACa                  >  1:1357706/1‑100 (MQ=255)                                                                                    gatctggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaa                 <  1:1857973/94‑1 (MQ=255)                                                                                    cattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaa                 >  1:1326521/7‑100 (MQ=255)                                                                                    cattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaa                 >  1:2074289/7‑100 (MQ=255)                                                                                    cattcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGACTGGGAAATAACCACCAGAACTGCGACaa                 >  1:2603269/7‑100 (MQ=255)                                                                                    aaaaaaGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaa                 >  2:37123/1‑100 (MQ=255)                                                                                     attcggcgcTCCTGCGGGAGCGCTGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt                >  1:2744964/6‑100 (MQ=255)                                                                                     attcggcgcTCCTGCGGGAGCGCCGAATTGATTAGCGACCGAACAGATCACGTTTTTTCGCAATAAACAACAAGAAAAACACCACCAGCACTGCGACaa                 >  1:1058340/6‑99 (MQ=255)                                                                                     attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAATTGCTGGGAAATCACCACCAGCACTGCGaga                  >  1:661885/6‑96 (MQ=255)                                                                                     attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt                >  1:2439125/6‑100 (MQ=255)                                                                                     attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAAt                >  1:2334292/6‑100 (MQ=255)                                                                                     attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAAACACCACAAGCACTGCGACAAt                >  1:1741839/6‑100 (MQ=255)                                                                                     attcggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCCCTTTAAACGGCTGGGAAATCACCACCAACACTGCGACaa                 >  1:2229220/6‑99 (MQ=255)                                                                                        cggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGTTTTAAACTGCTGGGAAATCACAACCAGCACTGCGACAATCa              >  2:2361605/3‑99 (MQ=255)                                                                                        aaGCGCTCCCGCAGGAGCGCTGAATGGATTAGTGTCCGAACAGGTCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACCGCGACAATc               <  1:1194902/98‑1 (MQ=255)                                                                                         ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACaaagatc            >  2:1857952/2‑95 (MQ=255)                                                                                         ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg            >  1:946715/2‑100 (MQ=255)                                                                                         ggcgcTCCTGCGGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg            >  2:1984179/2‑100 (MQ=255)                                                                                          gcgcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTCCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGt           >  2:1330423/1‑100 (MQ=255)                                                                                           cgcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTa          >  1:1811035/1‑100 (MQ=255)                                                                                           cgcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTa          <  1:1099269/100‑1 (MQ=255)                                                                                           cgcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACAACCAGCACTGCGACAATCAGGTa          >  2:575924/1‑100 (MQ=255)                                                                                           cgcTCCCGCAGGAGCACCGAATAGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTAGGAAAACACCACCAACACTGCGACAATCAGGTa          >  2:83968/1‑100 (MQ=255)                                                                                            gcTCCCGCAGTAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTCTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAgg            <  2:315152/97‑1 (MQ=255)                                                                                            gcTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTaa         >  1:2222724/1‑100 (MQ=255)                                                                                               cccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCg      >  1:232493/1‑100 (MQ=255)                                                                                               cccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCg      >  2:832301/1‑100 (MQ=255)                                                                                               cccGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGATTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGAGACAATCAGGTaa         >  2:1772819/1‑97 (MQ=255)                                                                                                ccGCAGGAGCGCCGAATGGATTAGCGACCGAAGAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGa     >  1:2885974/1‑100 (MQ=255)                                                                                                  gCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGACTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATa   >  2:2048787/1‑100 (MQ=255)                                                                                                   cAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATaa  >  2:2165198/1‑100 (MQ=255)                                                                                                    |                                                                                                   GCAGATGAAGCGAGCCATGCTTTCACCACGTTATACAACGCGAAGTTCTGATTTACTGAGGGTCAAATAAATATAGCGGCAGGAAAAAAGCGCTCCCGCAGGAGCGCCGAATGGATTAGCGACCGAACAGATCACGTTTTTTCGCTTTAAACGGCTGGGAAATCACCACCAGCACTGCGACAATCAGGTAAGCGATAA  >  NC_000913/1232385‑1232582

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