Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,042,604	0	A	C	45.0%	96.1 / 16.2	60	V182G (GTG→GGG)	yqfA	hemolysin III family HylIII inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (26/7);  new base C (0/27);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 8.34e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCG  >  NC_000913/3042456‑3042748                                                                                                                                                      |                                                                                                                                                 gcgcGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCAcc                                                                                                                                                >  2:693613/1‑151 (MQ=255) gcgcGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCAcc                                                                                                                                                >  2:474333/1‑151 (MQ=255)  cgcgTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCAccc                                                                                                                                               <  2:774108/151‑1 (MQ=255)         tATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGc                                                                                                                                          <  1:271024/149‑1 (MQ=255)         tATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCg                                                                                                                                         <  1:660531/150‑1 (MQ=255)          aTCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGa                                                                                                                                        <  2:612138/150‑1 (MQ=255)            cAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAAt                                                                                                                                      <  2:936713/150‑1 (MQ=255)                ccTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAAcc                                                                                                                                 >  2:233979/1‑151 (MQ=255)                  tACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAACCAc                                                                                                                               >  1:532124/1‑151 (MQ=255)                   aCAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAACCAcg                                                                                                                              >  1:482236/1‑151 (MQ=255)                    cAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAACCAcgc                                                                                                                             >  2:198797/1‑151 (MQ=255)                     aaaGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAACCAcgcc                                                                                                                            >  1:293336/1‑151 (MQ=255)                     aaaGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAACCAcgcc                                                                                                                            >  2:197954/1‑151 (MQ=255)                              cccATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGc                                                                                                                   >  2:746855/1‑151 (MQ=255)                                 aTTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGcc                                                                                                                  >  2:832103/1‑149 (MQ=255)                                      gCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAgg                                                                                                           >  1:1119129/1‑151 (MQ=255)                                           cccAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCCCCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTaa                                                                                                        <  2:247408/149‑1 (MQ=255)                                           ccaaaTATACAAATAGACCGCCAGAAAGTGGCACCCCCTCCCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAAGCCCCCCGAGCGAATAAACCCCGCCGCCTACCGCCAGTAAGGTAAc                                                                                                       <  1:359416/147‑1 (MQ=255)                                           ccaaaTATAAAAAAAGACCGCCAAAAGGTGGCACACCCTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAAGGCGTTTGCAGAC‑GTAGAAAATCCCCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACg                                                                                                      <  1:184138/148‑1 (MQ=255)                                            ccAATATACAAATAGATCCCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGc                                                                                                     <  2:383662/151‑1 (MQ=255)                                                        tAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCcgcg                                                                                          >  2:580584/1‑150 (MQ=255)                                                             cGCCAAAAGGGGGCACACACTCCCGCCGAGCACGAAGCCGTGCCAGTGGGCTTGTTTGTAGGGAAGGCGTTTGCAGAC‑GTAAAAAACCCCCCCGAGCAAAAAACCCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCtt                                                                                     <  2:646855/150‑1 (MQ=255)                                                               caaGAAAGGGGCACACCCTCCCCCCGAGCACAAAGCGGTGCCAGAGGGCAGGGTTGTATGGAATGCGTTTGCAGAT‑GTAAAAAATCCCCCCGAGCGAAAAACCCACGCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTaa                                                                                   <  2:305064/148‑1 (MQ=255)                                                                    aaaatgCACACACTACCGCTGATCACGAAGCCGTGACAGATGGCATGGATATCTGGGATGCGTTTGCAGAC‑GTAGAAAATCACCGCGAGGGAATAAACGACGCCGCCTACCACCGGGAAAATAACGCTGCCAGCCTCAAGATTAACTg                                                                                >  1:459716/6‑148 (MQ=255)                                                                       ggccacacaCTACCGCCGGGCACAAGGCCGTCCGAGAGGCCGTGTTTTTATGAATTGCGTTGGCAGCC‑GGAGAAAACCCCCCGGGGCAAAAAACCCACCCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCAttt                                                                          <  2:593897/148‑1 (MQ=255)                                                                                 cccgcccgAGCAGGAGGCGGTGCGAGAGGCCAGGTTTTTGTGGAATGGTTTGGCGGCC‑GTAAAAAACCACCCCGGCCAAATAACCCACGCCGCCCACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTac                                                                <  2:979534/146‑1 (MQ=255)                                                                                  aCCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCAAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTa                                                                 >  2:742456/1‑149 (MQ=255)                                                                                   ccgccgAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCCCCCCGAGCGAATAACCCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAAtt                                                                  <  2:563026/147‑1 (MQ=255)                                                                                       cgAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTaccacc                                                            >  2:56703/1‑149 (MQ=255)                                                                                        gAGCAGGAGGCCGGGCAAGAGGGCGTGTTGGTATGAAAGGCGTTGGCACAC‑TTAAAAAACCCCCCCGAGCAAAAAACCCCCGCCCCCTCCCGCCAGTAAGTTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTAccacca                                                           <  1:212482/149‑1 (MQ=255)                                                                                           cACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCTTTTGCAGAC‑GTAAAAAATCCCCCCGAGCAAATAACCCCCGCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGa                                                       <  2:127491/150‑1 (MQ=255)                                                                                           cACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAc                                                      >  1:84572/1‑151 (MQ=255)                                                                                            aCGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAc                                                      >  2:908664/1‑150 (MQ=255)                                                                                              gAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAcagc                                                   >  2:48150/1‑151 (MQ=255)                                                                                               aggCCGTGCAAGAGGGCAGGGTTGTAGGAAGGGCTTTGGCAGAG‑GTAAAAAATCCCCCCGAGCAAAAAACCCCCGCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTaccac                                                             <  1:990477/138‑1 (MQ=255)                                                                                               agcccggcccAGAGGGCATGGTTTTAGGGAATGCTTTTGCAGAG‑GTAAAAAACCCCCCCGAGCGAAAAACCCCCGCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAcagc                                                   <  2:811234/142‑1 (MQ=255)                                                                                               aaGCCGGGCAAGATGGCTTGGTTGTAGGGAATGCTTTTGCAGAC‑GTAAAAAATCCCCCCGAGCAAATAACCCCCGCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAcagc                                                   <  2:456981/150‑1 (MQ=255)                                                                                                  cggTGCCAGAGGGCTTGTTTGTATGGAATGCGTTTGCAGAC‑GTAAAAAATCCCCCCGAGCAAATAACCCCCGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGAcagcca                                                 <  2:1153510/147‑1 (MQ=255)                                                                                                     tGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGccc                                             >  2:418948/1‑150 (MQ=255)                                                                                                       ccAGAGGGCATGTTTTTAGGGAATGCTTTTGCAGAC‑TTAAAAAACCCCCCCGAGCAAAAAACCCCCCCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCAt                                           <  1:929414/150‑1 (MQ=255)                                                                                                            tGGCTTGTTTTTTTGAAATGCGTTTGCAGAG‑TTAAAAAACCCCCCCGAGCGAAAAACCCCCCCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGcc                                       <  2:538098/149‑1 (MQ=255)                                                                                                             ggCAGGTTTGTATGGAAGGTTTTTGGAGAC‑GTAAAAACCCCCCCGGGGCGAATAACCCCCCCCGCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCg                                         <  2:220300/146‑1 (MQ=255)                                                                                                             gccATGTTTGTAGGAATTGTTTTTGCAGCC‑GAAAAAAACCCCCCCGGGCAAAAAACCCCCCCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGcc                                       <  1:604190/146‑1 (MQ=255)                                                                                                             gccATGGTTGTATGGATTGCTTTGGCAGCC‑GTAAAAAACCCCCCCGAGCAAAAAACCCCCCCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATc                                          <  2:220745/143‑1 (MQ=255)                                                                                                             gcaaTGTTTGTATGAAATG‑GTTTGCTGGCGGTAGAAAACCCCCCGGGGCAAAAAACCCCCGCCGCCTTCCGCCAGTAAGGTAACGTTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGa                                    <  2:896066/148‑1 (MQ=255)                                                                                                                  gtttGTAGGGAAGGCTTTGGCAGAC‑GAAAAAAACCCCCCCGAGCGAAAAACCCCCGCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGt                               <  2:1038330/149‑1 (MQ=255)                                                                                                                     tGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCAc                            >  2:379448/1‑151 (MQ=255)                                                                                                                        aTGAAATGCGTTTGCAGAC‑GTAAAAAATCCCCCCGAGCGAAAAACCCACGCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCa                          <  2:50402/150‑1 (MQ=255)                                                                                                                          ggAATGCGTTTGCAGAC‑GTAAAAAACCCCCCCGAGCAAAAAACCCCCGCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCaga                        <  2:272543/150‑1 (MQ=255)                                                                                                                          gaaggccttttGCGGCC‑GAAAAAAACCCCCCCGGGCAAAAAACCCCCCCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCagag                       <  1:120614/143‑1 (MQ=255)                                                                                                                             aTGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGAt                     >  2:943541/1‑150 (MQ=255)                                                                                                                                 gTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAAt                >  1:817837/1‑151 (MQ=255)                                                                                                                                 gTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAAt                >  1:531710/1‑151 (MQ=255)                                                                                                                                  tttGCAGCC‑TAAAAAAACCCCCCCGAGCAAAAAACCCCCCCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATa                    <  1:174409/146‑1 (MQ=255)                                                                                                                                  ttggCAGAC‑GAAAAAAACCCCCCCGGGCAAAAAACCCCCCCCCCCTCCCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATaa                   <  1:531947/144‑1 (MQ=255)                                                                                                                                   tttggggcggTAAAAAATCCCCCCGAGCAAAAAAACCGCGCCCCCTCCCGCCATTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATaaa                  <  2:316134/139‑1 (MQ=255)                                                                                                                                         aC‑GTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATc        >  2:257710/1‑151 (MQ=255)                                                                                                                                            gaaaaaaaCCCCCCCGAGCAAAAAACCCCCCCCCCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATc        <  2:520207/145‑1 (MQ=255)                                                                                                                                                aaaaTCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTgcg  >  1:219930/1‑151 (MQ=255)                                                                                                                                                 aaaTCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTgcg  >  1:386178/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                 GCGCGTCTTATCAGGCCTACAAAGGCATACCCATTACGCCTGCCCAATATACAAATAGATCGCCAGAAAGTGGCACACACTACCGCCGAGCACGAAGCCGTGCCAGATGGCATGGTTGTATGGAATGCGTTTGCAGAC‑GTAGAAAATCACCCCGAGCGAATAAACCACGCCGCCTACCGCCAGTAAGGTAACGCTGCCCGCCGCGAGCTTAACTGCCATTTCATAAATTACCACCAGCGACAGCCAGCCCATCGCCAGATAGGTCACCAGAGATAAAATTTTGAATCGGTGCG  >  NC_000913/3042456‑3042748

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