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

CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGG  >  NC_000913/4378653‑4378948                                                                                                                                                           |                                                                                                                                                   ccccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ac                                                                                                                                                   >  1:603032/1‑151 (MQ=255) ccccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ac                                                                                                                                                   >  2:1137482/1‑151 (MQ=255)   ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑accg                                                                                                                                                 >  1:391685/1‑151 (MQ=255)      cGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCAATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑accgcca                                                                                                                                              <  2:563183/151‑1 (MQ=255)         gCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑accgccacc                                                                                                                                            <  1:447306/150‑1 (MQ=255)          cAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑accgccaccgc                                                                                                                                          >  2:659323/1‑151 (MQ=255)          cAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑accgccaccgc                                                                                                                                          >  2:1198359/1‑151 (MQ=255)             tAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCAt                                                                                                                                       >  1:437903/1‑151 (MQ=255)              aGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCAt                                                                                                                                       >  1:360285/1‑150 (MQ=255)                 gTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATAccc                                                                                                                                   >  2:839133/1‑151 (MQ=255)                             cACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTttgt                                                                                                                       >  1:27647/1‑151 (MQ=255)                             cACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTttgt                                                                                                                       >  1:24984/1‑151 (MQ=255)                                gCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGc                                                                                                                    >  1:477983/1‑151 (MQ=255)                                 ccAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGct                                                                                                                   >  1:233713/1‑151 (MQ=255)                                  aataaatttttttGCTGCCAAACCCAAACAAAAACAACTTTTGTTGTGGGGGGGGTTGTTTTTTGGCGTTGCCCGAATACCCCAAGAAAAAAAAAAAAGTTTT‑CCCT‑TA‑AAAAA‑TT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGctc                                                                                                                  <  2:464891/142‑1 (MQ=255)                                   aGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGctct                                                                                                                 >  2:927938/1‑151 (MQ=255)                                     tAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                               >  1:1115459/1‑151 (MQ=255)                                      aatttttttGCTGCCCAAACCAACCTCAAAAAGCTTTTGTTGGTGGAGGGGCTGTTTTTGGGCGTTCCCCCCATCCCCCCAGGAAAAGAAATAATTTTT‑CCCC‑TA‑AAAAA‑TC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtata                                                                                                              <  2:106006/146‑1 (MQ=255)                                      aatttttttGCTGACCAACCCTAACTAAAAAAATGTTTGTTGTGTGAGGGGCTGTTTTTGGGCGTTGCCCGCATACCCCCAGGAAAAGTAATAGGTTTT‑CCCT‑TA‑AAAAA‑TT‑CCCCCCCCCGCCATACCCGGGATCTTTTGTTGctc                                                                                                                  <  2:182537/142‑1 (MQ=255)                                        aTGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtata                                                                                                              >  1:19322/1‑149 (MQ=255)                                               gCTGACGAAACCAAACAAAAACACCGTTTGCTGGGTGAGGGGCTGTTTTTTGGCGATGCCCGCATGGCCCCAGAAAAAGTAAAAAGTTTT‑CCCC‑TG‑AAAAA‑TC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCggg                                                                                                        <  2:616074/148‑1 (MQ=255)                                               gCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGt                                                                                                      >  1:926002/1‑150 (MQ=255)                                                   aCCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtg                                                                                                 >  2:970809/1‑151 (MQ=255)                                                   aCCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtg                                                                                                 >  1:1203088/1‑151 (MQ=255)                                                    ccgaccttaccaaaaaCAACTTTTGTGGGGTGAGGGGTTTTTTTTTGGGGTTCCCCGC‑TACCCCCAAGTAAAAGAAATAGGTTT‑TCCC‑TT‑AAAAA‑TTCCCCGCAACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtg                                                                                                 <  1:174037/137‑1 (MQ=255)                                                    ccgacccaacccaaaaCACCGTTTGCTGTGTGACGGGCTGTTTTTGGGCGAGGCCCGAATACCCCCGGGAAAATTAATAGGTTTT‑CCCC‑TA‑AAAAT‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgt                                                                                                <  2:908089/138‑1 (MQ=255)                                                       aaccaaaccaAAAAAACGTTGGTGGTGTGCGGGGCTGTTTTTGGGCGTTGCCCGAATA‑CCCCAAGTAAAATAAAAAGGTTT‑TCCCTTA‑ATAAT‑TC‑CCCGCAACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgt                                                                                                <  1:355832/148‑1 (MQ=255)                                                            accaaaaaCAATGTTTGGTGTGCGCGGGTCTTTTTTTGGGAGTTGCCCGCATACCCCAAGAAAAAGTAATAAGTTTT‑TCCG‑AA‑AAAAT‑TC‑CCCCCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGc                                                                                              <  2:391686/140‑1 (MQ=255)                                                            aaCTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                         >  2:20133/1‑150 (MQ=255)                                                              cTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac                                                                                       >  1:1110416/1‑150 (MQ=255)                                                                 aaaCATCTTTTGCTGTGTGCCGGGTTGTTTTTGGGCGATCCCCCCATACCCCCGGTAAAAGTAAAAAGTTTT‑CCCC‑TA‑AAAAA‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAata                                                                                   <  2:562873/151‑1 (MQ=255)                                                                 aaaCAGCTTTTGCTGTGGGACGGGCTGTTTTTGGGCGATGCCCGCATAGCCCAAGGTAAAGAAATAGGTTTT‑CCCG‑TA‑AAAAA‑TC‑CCCCCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                         <  2:66874/145‑1 (MQ=255)                                                                 aaaCAACTTTTGCTGTGTGAGGGGTTGTTTTTTGGCGATGTCCGCATACCCCCAGGAAAAGAAATAGGTTTT‑CCCC‑TA‑AAAAA‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac                                                                                       <  2:1145890/147‑1 (MQ=255)                                                                     accttttGTTGTGTGCGGGGTTTTTTTTGGGCGTTGCCCGAAT‑CCCCCAAGAAAAAAAAATAGGTTT‑TCCC‑TT‑AAAAATTT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGt                                                                               <  2:12338/147‑1 (MQ=255)                                                                      acttttGTTGTGTGCGGGGTTGTTTTTGGGCGTTGCCCGAATACCCCCGGGAAAAGTAAAAGGTTTT‑CCCC‑TA‑AAAAA‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                              <  2:1225973/148‑1 (MQ=255)                                                                       cGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                             >  2:1180711/1‑151 (MQ=255)                                                                         tttGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAt                                                                           >  2:331451/1‑151 (MQ=255)                                                                            gCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttg                                                                        >  2:692117/1‑151 (MQ=255)                                                                             cTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgt                                                                       >  2:331424/1‑151 (MQ=255)                                                                              tgtgtgAGGGGCTTTTTTTGGGCGATCCCCGCATACCCCGAGATAAAGTAA‑AAGTTTTCCCCC‑TA‑ATAAA‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAttt                                                                         <  1:382156/148‑1 (MQ=255)                                                                               gtgggcggggCTCTTTTTGGGCAATCCCCGCATCGCCCCAGGAAAATAAATATGTTTT‑CCCC‑TA‑AAAAT‑TC‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttg                                                                        <  1:652381/141‑1 (MQ=255)                                                                                tgtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                     >  2:1133151/1‑150 (MQ=255)                                                                                   gaggggCTGTTTTTGGGGGTTTTCCGCAAACCCCGGGGAAAAGAAAGAGGTTTT‑CCCT‑TG‑AAAAT‑TC‑CCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAggg                                                                 <  2:620883/148‑1 (MQ=255)                                                                                   gACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAggg                                                                 >  1:620881/1‑151 (MQ=255)                                                                                     cGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                               >  1:559651/1‑151 (MQ=255)                                                                                                  ggCGATGCCCCCATCGCCCCA‑GAAAAATAAAAAAGTTT‑TCCC‑TT‑AAAAATTT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGa                                                   <  1:1226428/150‑1 (MQ=255)                                                                                                     gATGCCCCCATACCCCAAGAAAAAAAAATATGTTTT‑CCCT‑TA‑AAAAT‑TT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                <  1:577261/150‑1 (MQ=255)                                                                                                         cccacatccccccAGTAAAAATAAAAAGTTTT‑CCCC‑TA‑AAAAA‑TC‑CCCCCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAaga                                             <  2:873683/140‑1 (MQ=255)                                                                                                          ccGCATGCCCCCAGGAAAAGTAATAAGTTTT‑CCCT‑TA‑ATAAT‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAagag                                            <  2:130410/149‑1 (MQ=255)                                                                                                          ccGCATACCCCGAGATAAAATAATAAGTTTT‑CCCC‑TA‑AAAAA‑TC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAag                                              <  2:14002/147‑1 (MQ=255)                                                                                                              aTAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                       >  1:1015506/1‑150 (MQ=255)                                                                                                              aTACCCCAAGGAAAAGAAAAATGTTTT‑CCCC‑TG‑AAAAT‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGta                                      <  1:1214813/151‑1 (MQ=255)                                                                                                                  cccgAGGAAAAGAAATAAGTTTT‑ACCC‑TA‑ATAAT‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatt                                   <  2:739467/150‑1 (MQ=255)                                                                                                                    ccAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaat                                >  2:676305/1‑151 (MQ=255)                                                                                                                         aaaaGAAATAGGGTTT‑ACCT‑TA‑AAAAT‑TC‑ACCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAg                             <  1:467589/148‑1 (MQ=255)                                                                                                                          aaaGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg                          >  1:530025/1‑151 (MQ=255)                                                                                                                           aaagaaatatgTTT‑TCCC‑TT‑AAAAA‑TTCCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgc                           <  1:1229434/140‑1 (MQ=255)                                                                                                                           aaGTAA‑AATGTTT‑TCCC‑TGAATAAC‑TC‑ACGCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg                          <  2:127692/150‑1 (MQ=255)                                                                                                                           aaGAAAAAAGGTTT‑CCCC‑TA‑ATAAT‑TC‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgc                           <  1:812989/149‑1 (MQ=255)                                                                                                                               aatagGGTTT‑CCCT‑TA‑ATAAT‑TC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCagag                     <  1:373812/151‑1 (MQ=255)                                                                                                                               aataaGGTTT‑CCCC‑TG‑ATAAA‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCagag                     <  2:1120718/151‑1 (MQ=255)                                                                                                                               aataaGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCag                       >  1:739452/1‑149 (MQ=255)                                                                                                                               aaaaggtttt‑CCCC‑TA‑AAAAT‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCag                       <  2:1053595/142‑1 (MQ=255)                                                                                                                                aagaGGTTT‑TCCCTTG‑AAAAA‑TC‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCa                        <  2:77571/145‑1 (MQ=255)                                                                                                                                    ggTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtcg                >  2:512406/1‑151 (MQ=255)                                                                                                                                       tt‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGttt             >  2:62213/1‑151 (MQ=255)                                                                                                                                       tt‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGttt             >  2:188440/1‑151 (MQ=255)                                                                                                                                            cc‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGa          >  1:894812/1‑150 (MQ=255)                                                                                                                                                  aTAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt      >  2:694704/1‑150 (MQ=255)                                                                                                                                                   tAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACa       >  2:201076/1‑148 (MQ=255)                                                                                                                                                   tAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt      >  2:751469/1‑149 (MQ=255)                                                                                                                                                    aaa‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAgg   >  2:229861/1‑151 (MQ=255)                                                                                                                                                    aaa‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAgg   >  1:447043/1‑151 (MQ=255)                                                                                                                                                     aa‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg  <  1:722084/151‑1 (MQ=255)                                                                                                                                                           |                                                                                                                                                   CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTT‑ACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGG  >  NC_000913/4378653‑4378948

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