Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	599,576	0	T	G	41.8%	104.3 / 11.3	55	V288G (GTG→GGG)	cusA	copper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (0/32);  new base G (23/0);  total (23/32)
Fisher's exact test for biased strand distribution p-value = 5.36e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599427‑599725                                                                                                                                                      |                                                                                                                                                       aaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtg                                                                                                                                                       <  2:347184/151‑1 (MQ=255)   cGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCCTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGggggggg                                                                                                                                                     >  1:875893/1‑150 (MQ=255)    gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGggggt‑g                                                                                                                                                   >  2:835904/1‑151 (MQ=255)    gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtgggg                                                                                                                                                    >  2:174449/1‑149 (MQ=255)    gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGAGgggggggg                                                                                                                                                    >  2:447847/1‑149 (MQ=255)     cTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGggggggggg                                                                                                                                                   >  2:905685/1‑148 (MQ=255)        gacgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GGTcc                                                                                                                                               >  1:684636/1‑151 (MQ=255)         acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCt                                                                                                                                              <  1:710492/151‑1 (MQ=255)           gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGggggggggggg                                                                                                                                                     >  1:341164/1‑142 (MQ=255)               ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGa                                                                                                                                          <  2:1100049/149‑1 (MQ=255)               ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCACCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGAt                                                                                                                                         <  1:160442/150‑1 (MQ=255)                   tCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGG‑GGTCCTGCGAACCgg                                                                                                                                     >  1:530728/1‑150 (MQ=255)                     acaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCaa                                                                                                                                  <  2:1097437/151‑1 (MQ=255)                         cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCaaaaa                                                                                                                               <  2:224965/150‑1 (MQ=255)                                 aaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GATCCTGCGATCCGGCAAAAACGccc                                                                                                                          >  1:1112864/1‑147 (MQ=255)                                  aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGaa                                                                                                                     <  2:488277/151‑1 (MQ=255)                                    ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GGTTCCGCGATCCGGCAAAAAAGCCCCCGAAgg                                                                                                                   >  1:791981/1‑150 (MQ=255)                                     cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GGTCCTGCGGTCCGGCAAAAAAGCCCgggga                                                                                                                     >  1:54613/1‑144 (MQ=255)                                       gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGAt                                                                                                                <  1:681667/151‑1 (MQ=255)                                       gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGAt                                                                                                                <  1:758422/151‑1 (MQ=255)                                         gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATc                                                                                                               <  1:729317/150‑1 (MQ=255)                                         gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATc                                                                                                               <  2:675272/150‑1 (MQ=255)                                          aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGG‑GGTTCCGCGGTCCGGGAAAAAAGCCCGAGAAGTGAttg                                                                                                              >  2:158807/1‑148 (MQ=255)                                          aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GGTCCTGCGGTCCGGCAAAAAAGCCCGAGAAggg                                                                                                                  >  2:127322/1‑144 (MQ=255)                                          aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GGTACCGCGGTCCGGCCAAAAAGCCCGAGAAGTg                                                                                                                  >  1:1057409/1‑146 (MQ=255)                                           aaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgcc                                                                                                            <  1:815029/151‑1 (MQ=255)                                               ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCgccgcc                                                                                                         <  1:150107/150‑1 (MQ=255)                                               ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GAGCCTGCGGACCGGCAAAAACCGCCGAGGAACGATcggc                                                                                                            >  2:953618/1‑145 (MQ=255)                                                gCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGt                                                                                                       <  2:224099/151‑1 (MQ=255)                                                     cccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAgg                                                                                                  <  1:384342/151‑1 (MQ=255)                                                      ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGT‑GGTCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGGGGAgg                                                                                                  >  2:512056/1‑150 (MQ=255)                                                           tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACaaa                                                                                             <  1:535163/150‑1 (MQ=255)                                                            aTCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GGTCCCGCGGTCCGGCAAAAAAGCCCGCGAAGGGATCGCCGCCGGGAAGGACa                                                                                               >  1:809420/1‑147 (MQ=255)                                                                gcgcGATGTTGCGAAGGTCCAGATTGGTCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGa                                                                                        <  2:669792/150‑1 (MQ=255)                                                                   cgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACg                                                                                    <  1:1012461/151‑1 (MQ=255)                                                                      tGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGCGGAGGACAAACTGGAAACGCTg                                                                                 >  1:1033461/1‑151 (MQ=255)                                                                        ttGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GGTCCCGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGGCAAACTGGgaa                                                                                      >  1:176638/1‑144 (MQ=255)                                                                          gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGaaa                                                                              <  1:309122/150‑1 (MQ=255)                                                                            gAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGGCAAACTGGAAACGCTGaaa                                                                              >  2:774429/1‑148 (MQ=255)                                                                                  ccAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTg                                                                     <  1:540895/151‑1 (MQ=255)                                                                                  ccAGAGTGGCCCGGAGATGCCCCGGGGGATTGCCGAACTAAAGGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAACCGCCGGAGAAGTGATCGCCGCCGTGAAGGACAAGCTGGAAACGCTGAAAAGTAGTCTg                                                                     <  1:1087588/151‑1 (MQ=255)                                                                                     gATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGTGGGTACT‑CGCGCCGGCAAAAACGGCCCAGAAGGGCTCGCCGCCGTGGAGGGCCAAATCGAAACGCTGAAAAGAAGTATGCCg                                                                  >  2:598003/1‑151 (MQ=255)                                                                                       ttGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGGGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCgg                                                                 <  1:993549/150‑1 (MQ=255)                                                                                        tGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GATCCTGCGATCCGGCAAAAAAGCCCGCGAAGTGATCGCCGCCCGGAAGGACAAACTGGGAACGCTGAAAAGGAGTCTGCCgg                                                                 >  2:887353/1‑149 (MQ=255)                                                                                         ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGG‑GGTCCTTCGGTTCCGCAAAAACGCCCCAGCAGGGGTCGCCGCCGTGGAGGACACACTGGGAAAGCTGGAAAGTAGTCTGCCgg                                                                 >  1:124375/1‑148 (MQ=255)                                                                                          gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAgg                                                             <  2:623462/151‑1 (MQ=255)                                                                                          gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAgg                                                             <  2:1089266/151‑1 (MQ=255)                                                                                            ccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCg                                                           <  1:173097/151‑1 (MQ=255)                                                                                              ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTg                                                         <  1:32921/151‑1 (MQ=255)                                                                                                agaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GATCCTGCGATCCGGCAAAAAAGCCCGGGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTgg                                                        >  2:219385/1‑150 (MQ=255)                                                                                                  aTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGgaga                                                     <  2:112148/151‑1 (MQ=255)                                                                                                       ccGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGtt                                                <  1:651049/151‑1 (MQ=255)                                                                                                         ggggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GATCCTTCGAGCCGGCAAAAACGCCCGAGGAGGGGACGCCGCCGGGAAGGGCAAACTGGAAAAGCTGCAAAGAAGGCTGCCGGAAGGGGGGGGGAAAGtt                                                >  1:1069951/1‑149 (MQ=255)                                                                                                          gggCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGACGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTaca                                             <  2:124545/151‑1 (MQ=255)                                                                                                              aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GGTCCTTCGATCCGGCAAAAAAGCCCGAGAAGCGATCGCCGCCGTGAAGGACAAACTGGAAAAGCTGAAAAGAAGTCTGCCGGGAGGGGGGGGGGTAGtttaaaa                                           >  1:944619/1‑144 (MQ=255)                                                                                                                     aaCTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGG‑GATCCTGCGATCCGGCAAAAACGCCCGAGGAGTGATCGCCGCCGCGAAGGACAAAATGGAAACGCTGAAAAGTAGTCTGGCGGAAGGGGGGGGGGTAGGTACAACATACGAACg                                  >  1:955876/1‑151 (MQ=255)                                                                                                                         aaaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGca                                <  2:1108370/149‑1 (MQ=255)                                                                                                                           aCGGCGAAGGCGAAGTGGCGGGCGGGGGGGG‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAAAGCTGAAAAGGAGTCCGCCGGAAGGCGTGGGGGTAGTTACAACATAAGATCGcagcc                             >  1:480044/1‑150 (MQ=255)                                                                                                                                  aggcgaagTGGCGGGCGGGGGGGG‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAAAGCTGAAAAGTAGTCTGCCGGAAGGCGGGGAGATAGTTACAACATACGATCGCAGCCAGCTCAtt                     >  2:360767/1‑151 (MQ=255)                                                                                                                                   ggcgaagTGGCGGGCGGGGTGGT‑GATTCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTg                    <  2:49755/151‑1 (MQ=255)                                                                                                                                   ggcgaagTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCAtt                     <  1:7511/150‑1 (MQ=255)                                                                                                                                       aagTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACcg                <  2:1113002/151‑1 (MQ=255)                                                                                                                                                    ggtggt‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAc   >  2:1072928/1‑151 (MQ=255)                                                                                                                                                     gtggt‑gATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAcc  >  1:1070779/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                       AACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGT‑GATCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599427‑599725

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