Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,378,802	0	A	C	46.6%	85.8 / 17.5	58	V48G (GTA→GGA)	ampC	penicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (27/4);  new base C (0/27);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.23e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑C‑GTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTC  >  NC_000913/4378653‑4378950                                                                                                                                                             |                                                                                                                                                     ccccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑C‑GTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTac                                                                                                                                                     >  1:711133/1‑151 (MQ=255) ccccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑C‑GTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTac                                                                                                                                                     <  1:756094/151‑1 (MQ=255)           aaTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑C‑GTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTaccgccaccgcc                                                                                                                                           <  2:601356/151‑1 (MQ=255)            aTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑C‑GTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTaccgccaccgcc                                                                                                                                           >  2:546465/1‑150 (MQ=255)             tAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑C‑GTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCAt                                                                                                                                         >  2:889959/1‑151 (MQ=255)             tAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑C‑GTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCAt                                                                                                                                         >  2:503918/1‑151 (MQ=255)             tAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑C‑GTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCAt                                                                                                                                         >  2:310311/1‑151 (MQ=255)               gCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑C‑GTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCATAc                                                                                                                                       <  2:655936/151‑1 (MQ=255)                          aggcccGCCATAAATTGTTTTCCTGCCCAAACCTACCCAAAACAC‑C‑TTTTGTTGTGTGAGGGGTTTTTTTTGGGCGTTGCCCCCATACCCCAAGAT‑AAAG‑AAAAAAGTTTTCCC‑T‑GAATAAT‑TTCCCCCCACCGCCATACCCGGGATCttt                                                                                                                            <  1:49062/146‑1 (MQ=255)                          aaGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑C‑GTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCttt                                                                                                                            >  2:557796/1‑151 (MQ=255)                          aaGCACGCAATTAATTGTTTGGTTGCCCTAACTAAACTAAAACAA‑C‑TTTGGCTGTGGGAGGGGCTGTTTTTTGGGGAGGTCCGCATACCCCCAGAA‑AAAGAAATAAGGTTTACCT‑T‑A‑ATAAA‑TTCCCGCCACCGCCATACCCGGGATCt                                                                                                                              <  2:930153/149‑1 (MQ=255)                             cACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑C‑GTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTttgt                                                                                                                         >  2:296901/1‑151 (MQ=255)                               cgcaaatAATTGTTTGGCTGACCG‑ACCCAACCCAAAAAACC‑GTTTGTTGTGGGAGGGGCTTTTTTTTGGCGATGCCCAAATACCCCAAGAT‑AAAGAAAAAAGTTTTCCCC‑T‑G‑ATAAT‑TCCCCCCCACCGCCATACCCGGGATCTTttgt                                                                                                                         <  1:42851/143‑1 (MQ=255)                                      aatttttttGCTGCCCGACCATAACAAAAACAA‑C‑TTTGGCTGGGTGCGGGGTTGTTTTTTGGCGATGTCCGAATACCCCAAGAA‑AAAGTAATAGGTTTTCCCT‑T‑G‑ATAAA‑TCCCCGCCACCGCCATACCCGGGATCTTTTGTTGctct                                                                                                                   <  2:29412/143‑1 (MQ=255)                                      aatttttttGCTGACCGACCCAACCAAAAACAC‑C‑GTTTGCTGTGGGCGGGGCTTTTTTTTGGCGATGCCCCCATAGCCCAGGGA‑AAAGTAATAGGGTTTCCCT‑T‑A‑AAAAA‑TTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCta                                                                                                                  <  2:496242/144‑1 (MQ=255)                                      aaaTGTTTTGCTGACCGACCCTAACACAAACAC‑C‑GTTTGCTGTGTGACGGGCTTTTTTTGGGCGATGCCCGCATGCCCCCAGGA‑AAAAAAATAAGGTTTCCCC‑T‑A‑ATAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCta                                                                                                                  <  2:946122/149‑1 (MQ=255)                                           tttgGCTGGCCAACCCAAACAAAAACAC‑C‑TTTTGTTGTGTGGGGGGGTTTTTTTTGGCTATATCCGCATAGCCCCTGGT‑AAAGAAGTAGGTTTTCCCC‑T‑A‑AAAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATaa                                                                                                               <  1:327537/147‑1 (MQ=255)                                             ttGCTGACCGAACCTAACTCAAACAA‑C‑GTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCgggg                                                                                                         >  2:361310/1‑151 (MQ=255)                                                cagccataccatAACCCAAACAA‑C‑TTTTGTTGTG‑GGGGGGGTTTTTTTTGGGCTTTTCCGCAAAAGCCCCGGA‑AAAATAAAAAGGTTTTCCC‑G‑A‑AAAAATTTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTa                                                                                                       <  2:989362/139‑1 (MQ=255)                                                   cccaaccttaaccaaaaCAA‑C‑TTTTGCTGGGTGCGGGGTTGTTTTTTGGCGATGCCCGCACACCCCCAGAA‑AAAGAAAAATGTTTTCCCT‑T‑A‑ATAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtg                                                                                                   <  1:861078/137‑1 (MQ=255)                                                   aCCGAACCTAACTCAAACAA‑C‑GTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAtt                                                                                                    >  1:314413/1‑150 (MQ=255)                                                   aCCGAACCTAACTCAAACAA‑C‑GTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtg                                                                                                   >  1:710104/1‑151 (MQ=255)                                                    ccGAACCAAACCAAAACAC‑C‑TTTTGCTGTGTGGCGGGCTGTTTTTTGGCGTTGTCCGCATACCCCCAGGT‑AAAGTAATAAGTTTTCCCT‑T‑G‑AAAAT‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTa                                                                                                       <  2:756998/146‑1 (MQ=255)                                                       aaCCTAACTAAAACAG‑T‑TTTTGTGGTGTGAGGGGTTGTTTTTGGGTGTTGTCCGCTTAG‑CCCAGGTAAAATTAATAAGTTTTCCCC‑T‑G‑ATAAA‑TTCCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCggg                                                                                                          <  2:344902/140‑1 (MQ=255)                                                        aCCTAACTCAAACAA‑C‑GTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCg                                                                                               >  1:362349/1‑150 (MQ=255)                                                         ccTAACTCAAACAA‑C‑GTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGAt                                                                                             >  1:652774/1‑151 (MQ=255)                                                             aactaaaCAA‑CGGTTTGCGGTGGGACGGGCTTTTTTTTGGGGATCCCCGAAGAGCCCCAGAA‑AAAGAAAAAGTTTTTCCCC‑T‑A‑ATAAA‑TCCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGca                                                                                          <  2:158357/147‑1 (MQ=255)                                                                         cttttCCTGTGGGACGGGCTTTTTTTGGGCGATGCCCGCATACCCCCGGGT‑AAAGTAATAAGTTTTACCC‑T‑G‑ATAAT‑TCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                                <  2:581665/148‑1 (MQ=255)                                                                           tttGTTGGGTGCCGGTTTGCTTTTTGGCGATGCCCGAATACCCCCGGTA‑AAAGAAAAAAGGTTTCCCC‑T‑A‑ATAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                                <  2:441156/148‑1 (MQ=255)                                                                           tttGGTTTGTGCGGGGCTTTTTTTTGGCGATGTCCCCATCGCCCCGGTA‑AAATTAATAA‑GTTTTCCC‑T‑GTAAAAA‑TTCCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                                <  1:799398/148‑1 (MQ=255)                                                                           tttGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAt                                                                             >  2:382210/1‑151 (MQ=255)                                                                            ttGCTGTGTGAGGGGCTTTTTTTGGGCGATGCCCGCATACCCCAAGGA‑AAAGAAATAATGTTTCCCC‑T‑G‑ATAAA‑TCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAt                                                                             <  1:228101/150‑1 (MQ=255)                                                                             tGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGCAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCATACCCGGGATATTTTGTTGCTCTATAAGCGGGGTACTTGCGCGATGCACAATATTGTTGAttt                                                                           >  2:711254/1‑151 (MQ=255)                                                                                 gtgtgAGGGGCTGTTTTTTGGCGATGCCCGAATACCCCCAGGA‑AAAGAAAAAAGGTTTCCCT‑T‑G‑AAAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgtt                                                                        <  2:14805/150‑1 (MQ=255)                                                                                     gACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAggg                                                                   >  2:203160/1‑151 (MQ=255)                                                                                        gggTTGTTTTTGGGCGATGCCCGCACGCCCCCAGGA‑AAAGTAA‑AAGGTTTTCCCTT‑A‑ATAAA‑TCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgggg                                                                  <  1:523087/149‑1 (MQ=255)                                                                                        gggCTTTTTTTGGGCGATGCCCGCATAGCCCCAGAA‑AAAGTAATAAGTTTTACCC‑T‑A‑AAAAA‑TCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAggg                                                                   <  1:910743/148‑1 (MQ=255)                                                                                        gggCTGTTTTTGGGCTATGTCCGCATAGCCCCAGGT‑AAAGTAAAAAGTTTTCCCC‑T‑A‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgca                                                                <  2:237404/151‑1 (MQ=255)                                                                                        gggCTGTTTTTGGGCGATGCCCGCATGCCCCCAGGA‑AAAGTAATAGGTTTTCCCC‑T‑G‑ATAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                                 <  2:163075/150‑1 (MQ=255)                                                                                        gggCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgg                                                                    >  1:523907/1‑147 (MQ=255)                                                                                        gggCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                                 >  2:36366/1‑150 (MQ=255)                                                                                          gctgctTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagc                                                              >  1:931579/1‑151 (MQ=255)                                                                                                        catgcctaccaaccccAGGT‑AAAAATAAATGGTTTTCCC‑TAG‑ATAAA‑GCCCCACCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCaa                                                 <  1:610196/139‑1 (MQ=255)                                                                                                         tgcccGCATAGCCCCAGGA‑AAAGTAATAAGTTTTCCCC‑T‑G‑ATAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                  <  2:524677/145‑1 (MQ=255)                                                                                                            ccGCATAGCCCCAGGA‑AAAGTAATAAGTTTTCCCC‑T‑G‑ATAAA‑TCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAagag                                              <  2:418851/149‑1 (MQ=255)                                                                                                            ccGCATACCCCCAGGA‑AAAGTAA‑AAGGTTTTCCC‑G‑GAATAAA‑TTACCCCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAagag                                              <  1:947936/149‑1 (MQ=255)                                                                                                                aTAGCCCCAGTA‑AAAGTAATAGGTTTTCCCA‑T‑A‑AAAAT‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGta                                        <  1:791116/151‑1 (MQ=255)                                                                                                                aTAGCCCAAGAA‑AAAAAAATAGGGTTTACCC‑T‑A‑ATAAA‑TTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAagag                                              <  1:272887/145‑1 (MQ=255)                                                                                                                          gT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                              >  2:669594/1‑151 (MQ=255)                                                                                                                             aaaGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg                            >  2:981903/1‑151 (MQ=255)                                                                                                                             aaaGAAATAGGTTTTCCCC‑T‑G‑AAAAA‑TCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                              <  1:401304/149‑1 (MQ=255)                                                                                                                               aGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg                            >  1:549261/1‑149 (MQ=255)                                                                                                                               aGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg                            >  1:279014/1‑149 (MQ=255)                                                                                                                                gTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                           >  1:990352/1‑149 (MQ=255)                                                                                                                                      aGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtc                   >  2:664958/1‑151 (MQ=255)                                                                                                                                       gttttCCCC‑T‑A‑ATAAT‑TTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtc                   <  2:988883/148‑1 (MQ=255)                                                                                                                                             ccc‑T‑G‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg             >  1:928911/1‑150 (MQ=255)                                                                                                                                             ccc‑T‑A‑ATAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtttt              <  1:611604/149‑1 (MQ=255)                                                                                                                                                   g‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACa         >  1:284448/1‑150 (MQ=255)                                                                                                                                                   g‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACa         >  2:29440/1‑150 (MQ=255)                                                                                                                                                       aaa‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAgg     >  1:1022368/1‑151 (MQ=255)                                                                                                                                                           ttACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTc  >  1:876859/1‑151 (MQ=255)                                                                                                                                                             |                                                                                                                                                     CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑C‑GTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCC‑T‑G‑ATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTC  >  NC_000913/4378653‑4378950

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