Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,236,716	C→T	6.0%	N59N (AAC→AAT)	ygjQ →	DUF218 superfamily protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,236,716	0	C	T	6.0%	206.3 / 2.7	66	N59N (AAC→AAT)	ygjQ	DUF218 superfamily protein
Reads supporting (aligned to +/- strand):  ref base C (30/32);  new base T (2/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

CGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCAT  >  NC_000913/3236570‑3236858                                                                                                                                                   |                                                                                                                                               cGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGCCCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGt                                                                                                                                              >  6:105187/1‑149 (MQ=255) cGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGt                                                                                                                                              <  5:164173/149‑1 (MQ=255) cGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGt                                                                                                                                              >  2:172133/1‑149 (MQ=255)   tATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTgg                                                                                                                                            <  3:228128/149‑1 (MQ=255)     tttGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGAGCATCCTTATCGCAGCTCGCTTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTggga                                                                                                                                          <  6:75980/149‑2 (MQ=255)          ttttCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGCTCCGGCGCGTAACGTGGGTTTAt                                                                                                                                      <  1:142764/148‑1 (MQ=255)                 gCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTg                                                                                                                                 <  1:172133/146‑1 (MQ=255)                     gACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGa                                                                                                                          <  8:30355/149‑1 (MQ=255)                          tAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAATGTGGGTTTATTGCTGGGGGCGAGGCCg                                                                                                                     >  6:150704/1‑149 (MQ=255)                          tAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAATGTGGGTTTATTGCTGGGGGCGAGGCCg                                                                                                                     >  6:147175/1‑149 (MQ=255)                          tAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCg                                                                                                                     >  7:15315/1‑149 (MQ=255)                          tAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCg                                                                                                                     >  8:195729/1‑149 (MQ=255)                          tAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCGGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCg                                                                                                                     >  2:52104/1‑149 (MQ=255)                            aaGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCggg                                                                                                                   <  5:76555/149‑1 (MQ=255)                            aaGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCggg                                                                                                                   <  3:108427/149‑1 (MQ=255)                            aaGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCg                                                                                                                     <  7:56749/147‑1 (MQ=255)                              gATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGt                                                                                                                  >  6:63081/1‑148 (MQ=255)                               aTTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTaa                                                                                                                >  1:107039/1‑149 (MQ=255)                                ttGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAAt                                                                                                               >  8:242023/1‑149 (MQ=255)                                        ttACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACtt                                                                                                       >  7:90533/1‑149 (MQ=255)                                           cTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAgg                                                                                                                        >  3:208541/1‑129 (MQ=255)                                           cTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAgg                                                                                                                        <  4:208541/129‑1 (MQ=255)                                           cTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGGCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTAc                                                                                                    <  6:205663/149‑1 (MQ=255)                                            tGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTAcc                                                                                                   >  4:63190/1‑149 (MQ=255)                                              ttACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGCGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTAccc                                                                                                  >  3:34443/1‑148 (MQ=255)                                                     tcgCAGGAGCGACCATCCTTCCCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGTGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTAtt                                                                                          <  5:199214/147‑1 (MQ=255)                                                       gCAGGAGCGCCCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTACTCGCTACTTTACCCGCCGTATTg                                                                                         >  4:28320/1‑148 (MQ=255)                                                       gCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCAGGAGCGATGTCAACGCCGTTCCGGCCACTAACGTGGCATTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGa                                                                                        >  7:173885/1‑149 (MQ=255)                                                         aGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATa                                                                                      <  7:10543/149‑1 (MQ=255)                                                          ggAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATAc                                                                                     >  6:54563/1‑149 (MQ=255)                                                                    tCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCg                                                                           >  1:167442/1‑149 (MQ=255)                                                                      cTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCCGCAGCGtt                                                                         <  5:63081/149‑1 (MQ=255)                                                                      cTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGtt                                                                         <  7:195729/149‑1 (MQ=255)                                                                         aTCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCg                                                                                                                           <  3:16855/96‑1 (MQ=255)                                                                         aTCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCg                                                                                                                           >  4:16855/1‑96 (MQ=255)                                                                           cGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTAc                                                                     >  4:40476/1‑148 (MQ=255)                                                                             ccGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCAc                                                                  <  2:167442/149‑1 (MQ=255)                                                                               gATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGc                                                                >  7:59410/1‑149 (MQ=255)                                                                               gATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGc                                                                <  2:107039/149‑1 (MQ=255)                                                                                atagcgTAATGGTCAACGACAGCAAACAGCTGACATGGAGCGATGTCAACGACGTTCCGGCGAGTAACGTTGGTTTATTGCTGGGGGAGATGCCGGGTAATAGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTTTACCACGcc                                                               >  5:214687/4‑149 (MQ=255)                                                                                     gTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAAGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTAc                                                                     >  6:186524/1‑138 (MQ=255)                                                                                     gTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAAGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTAc                                                                     <  5:186524/138‑1 (MQ=255)                                                                                      tAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAATGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGCCGGaaaa                                                         <  5:150704/149‑1 (MQ=255)                                                                                      tAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAATGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGCCGGaaaa                                                         <  5:147175/149‑1 (MQ=255)                                                                                      tAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGCCGGaaaa                                                         <  8:15315/149‑1 (MQ=255)                                                                                      tAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGCCGGaa                                                           >  3:212256/1‑147 (MQ=255)                                                                                       aaTGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAAt                                                                                                               <  3:168243/94‑1 (MQ=255)                                                                                           gTCAACGCCAGTAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGCCGGAAAAGTGa                                                     >  6:233247/1‑148 (MQ=255)                                                                                           gTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGCCGGAAAAGTGaa                                                    <  5:54563/149‑1 (MQ=255)                                                                                              aacgaatgCACACAGCTGTCCTGGATCTATGTCAAGGCCTTTCCTGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGGTACTTTACCCGTCGTATTGATACCGCAGCAGCGTTGTACCACGCCGGAAAAGTGAAAt                                                  <  8:173885/141‑1 (MQ=255)                                                                                               aCGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGCCGGAAAAGTGAAATgg                                                >  7:79651/1‑149 (MQ=255)                                                                                               aCGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGCCGGAAAAGTGAAATgg                                                <  6:216666/149‑1 (MQ=255)                                                                                                      cAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTgg                                         >  4:157365/1‑149 (MQ=255)                                                                                                       aaaCAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTgg                                         >  3:213596/1‑148 (MQ=255)                                                                                                           aGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGt                                    >  5:158611/1‑149 (MQ=255)                                                                                                             cTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTgg                                  <  5:76997/149‑1 (MQ=255)                                                                                                                    ggAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAAc                           >  7:135337/1‑149 (MQ=255)                                                                                                                     gAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACg                          <  5:105187/149‑1 (MQ=255)                                                                                                                               aaCGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACGACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTaaaaa                <  8:59410/149‑1 (MQ=255)                                                                                                                                  gCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTa             <  3:63190/149‑1 (MQ=255)                                                                                                                                   ccGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTAc            <  3:51297/149‑1 (MQ=255)                                                                                                                                    cGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGCCGGAAAAGTGAAATGGctgctg                                          <  6:141710/118‑1 (MQ=255)                                                                                                                                    cGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGCCGGAAAAGTGAAATGGctgctg                                          >  5:141710/1‑118 (MQ=255)                                                                                                                                    cGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACg           <  5:233247/149‑1 (MQ=255)                                                                                                                                         cGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGAAACCGCAGCAGCGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGaa      <  7:88931/149‑1 (MQ=255)                                                                                                                                             gcgTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACAGCAGCAGCGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGAAAGTGGCGATAACTGACGTAAAAATTACGATGAAGCAt  >  3:25958/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                               CGTATTTGTTTTTCTCGCCGGACTCTAAAGATTGCCTGCTTACTGTTACTTGTGGCAGGAGCGACCATCCTTATCGCCGATCGCGTAATGGTCAACGCCAGCAAACAGCTGACCTGGAGCGATGTCAACGCCGTTCCGGCGCGTAACGTGGGTTTATTGCTGGGGGCGAGGCCGGGTAATCGCTACTTTACCCGACGTATTGATACCGCAGCAGCGTTGTACCACGCCGGAAAAGTGAAATGGCTGCTGGTAAGTGGCGATAACGGACGTAAAAATTACGATGAAGCAT  >  NC_000913/3236570‑3236858

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