Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,378,802	0	A	C	45.6%	91.0 / 10.6	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/26);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 3.11e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTG  >  NC_000913/4378655‑4378952                                                                                                                                                       |                                                                                                                                                         ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTACCC‑TGATAAATT‑a‑cc‑g                                                                                                                                                     <  1:953541/151‑1 (MQ=255)    cGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTACCC‑TGATAAATT‑g‑cc‑gcc                                                                                                                                                   <  1:913470/150‑1 (MQ=255)     caGCACTACCGCCGCCACCAACCCCCCCCGCCACTGCTTTGCTGACCGCCCCTCCCTCCACCACCGTTTGCCGCGTGACGGGCTGCTTTTTGGCGATGCCCGCATACCCCCAGG‑TAAAGTAATAAGGCTTACCC‑TGATAAATT‑a‑cc‑gccac                                                                                                                                                 >  1:544077/2‑151 (MQ=255)         aaTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTACCC‑TGATAAATT‑a‑cc‑gccaccgcc                                                                                                                                             >  2:808017/1‑151 (MQ=255)                 cGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCgg                                                                                                                                     >  1:143659/1‑151 (MQ=255)                     accaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATc                                                                                                                                 >  1:188162/1‑151 (MQ=255)                        aaccACGCCAGAAAATTTTTTGCTGACCGACCCTAACTCAAACAACGTTTGTTGTGTGACGGGCTGTTTTTGGGCGATGTCCGCATAGCCCCAGG‑AAAAGTAATAAGGTTTCCCC‑TAAAAAATT‑C‑CC‑GCCACCGCCATACCCGGGATCtt                                                                                                                               <  2:275206/147‑1 (MQ=255)                        aaGCCCGCCAGTAAATTTTTTGCGGCCCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGTTTTTTGGCGATGTCCGCATACCCCCAGG‑TAAAGTAATAAGTTTTCCCC‑TGATAAATT‑C‑CC‑GCCACCGCCATACCCGGGATCttt                                                                                                                              <  1:704699/151‑1 (MQ=255)                           cccGCAAGAAATTGTTTTGCTGACCGACCCAAACACAAACAACTTTTGGGTGGGGAGGGGTTTTTTTTTGGCATTGTCCGCATCCCCCCGGG‑TAAAGTAATAAGTTTTCCCC‑TGA‑AAATT‑ACCC‑CCCACCGCCATACCCGGGATCtttt                                                                                                                             <  1:899583/147‑1 (MQ=255)                           cACGCCAGTAAATGTTTTGCTGCCCAACCCTAACTCAAACAACTTTTGCTGTGTGAGGGGCTGTTTTTTGGCGATTTCCGCATACCCCCAGG‑TAAAGTAATAGGGTTTCCCC‑TAAAAAATT‑C‑CC‑GCCACCGCCATACCCGGGATCTTttgt                                                                                                                           <  2:641555/151‑1 (MQ=255)                            aCGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTttgtt                                                                                                                          >  2:319813/1‑151 (MQ=255)                             ccccgttaatttttttGCGGCCCAAACCTAACTAAAACAACTTTTGCTGTGTGCGGGGCTGTTTTTGGGCTTTGCCCCAACACCCCCAGG‑TAAGTAAAAAATGTTTCCCT‑TGAAAATTT‑C‑CC‑GCCACCGCCATACCCGGGATCTTttgtt                                                                                                                          <  1:879155/138‑1 (MQ=255)                             cGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTttgtt                                                                                                                          >  2:250984/1‑150 (MQ=255)                                       ttttttGCTGCCCGACCCTAACTAAAACAACGTTTGCTGTGTGGCGGGCTTTTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAATTAATAAGTTTTCCCT‑TGAAAATTT‑C‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAg                                                                                                                <  2:881205/148‑1 (MQ=255)                                              cTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTa                                                                                                         >  2:619014/1‑150 (MQ=255)                                                  cgaacccTAACCAAAACAACGTTTGTTGTGGGGCGGGTTTTTTTTTGCGGTGGTCCGAATGCCCCCAGG‑AAAATAAATAAGTTTTACCC‑‑GAAAAATT‑T‑CCGGCAACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtg                                                                                                     <  1:620558/145‑1 (MQ=255)                                                  ccGAACCTACCTCAAACACCTTTTGTTGTGTGACGGGCTGTTTTTTGGCGATGCCCGCATACCCCAAGG‑AAAAGTAAAAAGTTTTCCCC‑TGAAAAATT‑C‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgt                                                                                                    <  2:612102/151‑1 (MQ=255)                                                  ccGAACCTAACTAAAACACCGTTTGCTGTGTGACGGGTTGTTTTTTGGCGATGTCCGCATAGCCCCAGG‑AAAAGTAATAAGTTTTCCCC‑TGATAAATT‑C‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgt                                                                                                    <  2:275120/151‑1 (MQ=255)                                                  ccGAACCTAACTAAAACAACGTTTGCGGTGTGACGGGCTGTTTTTGGGCGATGCCCGAATAGCCCCAGG‑AAAAGAAATAAGGTTTCCCC‑TAATAATTT‑C‑CC‑CCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAACtgt                                                                                                    <  2:612129/151‑1 (MQ=255)                                                                 acaccGTTTGCTGTGTGCGGGGCTGCTTTTTGGCGGTGTCCGCAGGGCCCAAGG‑AAAAGTAATGGGGTTAGCCG‑AGATAAATT‑C‑CC‑GCCACCGCCATACCCGGGATCTTTTATTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATc                                                                                     <  2:51987/151‑1 (MQ=255)                                                                 acaacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATc                                                                                     >  1:271919/1‑151 (MQ=255)                                                                  caacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCg                                                                                    >  1:225993/1‑151 (MQ=255)                                                                   accttttGCTGTGTGAGGGGCTGTTTTTTGGGGATGTCCGCATACCCCCAGG‑TAAAGAAATAAGTTTTACCC‑TAATAAATT‑C‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGt                                                                                   <  2:357579/147‑1 (MQ=255)                                                                    acGTTTGTTGTGTGGGGGGCTGTTTTTTGGCGATTCCCCCATACCCCAAGG‑TAAAGTAATAAGGTTTCCCT‑TAATAAATT‑C‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                                  <  2:325051/151‑1 (MQ=255)                                                                         tGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAttt                                                                             >  1:848261/1‑151 (MQ=255)                                                                             gggtggggggTTGTTTTTTGGCTTTTCCCGCAAACCCCCAGG‑TAAA‑TAAAAAGGTTTTCCCTTGAAAATTT‑C‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgtt                                                                          <  1:1013698/143‑1 (MQ=255)                                                                              tgtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                         >  2:780743/1‑150 (MQ=255)                                                                                 gACGGGTTGTTTTTTGGCGATGTCCGCATACCCCCAGA‑AAAAGTAATAAGTTTTCCCC‑TAATAATTT‑C‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAggg                                                                     <  2:291839/151‑1 (MQ=255)                                                                                  aCGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgggg                                                                    >  1:443286/1‑151 (MQ=255)                                                                                    gggTTGCTTTTGGGCGTGGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTACCC‑TGATAATTT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                                 <  1:811073/136‑1 (MQ=255)                                                                                     ggCTGTTTTTGGGGGATGTCCGTATACCCCCAGG‑AAAAGAAAAAATGTTT‑CCC‑TTAAAAATTAA‑CC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcag                                                                 <  2:366346/151‑1 (MQ=255)                                                                                     ggCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                                   >  1:24694/1‑149 (MQ=255)                                                                                       ctgctTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagc                                                                >  1:312986/1‑150 (MQ=255)                                                                                        tgctTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagca                                                               >  2:792981/1‑150 (MQ=255)                                                                                            ttttGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtg                                                           >  1:833990/1‑150 (MQ=255)                                                                                                ggCGTTGTCCGCAT‑GCCCCAGGAAAAATTAATAGGTTTTCCCT‑TAATAAATT‑C‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                        <  1:182945/149‑1 (MQ=255)                                                                                                        ccGCATAGCCCCAGG‑AAAAGTAATAAGTTTTACCC‑TAAAAAATT‑C‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAagag                                                <  2:1014855/149‑1 (MQ=255)                                                                                                        ccGCATACCCCAAGG‑TAAAGAAATAAGTTTTCCCC‑TGATAATTT‑C‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAaga                                                 <  2:517204/148‑1 (MQ=255)                                                                                                        ccGCATA‑CCCCAAG‑GAAAATAAAAAGGTTTTCCCTTAAAAAATC‑C‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAgg                                               <  1:85945/150‑1 (MQ=255)                                                                                                                ccccAGA‑AAAAAAAATAAGTTTTCCCT‑TGAAAATTT‑C‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatta                                      <  1:967878/151‑1 (MQ=255)                                                                                                                     gg‑TAAAGTAATAAGTTTTCCCC‑TAAAAAATT‑C‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaataa                                  <  2:373798/150‑1 (MQ=255)                                                                                                                     gg‑TAAAGTAATAAGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaata                                   >  2:785347/1‑149 (MQ=255)                                                                                                                      g‑TAAAGTAATAAGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                                >  2:695614/1‑151 (MQ=255)                                                                                                                        tAAAGTAATAAGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                                >  2:859716/1‑150 (MQ=255)                                                                                                                           aGTAATAAGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                             >  2:845729/1‑150 (MQ=255)                                                                                                                            gTAATAAGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                             >  2:570868/1‑149 (MQ=255)                                                                                                                             taataaGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                             >  2:430181/1‑148 (MQ=255)                                                                                                                             taataaGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                             >  1:184063/1‑148 (MQ=255)                                                                                                                              aataaGGTTTCCTC‑TAATAATTT‑C‑CC‑GCCGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCaga                          <  1:304558/150‑1 (MQ=255)                                                                                                                              aataaGGTTTCCCT‑TAATAAATT‑C‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCa                            <  2:417733/148‑1 (MQ=255)                                                                                                                              aaaaaGTTTTCCCC‑TAAAAATTT‑A‑CC‑CCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCaga                          <  2:839598/147‑1 (MQ=255)                                                                                                                               ataaGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGc                        >  2:44245/1‑151 (MQ=255)                                                                                                                                taaGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcg                       >  1:35132/1‑151 (MQ=255)                                                                                                                                 aaGTTTTCCCC‑TAAAAAATT‑C‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcg                       <  2:381233/150‑1 (MQ=255)                                                                                                                                 aaGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgt                      >  2:850946/1‑151 (MQ=255)                                                                                                                                   ggTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtcg                    >  2:55305/1‑151 (MQ=255)                                                                                                                                      ttcccc‑TGATAATTT‑C‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtt                  <  2:386208/147‑1 (MQ=255)                                                                                                                                         cccttaaaaaaTT‑C‑CC‑CCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg               <  1:53183/143‑1 (MQ=255)                                                                                                                                          ccttGAAAATTT‑C‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtttt                <  2:603536/146‑1 (MQ=255)                                                                                                                                               aTAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACa           >  1:227843/1‑149 (MQ=255)                                                                                                                                               aTAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt          >  1:5445/1‑150 (MQ=255)                                                                                                                                                 aaaTT‑C‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGaa             <  1:818985/145‑1 (MQ=255)                                                                                                                                                       a‑cc‑gccaccgccATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTg  >  1:738446/1‑151 (MQ=255)                                                                                                                                                       a‑cc‑gccaccgccATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTg  >  2:984567/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                         CCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTACCC‑TGATAAATT‑A‑CC‑GCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTG  >  NC_000913/4378655‑4378952

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