Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,231,864	C→T	5.9%	R227C (CGT→TGT)	umuC →	translesion error‑prone DNA polymerase V subunit; DNA polymerase activity

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,231,864	0	C	T	5.9%	223.3 / 3.1	68	R227C (CGT→TGT)	umuC	translesion error‑prone DNA polymerase V subunit; DNA polymerase activity
Reads supporting (aligned to +/- strand):  ref base C (31/33);  new base T (2/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

TCTCCCCGTGGATGACGTCTGGGGGATTGGACGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCC  >  NC_000913/1231719‑1232001                                                                                                                                                  |                                                                                                                                          tctcCCCGTGGATGACGTCTGGGGGATTGGACGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCg                                                                                                                                          <  8:257411/147‑1 (MQ=255) tctcCCCGTGGATGACGTCTGGGGGATTGGACGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCg                                                                                                                                          >  7:257411/1‑147 (MQ=255)     cccGTGGATGACGTCTGGGGGATTGGACGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGa                                                                                                                                       <  6:212530/146‑1 (MQ=255)     cccGTGGATGACGTCTGGGGGATTGGACGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAAc                                                                                                                                     >  4:294551/1‑148 (MQ=255)         tGGATGACGTCTGGGGGATTGGACGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTgcg                                                                                                                                 >  2:312046/1‑148 (MQ=255)               aCGTCTGGGGGATTGGACGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGa                                                                                                                           <  7:62822/148‑1 (MQ=255)                   cggggggCTTGGCCGGCGGCTCACCAAAAAACCGGACCCGCAGGGGAACCAAACCCGTCGCGAATTGGCGGAAACCGCTATCCCGGTTATCCGGAAACCATTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAAccc                                                                                                                       <  1:262467/146‑1 (MQ=255)                     gggggATTGGACGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGa                                                                                                                                       <  1:241778/130‑1 (MQ=255)                     gggggATTGGACGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGa                                                                                                                                       >  2:241778/1‑130 (MQ=255)                     gggggATTGGACGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTg                                                                                                                     >  8:258234/1‑148 (MQ=255)                      ggggATTGGACGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTg                                                                                                                     <  3:156546/147‑1 (MQ=255)                       gggATTGGACGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGtt                                                                                                                   <  2:168699/148‑1 (MQ=255)                              gACGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACt                                                                                                            >  4:277089/1‑148 (MQ=255)                                     gATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCa                                                                                                               >  6:277547/1‑138 (MQ=255)                                     gATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCa                                                                                                               <  5:277547/138‑1 (MQ=255)                                     gATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAagag                                                                                                     <  5:124811/148‑1 (MQ=255)                                       tCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGtt                                                                                                   >  6:277721/1‑148 (MQ=255)                                            aaaaaaCTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCAc                                                                                              >  5:310366/1‑148 (MQ=255)                                              aaaaCTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTgcg                                                                                                                                 <  3:13603/111‑1 (MQ=255)                                              aaaaCTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTgcg                                                                                                                                 >  4:13603/1‑111 (MQ=255)                                              aaaaCTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGAGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACcg                                                                                            <  5:210417/148‑1 (MQ=255)                                               aaaCTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACcga                                                                                           >  5:107260/1‑148 (MQ=255)                                                  cTGGACGCGAGGGGGATCAACACCGTGCTCGATTTGGCGGATACAGATATCCGGATTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACcgacga                                                                                        <  5:216047/148‑1 (MQ=255)                                                             ggggATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAAtt                                                                             >  6:1423/1‑148 (MQ=255)                                                              gggATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAAtt                                                                             <  7:7954/147‑1 (MQ=255)                                                              gggATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTa                                                                            >  6:42860/1‑148 (MQ=255)                                                                  tCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTg                                                                        >  2:241678/1‑148 (MQ=255)                                                                    aaaaCCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGtt                                                                      <  6:107260/148‑1 (MQ=255)                                                                      aaCCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTcc                                                                    <  1:312046/148‑1 (MQ=255)                                                                                gATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGttt                                                           <  3:294551/147‑1 (MQ=255)                                                                                gATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGttt                                                           <  5:89252/147‑1 (MQ=255)                                                                                     ggCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTAGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGa                                                      <  8:167549/147‑1 (MQ=255)                                                                                        ggATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGc                                                  <  1:241678/148‑1 (MQ=255)                                                                                           tACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCAt                                                >  5:64446/1‑147 (MQ=255)                                                                                             cAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCAc                                             >  2:21892/1‑148 (MQ=255)                                                                                             cAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCAc                                             >  8:168769/1‑148 (MQ=255)                                                                                             cAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCAc                                             >  3:278626/1‑148 (MQ=255)                                                                                                atatCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGa                                          <  4:278626/148‑1 (MQ=255)                                                                                                atatCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGa                                          >  5:108159/1‑148 (MQ=255)                                                                                                  atCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGAGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGa                                                      >  4:148301/1‑134 (MQ=255)                                                                                                  atCAGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTTCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGa                                                      <  3:148301/134‑1 (MQ=255)                                                                                                    ccGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGt                                                                                                    >  7:178212/1‑86 (MQ=255)                                                                                                    ccGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGt                                                                                                    <  8:178212/86‑1 (MQ=255)                                                                                                        tttATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGt                                  >  2:221802/1‑148 (MQ=255)                                                                                                           aTCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCg                                <  5:1423/147‑1 (MQ=255)                                                                                                           aTCCGTAAACATTTTAATGGCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGATTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCg                                <  2:42113/147‑1 (MQ=255)                                                                                                             ccGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATg                             >  1:199214/1‑148 (MQ=255)                                                                                                                tAAACATTTTAATGTCGTGCTCGAAAGAACGGTGTGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCgg                          >  4:136079/1‑148 (MQ=255)                                                                                                                tAAACATTTTAATGTCGTGCTCGAAAGAACGGTGTGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCgg                          >  6:123822/1‑148 (MQ=255)                                                                                                                 aaaCATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTAt                                      >  7:108684/1‑135 (MQ=255)                                                                                                                 aaaCATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTAt                                      <  8:108684/135‑1 (MQ=255)                                                                                                                 aaaCATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCg                                <  1:216344/141‑1 (MQ=255)                                                                                                                 aaaCATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCg                                >  2:216344/1‑141 (MQ=255)                                                                                                                  aaCATTTTAATGTCGTGCTCGAAAGAACGGTGTGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGc                         <  3:136079/147‑1 (MQ=255)                                                                                                                  aaCATTTTAATGTCGTGCTCGAAAGAACGGTGTGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGc                         <  5:123822/147‑1 (MQ=255)                                                                                                                  aaCATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGc                         <  1:21892/147‑1 (MQ=255)                                                                                                                       tttAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGc                                                  <  7:121945/117‑1 (MQ=255)                                                                                                                       tttAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGc                                                  <  4:89508/117‑1 (MQ=255)                                                                                                                       tttAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGc                                                  >  8:121945/1‑117 (MQ=255)                                                                                                                       tttAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGc                                                  >  3:89508/1‑117 (MQ=255)                                                                                                                          aaTGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGa                                                      >  4:7640/1‑110 (MQ=255)                                                                                                                          aaTGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGa                                                      <  3:7640/110‑1 (MQ=255)                                                                                                                            tGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGt              <  3:277089/148‑1 (MQ=255)                                                                                                                             gTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTg               >  8:145626/1‑146 (MQ=255)                                                                                                                              tCGTGCTCGATAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACGGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTa             <  6:310366/147‑1 (MQ=255)                                                                                                                                   cTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTAc        <  5:42860/147‑1 (MQ=255)                                                                                                                                     cGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCt     >  4:1166/1‑148 (MQ=255)                                                                                                                                        aaGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGcc  >  2:82735/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                          TCTCCCCGTGGATGACGTCTGGGGGATTGGACGGCGGATCAGCAAAAAACTGGACGCGATGGGGATCAAAACCGTTCTCGATTTGGCGGATACAGATATCCGGTTTATCCGTAAACATTTTAATGTCGTGCTCGAAAGAACGGTGCGTGAACTGCGCGGCGAACCCTGTTTGCAACTGGAAGAGTTTGCACCGACGAAGCAGGAAATTATCTGTTCCCGCTCGTTTGGTGAACGCATCACGGATTATCCGTCGATGCGGCAGGCCATTTGTAGTTACGCTGCC  >  NC_000913/1231719‑1232001

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