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

CCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCT  >  NC_000913/4378654‑4378951                                                                                                                                                          |                                                                                                                                                      cccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTacc                                                                                                                                                     <  1:1215845/151‑1 (MQ=255)  ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTaccg                                                                                                                                                    >  1:321798/1‑151 (MQ=255)    tCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTaccgcc                                                                                                                                                  >  2:250634/1‑151 (MQ=255)     cGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTaccgcca                                                                                                                                                 <  1:397790/151‑1 (MQ=255)      gAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTaccgccac                                                                                                                                                >  2:964027/1‑151 (MQ=255)        gCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTaccgccaccg                                                                                                                                              >  2:402266/1‑151 (MQ=255)          aaTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTaccgccaccgcc                                                                                                                                            >  2:326788/1‑151 (MQ=255)            tAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCAt                                                                                                                                          >  2:744284/1‑151 (MQ=255)               cGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATAcc                                                                                                                                       >  2:1172908/1‑151 (MQ=255)                  cGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCgg                                                                                                                                    >  2:1156087/1‑151 (MQ=255)                   gCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTATCGCCACCGCCATACCCggg                                                                                                                                   >  2:210980/1‑151 (MQ=255)                    ccaccaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCggg                                                                                                                                   >  1:493855/1‑150 (MQ=255)                      accaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATc                                                                                                                                >  1:462689/1‑151 (MQ=255)                         aaGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCttt                                                                                                                             >  1:156304/1‑151 (MQ=255)                            cACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTttgt                                                                                                                          >  1:689659/1‑151 (MQ=255)                              cacaaGTAAGTTTTTTGCGGCCAAAACCAAACTAAAACATCTTTTGTGGTGGGAGGGGCT‑TTTTTTGGGGATTGCCCCCATACCCCCAGG‑AAAA‑GTAATAAGT‑TTTCCCTTAAAAAT‑TTCCCGCCACCGCCATACCCGGGATCTTttgttg                                                                                                                        <  2:45334/147‑1 (MQ=255)                              cGCCAGTAATTTTTTTGCTGACCGACCCTAACCAAAACAACTTTTGCTGTGTGGCGGGCTGTTTTTTGGCGA‑TGCCCGAATAGCCCCAGG‑AAAA‑GTAATAAGT‑TTTCCCCTAAAAAA‑TTCCCGCCACCGCCATACCCGGGATCTTttgt                                                                                                                          <  1:25618/149‑1 (MQ=255)                                caattaatttttttGCTGCCCGACCCAAACCCAACCACCGTTTGCTGTGGGACGGGCTTTTTTTTGGGGA‑TGCCCGCATACCCCCGGA‑AAAA‑ATAATAAGT‑TTACCCCTAAAAAA‑TTCCCGCCACCGCCATACCCGGGATCTTttgttg                                                                                                                        <  1:334982/139‑1 (MQ=255)                                  aGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGctct                                                                                                                    >  2:307524/1‑151 (MQ=255)                                    tAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                                  >  1:818626/1‑151 (MQ=255)                                    tAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                                  >  1:375546/1‑151 (MQ=255)                                     aaaTGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                                  >  2:1128332/1‑150 (MQ=255)                                        ttttttGCTGACCGACCCAAACCAAAACAACTTTTGCTGTGTGGGGGGCTGTTTTTTGGCGA‑TGCCCGCAAACCCCCAGG‑AAAA‑TAAATAAGT‑TTTCCCTTGAAAAT‑TCCCCGCCGCCGCCATACCCGGGATCTTTTGTTGCTCTATaa                                                                                                                <  2:653795/147‑1 (MQ=255)                                             tGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGt                                                                                                         >  1:252389/1‑151 (MQ=255)                                                tGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAt                                                                                                      >  2:700949/1‑151 (MQ=255)                                                tGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAt                                                                                                      >  1:663871/1‑151 (MQ=255)                                                   ccaacccaaaCTCAAACACCGTTTGCTGTGTGCGGGGCTGTTTTTGGGCGA‑TGTCCGCATACCCCCAGG‑AAAA‑GTAATAAGT‑TTTCCCCTGATAAT‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgt                                                                                                   <  1:859677/143‑1 (MQ=255)                                                           aaCTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                            >  2:861473/1‑150 (MQ=255)                                                              tCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcaca                                                                                         >  2:918446/1‑150 (MQ=255)                                                                aaaaatcttttggggtgtgAGGGGCTGTTTTTGGGGGA‑GGTCCACATAGCCCCAGG‑AAAA‑G‑AAAAAGG‑TTTTCCCTTAAAAATTTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAata                                                                                      <  2:306395/137‑1 (MQ=255)                                                                aaaaaaCGTTTGCTGTGTGCCGGGTTTTTTTTTGGGGA‑TGCCCCCATACCCCCGGG‑AAAA‑GAAATAAGT‑TTTCCCGTAAAAAA‑TCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                       <  1:812725/150‑1 (MQ=255)                                                                aaaCACCTTTTGCTGTGTGCCGGGTTTTTTTTGGGGGT‑TGCCCGCATCCCCCCAGG‑AAAA‑GTAAAAAGT‑TTTCCCTTAATAAT‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGca                                                                                           <  1:1075834/146‑1 (MQ=255)                                                                aaaCAACGTTTGCTGTGTGACGGGCTGTTTTTTGGCGA‑TGTCCGCATCCCCCCAGG‑AAAA‑GAAAAAAGT‑TTTCCCCTGATAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACaa                                                                                        <  1:868907/149‑1 (MQ=255)                                                                aaaCAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                       >  2:642660/1‑150 (MQ=255)                                                                aaaCAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                       >  2:436423/1‑150 (MQ=255)                                                                 aacaacTTTTGCTGTGTGACGGGCTGTTTTTTGGCGA‑TGTCCGAATACCCCAAGG‑AAAA‑TTAAAAGGG‑TTTCCCTTAATAAT‑TTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAata                                                                                      <  1:885701/150‑1 (MQ=255)                                                                               tgtgACGGGTTGTTTTTTGGCGA‑TGCCCGCATACCCCCGGG‑TAAA‑GTAATAGGG‑TTTACCCTGAAAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttg                                                                           <  1:77827/147‑1 (MQ=255)                                                                                gtgAGGGGCTGTTTTTTGGCGT‑TGCCCGCATA‑CCCCAGGTTAAA‑TTAATAAGT‑TTTCCCTTAATAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAttt                                                                            <  1:842072/145‑1 (MQ=255)                                                                                gtgACGGGTTTTTTTTGGGCGA‑TGCCCGCATAGCCCAAGG‑TAAA‑TTAATAAGT‑TTTCCCCTGAAAAT‑TTCCCGCCGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAg                                                                      <  1:515511/151‑1 (MQ=255)                                                                                gtgACGGGCTGTTTTTTGGCGA‑TGTCCGAATACCCCCAGG‑AAAA‑GTAATAAGG‑TTTCCCCTAATAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgtt                                                                         <  1:23282/148‑1 (MQ=255)                                                                                     gggCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                                  >  2:36839/1‑150 (MQ=255)                                                                                     gggCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                                  >  1:269765/1‑150 (MQ=255)                                                                                     ggcctgttTTTGGGCGT‑TGTCCGCTTAGCCCAAGG‑TAAA‑GAAATAAGT‑TTTCCCTTGAAAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                                  <  1:965134/147‑1 (MQ=255)                                                                                         tttttttgGGCGA‑TCCCCGCATA‑CCCCCGG‑TAAAGGTAAT‑AGT‑TTTTCCCTTAAAAATTTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagc                                                               <  2:530038/146‑1 (MQ=255)                                                                                              tttGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgtg                                                        >  2:619635/1‑151 (MQ=255)                                                                                              tttGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgtg                                                        >  2:18253/1‑151 (MQ=255)                                                                                               ttGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                       >  2:116853/1‑151 (MQ=255)                                                                                               ttGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                       >  2:138426/1‑151 (MQ=255)                                                                                                  gCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGc                                                    >  1:747492/1‑151 (MQ=255)                                                                                                    gA‑TGCCCGCATACCCCCGGG‑AAAA‑GAAATAAGG‑TTTCCCCGAAAAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCaa                                                  <  1:449133/151‑1 (MQ=255)                                                                                                    gA‑TGCCCGAATACCCCAAGA‑AAAA‑GAAATAGGT‑TTTCCCCTAAAAAT‑TTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCaa                                                  <  1:100313/151‑1 (MQ=255)                                                                                                    cgttGTCCGCATACCCCCAGG‑AAAA‑GTAATAGGT‑TTTCCCCTGATAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                   <  1:256339/148‑1 (MQ=255)                                                                                                      atccccGCATACCCCCAGG‑TAAA‑AAAAAAGGG‑TTTCCCCTAAAAAA‑TTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCaa                                                  <  2:679790/146‑1 (MQ=255)                                                                                                          ccccatccccccAGA‑AAAA‑GTAATAAGT‑TTTCCCCTAAAAAT‑TCACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAaga                                                <  2:321792/140‑1 (MQ=255)                                                                                                          ccGCATACCCCCAGG‑AAAA‑AAAAAAGGT‑TTTCCCCTAATAAT‑TTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAgg                                              <  1:198188/150‑1 (MQ=255)                                                                                                             cATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                          >  1:447222/1‑151 (MQ=255)                                                                                                              aTAGCCCGAGG‑AAAA‑ATAATAAGT‑TTTCCCTTGAAAAT‑TTCCCGCCGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                          <  1:865817/150‑1 (MQ=255)                                                                                                              aTAGCCCAAGG‑AAAA‑GTAATAAGT‑TTTCCCCTAATAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGta                                         <  2:1176913/151‑1 (MQ=255)                                                                                                               ataccccAAG‑GAAA‑ATAAAAAGTTTTTCCCCTAAAAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCg                                            <  2:696951/145‑1 (MQ=255)                                                                                                                aGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaa                                        >  1:528496/1‑150 (MQ=255)                                                                                                                  ccccAGA‑AAAA‑GTAATAGGT‑TTTACCTTGAAAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                          <  2:492555/146‑1 (MQ=255)                                                                                                                            aa‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg                             >  1:1221509/1‑150 (MQ=255)                                                                                                                               gTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                            >  1:762973/1‑149 (MQ=255)                                                                                                                                 taacacgttttCCCCTGATAAA‑TCCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCa                           <  1:397315/141‑1 (MQ=255)                                                                                                                                 aataaGT‑TTTCCCCTAAAAAT‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCa                           <  2:375538/148‑1 (MQ=255)                                                                                                                                 aataaGG‑TTTCCCATAAAAAA‑TCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCaga                         <  1:577623/150‑1 (MQ=255)                                                                                                                                   taaGT‑TTTTCCCTGAAAAA‑TTCCCGCAACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCaga                         <  2:745161/148‑1 (MQ=255)                                                                                                                                         tttccctttttAAT‑TTCCCGCCGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtcgt                  <  1:452499/140‑1 (MQ=255)                                                                                                                                         tttCCCCTGAAAAA‑TTCCCGCCGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGc                       <  2:76106/145‑1 (MQ=255)                                                                                                                                           tccccTAAAAAA‑TTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGttt                <  1:1180068/148‑1 (MQ=255)                                                                                                                                             ctatGAAAAA‑TTCCCGCCGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGttt                <  2:809217/145‑1 (MQ=255)                                                                                                                                                  aTAAA‑TCCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGa             <  1:714379/146‑1 (MQ=255)                                                                                                                                                        ttACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGAGCAGAGCGTCGGTTTGAACATAgggac   >  2:1200766/1‑149 (MQ=255)                                                                                                                                                         tACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGCCGTTTTGAACATAGGGTCt  >  1:1152651/1‑151 (MQ=255)                                                                                                                                                          acccccaccgccATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGAAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGt    >  1:308744/1‑148 (MQ=255)                                                                                                                                                          |                                                                                                                                                      CCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGA‑TGTCCGCATAGCCCCAGG‑TAAA‑GTAATAAGG‑TTTACCCTGATAAA‑TTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCT  >  NC_000913/4378654‑4378951

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