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

AAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599426‑599725                                                                                                                                                       |                                                                                                                                                        aaaCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggt                                                                                                                                                         <  2:137244/151‑1 (MQ=255)     gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggg                                                                                                                                                     >  1:153986/1‑149 (MQ=255)     gCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggg                                                                                                                                                     >  2:638732/1‑149 (MQ=255)          acgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCt                                                                                                                                               <  2:757762/151‑1 (MQ=255)            gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGc                                                                                                                                             <  2:569136/151‑1 (MQ=255)                   aTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGG‑TTCTGCGGTcccg                                                                                                                                      >  1:61436/1‑149 (MQ=255)                     cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGc                                                                                                                                     <  2:56403/150‑1 (MQ=255)                     cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCa                                                                                                                                    <  1:945912/151‑1 (MQ=255)                     cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGG‑TCCTGCGGTCCgg                                                                                                                                      >  1:411661/1‑149 (MQ=255)                        aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCaaaa                                                                                                                                 <  2:780690/151‑1 (MQ=255)                          cGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCaaaaa                                                                                                                                <  2:694881/150‑1 (MQ=255)                           gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACg                                                                                                                              <  1:75748/151‑1 (MQ=255)                           gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGA‑TCCTGCGATCCGGCaaaaa                                                                                                                                >  1:638825/1‑149 (MQ=255)                               tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGGGAAGTGGCGGGGGGGGGGGGGG‑TGCTGCGCGCCGGCaaaaaagacc                                                                                                                           >  2:605143/1‑145 (MQ=255)                                aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATTCCTG‑GATCCGGCAAAAACGccc                                                                                                                           >  1:465228/1‑149 (MQ=255)                                   aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGG‑GCCTGCGACCCCGCCAAAACCCCCgag                                                                                                                        >  2:477382/1‑149 (MQ=255)                                   aGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGA‑TCCTGCGATCCGGCAAAAACGccc                                                                                                                           >  1:714420/1‑146 (MQ=255)                                    gCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGG‑TCCTGCGGTCCGGCaaaaaacgcc                                                                                                                           >  2:460128/1‑140 (MQ=255)                                      cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGG‑ACCTGCGGTCCGGCaaaaaacccccggaaggg                                                                                                                   >  1:2650/1‑138 (MQ=255)                                         tGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGG‑TTCTGCGGTCCGGCCAAAAAGCCCggg                                                                                                                        >  1:497587/1‑141 (MQ=255)                                          gAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGA‑TCCTGCGATCCGGCAAAAACGCCCGGGGAGTGATcg                                                                                                               >  1:403617/1‑151 (MQ=255)                                           aaaaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGA‑TCCTGCGATCCCGGAAAAAAGCCCGAGAAggggt                                                                                                                 >  1:110106/1‑144 (MQ=255)                                           aaaaTGGCGCGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGA‑TCCCGCGAACCGGCAAAAAAACCCCAGGAGTggt                                                                                                                 >  2:288318/1‑146 (MQ=255)                                             aaTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATcgccg                                                                                                            <  2:565082/151‑1 (MQ=255)                                                ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGccgccg                                                                                                         <  2:158851/151‑1 (MQ=255)                                                     gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAg                                                                                                    <  2:85603/151‑1 (MQ=255)                                                            tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGG‑TCCTGCGATCCGGCAAAAAAGCCCGAGAAGTGATCGCCGCCGTGAAGGACaaa                                                                                              >  1:287541/1‑150 (MQ=255)                                                                 gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACt                                                                                            <  1:253739/147‑1 (MQ=255)                                                                 gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTgg                                                                                          <  2:594445/149‑1 (MQ=255)                                                                  cgcgATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGaaa                                                                                       <  1:59696/151‑1 (MQ=255)                                                                     gATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑Gc                                                                                   <  2:223000/151‑1 (MQ=255)                                                                       tGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGGGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTg                                                                                 >  1:268362/1‑151 (MQ=255)                                                                         ttGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTg                                                                                 >  1:590301/1‑149 (MQ=255)                                                                          tGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGA‑GCCTGCGCTCCGGCGAAAACGCCCGAGGAGGGAACGCCGCCGTTGAGGACTAACTGGaaaaggttaa                                                                                >  1:598215/1‑143 (MQ=255)                                                                           gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGaaaa                                                                             <  1:547008/151‑1 (MQ=255)                                                                                 gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGG‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAagtaggc                                                                       >  1:118335/1‑149 (MQ=255)                                                                                  tCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGG‑TCCCGCGATCCGGCAAAAAAGCCCCGGAAGGGAACGCCGCCGGGGAGGGCAAAATGGAAAA‑GCCCAAAagg                                                                           >  2:130541/1‑145 (MQ=255)                                                                                       aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAAGTAGTCTGCCgg                                                                 <  1:504415/151‑1 (MQ=255)                                                                                          ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGAGGCGGTGGGGGGGGGGG‑GCTTGTGAACCGGCAAAAACGACCGGCGAGGGATGGGCGGCGAGAAAGGCAAAATAGAAACAGCGGaaaacagaacagcggg                                                                 >  2:747939/1‑137 (MQ=255)                                                                                            cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAAGTAGTCTGCCGGAAGGc                                                            <  1:171779/151‑1 (MQ=255)                                                                                            cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCATGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAAGTAGTCTGCCGGAAGGc                                                            <  1:720189/151‑1 (MQ=255)                                                                                               ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGG‑TCCTGCGATCCGGCAAAAACGCCCGGGAAGTGATCGCCGCCGGGAAGGACAAACTGGGAAA‑GCTGAAAAGAAGTCTGCCGGAAgg                                                             >  1:397517/1‑147 (MQ=255)                                                                                                 agaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGGGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAAGTAGTCTGCCGGAAGGCGTgg                                                        >  1:391157/1‑150 (MQ=255)                                                                                                  gaTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAAGTAGTCTGCCGGAAGGCGTGga                                                       <  1:945107/150‑1 (MQ=255)                                                                                                   aTGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGG‑TCCTGCGATCCGGCAAAAAAGCCCGAGAAAGGATCGCCGCCCGGGAGGGCAAAATGGAAAC‑GCTGGCAAGTAGTCTGCCGGAAGGCgggg                                                        >  2:942134/1‑148 (MQ=255)                                                                                                     gcgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATa                                                   <  1:437256/151‑1 (MQ=255)                                                                                                              cATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaaca                                          <  2:113588/151‑1 (MQ=255)                                                                                                               aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACAt                                         <  2:403456/151‑1 (MQ=255)                                                                                                               aTTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCCGCGCGCCGGGAAAAAAGGCCGAGAAGGGGGCGCCGCCCGGAAGGGAAAACTGGAAAA‑GCTGAAAAGGAGTCTGCGCGAAGGCGGGGGGGTGGTTacaaaa                                          >  2:526353/1‑148 (MQ=255)                                                                                                                  gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATAc                                       <  2:933883/150‑1 (MQ=255)                                                                                                                  gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATAc                                       <  1:806276/150‑1 (MQ=255)                                                                                                                     gAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGA‑TCCTGCGGTCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAAGTAGTCTGCCGGAAGGCGGGGGGATAGTTACAACATACg                                      >  1:702491/1‑148 (MQ=255)                                                                                                                      aaCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATc                                   <  2:61427/150‑1 (MQ=255)                                                                                                                             cGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagccag                           <  1:288433/151‑1 (MQ=255)                                                                                                                                     gcgaagTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGGAAC‑GCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGa                   <  2:102192/151‑1 (MQ=255)                                                                                                                                      cgaagTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAc                  <  1:484771/151‑1 (MQ=255)                                                                                                                                               gggcgggGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc         >  1:59569/1‑151 (MQ=255)                                                                                                                                                ggcgggGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCg        <  2:116939/151‑1 (MQ=255)                                                                                                                                                   ggggTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc      >  2:880783/1‑150 (MQ=255)                                                                                                                                                   ggggTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACa     <  1:171731/151‑1 (MQ=255)                                                                                                                                                      gtggtgA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAcc  <  1:59695/151‑1 (MQ=255)                                                                                                                                                       tggtgA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAAcc  >  2:502469/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                        AAACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGA‑TCCTGCGATCCGGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGACAAACTGGAAAC‑GCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACAACC  >  NC_000913/599426‑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: ‑