Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	944,073	C→T	5.9%	T14M (ACG→ATG)	dmsC →	dimethyl sulfoxide reductase, anaerobic, subunit C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	944,073	0	C	T	5.9%	217.6 / 2.9	68	T14M (ACG→ATG)	dmsC	dimethyl sulfoxide reductase, anaerobic, subunit C
Reads supporting (aligned to +/- strand):  ref base C (33/31);  new base T (2/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

CCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATGGGCATTGGCTTTATCGCCTCTAT  >  NC_000913/943927‑944220                                                                                                                                                   |                                                                                                                                                    ccGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACg                                                                                                                                                    <  2:323640/148‑1 (MQ=255)      aGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTc                                                                                                                                              <  4:37736/149‑1 (MQ=255)                      ccGTATATTGTGATCAAACCGAATGCCAATAGCCGCCCGACCGTGTATACCACTGGCTATCTGGCAAAGCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAg                                                                                                                              <  2:216056/149‑1 (MQ=255)                       cGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAgg                                                                                                                             <  5:396353/149‑1 (MQ=255)                           tatTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCATGGTCTTCGGGCAATGTGTAGCAggtggt                                                                                                                         >  8:65113/1‑149 (MQ=255)                           tatTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCATGGTCTTCGGGCAATGTGTAGCAggtggt                                                                                                                         >  6:222756/1‑149 (MQ=255)                              tgtgATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGttt                                                                                                                       >  8:161669/1‑148 (MQ=255)                                tgATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTa                                                                                                                     <  2:405965/148‑1 (MQ=255)                                  aTCAAACCCAATGCCAATAGGCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATc                                                                                                                   <  7:76471/148‑1 (MQ=255)                                     aaaCCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTc                                                                                                               >  1:361003/1‑149 (MQ=255)                                         ccAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTAt                                                                                                                    >  1:413286/1‑140 (MQ=255)                                         ccAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTAt                                                                                                                    <  2:413286/140‑1 (MQ=255)                                              gCCAATAGCCGCCCGCCCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCCCGGTCTTCGGGCAATGTGTAGCAGGGGGTTTTATCGTTCTGGCTTTgg                                                                                                      >  6:329105/1‑149 (MQ=255)                                                 aaTAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATgtgt                                                                                                                                   <  6:209570/117‑1 (MQ=255)                                                 aaTAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATgtgt                                                                                                                                   >  5:209570/1‑117 (MQ=255)                                                        gccCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCa                                                                                             >  2:141529/1‑148 (MQ=255)                                                           cGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAgg                                                                                         >  4:44141/1‑149 (MQ=255)                                                                 gggATACAACTGGCTATCTGGCAAACCCTAAGGAGGAGTGAGATGGGAAGTGGATGGGATGAATGGCCGCTGATGATCTTCACTGTCTTCGGGCAATGAGTATCAGGTGGTTTTATCGTGATGGCTTTGGCGCTGCTCAAAGGCGACCt                                                                                   >  2:354414/1‑149 (MQ=255)                                                                   gATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGc                                                                                 <  1:141529/149‑1 (MQ=255)                                                                      aCCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGgcgc                                                                                                   >  8:248810/1‑128 (MQ=255)                                                                      aCCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGgcgc                                                                                                   <  7:248810/128‑1 (MQ=255)                                                                      aCCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAg                                                                              >  5:126408/1‑149 (MQ=255)                                                                       cgaCTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTgg                                                                                                            <  6:60485/116‑1 (MQ=255)                                                                         aCTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAg                                                                              >  3:261327/1‑146 (MQ=255)                                                                              cTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGcc                                                                      <  6:126408/149‑1 (MQ=255)                                                                               tATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAAtgt                                                                                                                                     <  6:146658/85‑1 (MQ=255)                                                                               tATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAAtgt                                                                                                                                     >  5:146658/1‑85 (MQ=255)                                                                               tATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGcc                                                                      >  4:314981/1‑148 (MQ=255)                                                                               tATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGcc                                                                      >  2:188497/1‑148 (MQ=255)                                                                                aTCTGGCAAACCCGAAGGAGGTGTTAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGAGCTGCTCAAAGGCTCCCTGCGAGAAGAAGcc                                                                      >  6:73172/1‑147 (MQ=255)                                                                                aTCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGcc                                                                      >  5:158467/1‑147 (MQ=255)                                                                                 tCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCATGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCca                                                                    <  7:65113/148‑1 (MQ=255)                                                                                 tCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCATGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCca                                                                    <  5:222756/148‑1 (MQ=255)                                                                                     gCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCagcagc                                                               >  3:155080/1‑149 (MQ=255)                                                                                      cAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCg                                                              <  2:361003/149‑1 (MQ=255)                                                                                       aaaCCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAggtg                                                                                                                           <  7:236188/87‑1 (MQ=255)                                                                                       aaaCCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAggtg                                                                                                                           >  8:236188/1‑87 (MQ=255)                                                                                       aaaCCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCgt                                                             >  4:324695/1‑149 (MQ=255)                                                                                       aaaCCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCgt                                                             >  4:251716/1‑149 (MQ=255)                                                                                          cccGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCgtgt                                                           >  1:39192/1‑148 (MQ=255)                                                                                            cGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGtt                                                                                                                        >  6:308924/1‑85 (MQ=255)                                                                                            cGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGtt                                                                                                                        <  5:308924/85‑1 (MQ=255)                                                                                               aggaggTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGc                                                     <  3:251716/149‑1 (MQ=255)                                                                                                       tgAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCcag                                                                   >  1:42724/1‑127 (MQ=255)                                                                                                       tgAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCcag                                                                   <  2:42724/127‑1 (MQ=255)                                                                                                       tgAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGt                                             >  1:36342/1‑149 (MQ=255)                                                                                                            tGGGATGTGGATGGCATGAATGGCCGCTGATGATCTTCACTGTATTCTGGCAATGTGTAGCAGGTGGGTTTATCGTTCTGGCTGTGGGGCTGCTCAAAGGAGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTgg                                         <  2:260287/148‑1 (MQ=255)                                                                                                             gggAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGtt                                       >  2:348945/1‑149 (MQ=255)                                                                                                              ggAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGtt                                       <  1:348945/148‑1 (MQ=255)                                                                                                               gAAGTGGCTGGCATGAATGGCCGCTGATGTTCTTCTCGGGCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGcta                                                                               <  3:218883/107‑3 (MQ=255)                                                                                                                  gTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATgg                                  <  3:44141/149‑1 (MQ=255)                                                                                                                     gATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTAAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGt                                >  5:110427/1‑148 (MQ=255)                                                                                                                          cATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGAt                          >  3:49173/1‑149 (MQ=255)                                                                                                                            ttaaTGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCaaat                                                                                          >  3:36675/3‑82 (MQ=255)                                                                                                                            ttaaTGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAg                                                                                          <  4:36675/81‑1 (MQ=255)                                                                                                                            tGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATgg                        <  2:39192/149‑1 (MQ=255)                                                                                                                              aaTGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCa                                                <  8:381183/123‑1 (MQ=255)                                                                                                                              aaTGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCa                                                >  7:381183/1‑123 (MQ=255)                                                                                                                                 ggCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATGGGCAtt                   >  6:335871/1‑149 (MQ=255)                                                                                                                                   ccGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATGGGCATTgg                 <  4:49173/149‑1 (MQ=255)                                                                                                                                         atgatCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGt                                        >  6:254123/1‑120 (MQ=255)                                                                                                                                         atgatCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGt                                        <  5:254123/120‑1 (MQ=255)                                                                                                                                         atgatCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATGGGCATTGGCTTTAt           <  5:335871/149‑1 (MQ=255)                                                                                                                                          tgatCTTCACGGTATGCGCTCAATTTGTAGCAGGTGGTTTTATCGTTCTTGCTTTGGCGCTGCTCAAAGGCGACATGCGAGCAGAAGCACAGCATCGTGTTATCGCCTTCATGTTTGTTTTATGGGTGCTGATGGGCATTGGCTTTATc          <  5:73172/149‑1 (MQ=255)                                                                                                                                             tCTTCACGGTCTTCGGGCATTGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGACGGCATGTTTGGTTTATGGGTGCTGATGGGCATTGTCTTTATCGcc       >  6:95832/1‑149 (MQ=255)                                                                                                                                                 cACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATGGGCATTGGCTTTATCGc        <  3:324695/144‑1 (MQ=255)                                                                                                                                                  aCGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATGGGCATTGGCTTTATCGCCTCTa   <  4:155080/148‑1 (MQ=255)                                                                                                                                                  aCGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATGGGCATTGGCTTTATCGCCTCTAt  >  3:343958/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATGGGCATTGGCTTTATCGCCTCTAT  >  NC_000913/943927‑944220

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