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

CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGT  >  NC_000913/4378653‑4378949                                                                                                                                                          |                                                                                                                                                     ccccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ac                                                                                                                                                     >  2:229422/1‑151 (MQ=255)   ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑a                                                                                                                                                      >  2:112038/1‑148 (MQ=255)      cGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑accg‑c                                                                                                                                                 <  1:374168/149‑1 (MQ=255)       gAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑accg‑ccac                                                                                                                                              >  2:542777/1‑151 (MQ=255)        aGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑accg‑ccacc                                                                                                                                             >  1:51123/1‑151 (MQ=255)         gCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑accg‑ccaccg                                                                                                                                            <  2:832306/151‑1 (MQ=255)             tAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCAt                                                                                                                                        >  1:931025/1‑151 (MQ=255)             tAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCAt                                                                                                                                        >  2:67722/1‑151 (MQ=255)                cGTCGCCACCAAGCACGCCAGTAAATGGTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATAcc                                                                                                                                     >  1:29773/1‑151 (MQ=255)                 gTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATAccc                                                                                                                                    >  1:825494/1‑151 (MQ=255)                     ccaccaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCggg                                                                                                                                 >  2:194166/1‑150 (MQ=255)                          aaGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCttt                                                                                                                           >  2:950931/1‑151 (MQ=255)                              aCGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTttgtt                                                                                                                       >  2:243657/1‑151 (MQ=255)                               cGCAAGTAATTGTTTTGCGGCCCGAACCTAACTAAAACAACTTTTGCTGTGTGACGGGCTGTTTTTGGGCGATGTCCGCATACCCCCAGTAAAAGT‑AATAGGGTTTCCCCGA‑ATAAA‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTttgt                                                                                                                        <  2:761085/149‑1 (MQ=255)                                   aGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGctct                                                                                                                  >  1:599407/1‑151 (MQ=255)                                      aattGTTTTGCCCACCGAACAAAACTCAAACAATTTCTGTTGTGTGTCGGGCTTTTTTTTGGCGCTGTCCGCACCGCCCCTGCTAAACT‑AATATGTTTTACCTTA‑ATATT‑TC‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGctct                                                                                                                  <  1:68300/145‑1 (MQ=255)                                       aaTGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCtata                                                                                                               >  2:922166/1‑150 (MQ=255)                                       aaTGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTAACCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCtata                                                                                                               >  1:288854/1‑150 (MQ=255)                                         ttttttGCTGCCCAACCCAACCTAAAACAACGTTTGTTGTGTGCCGGGCTGTTTTTTGGCAATGCCCCCTTACCCCCGGGAAAAGT‑AAAAAGTTTTCCCTTG‑AAAAT‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAg                                                                                                             <  2:920363/148‑1 (MQ=255)                                         ttttttGCTGCCCAAACCAAACTAAAACAACGTTGGCTGTGTGGGGGGTTTTTTTTGGGCAATGCCCGCATACCCCAAGGAAAATA‑AATAATTTTTCCCCTG‑ATAAA‑TC‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAg                                                                                                             <  1:711024/148‑1 (MQ=255)                                          gTTTTCCGGCCCGAACCTAACACAAACAAGTTTTTTTGTGAGAGGGGTTTTTTTTGGGCGATTCCCGCAGCCCCCCAGGAAAAGA‑AAAAAGTTTTCCCTTA‑AAAAA‑TC‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGc                                                                                                            <  1:531232/150‑1 (MQ=255)                                           ttttGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCgg                                                                                                          >  2:570609/1‑151 (MQ=255)                                            tttGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCgtg                                                                                                         >  2:347573/1‑149 (MQ=255)                                                cggcccGACCCTAACTCAAACAACTTTTGCGGGGGGAGGGGTTTTTTTTGGGCGATGCCCGAATACCCCAAGGTAAATT‑AAAAAGTTTTCCCCTA‑AAAAA‑TT‑CCCG‑CCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGt                                                                                                       <  2:178996/145‑1 (MQ=255)                                                  gACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAtt                                                                                                   >  1:1024696/1‑151 (MQ=255)                                                   aCCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtg                                                                                                  >  2:331539/1‑151 (MQ=255)                                                    ccaaaCCAAACAAAAACAACGTTTGCTGTGGGCGGGGCTGTTTTTTGGGGATTCCCGCAGACCCCCAGGAAAAGT‑AAT‑AGGTTTCCCCTT‑AAAAATTT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtg                                                                                                  <  1:895756/147‑1 (MQ=255)                                                    ccGACCCTAACTCAAACACCGTTTGCTGGGTGACGGGCTGTTTTTTGGCGATGCCCGCATAGCCCCGGGTAAATT‑AATAAGGTTTCCCTTG‑ATAAA‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTaa                                                                                                     <  2:925213/147‑1 (MQ=255)                                                     cGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgtg                                                                                                >  2:317989/1‑151 (MQ=255)                                                      gAACCTAACCAAAACAACTTTTGCTGTGTGGGGGGCTGTTTTTTGGCGATGTCCGCATAGCCCCAGGAAAAGA‑AATAAGGTTTCCCCTA‑ATAAA‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAtt                                                                                                   <  2:134412/147‑1 (MQ=255)                                                       acccTAACAAAAACATTTTTTGTTGTGTGACGGGTTTTTTTTTGGCGTTGTCCGCATGGCCCAAGTAAAAGT‑AATAAGTTTTCCCTTA‑ATAAT‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgtg                                                                                                <  1:956030/147‑1 (MQ=255)                                                         cttAATTCAAACATCGTTTGTGGTGTGCCGGTCTGTTTTTGGGCGTTTTCCGCATACCCCCAGGTAATTA‑AATAA‑GTTTTCCCTGAAAAAA‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAtt                                                                                                   <  2:214649/142‑1 (MQ=255)                                                            aaccaaaaCAACTTTTTCTGTGTGGCGGGTTGTTTTTTGGCGATCCCCGCATAGCCCCAGGAAAAGT‑AATAAGTTTTCCCCTA‑AAAAA‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGa                                                                                             <  1:819394/142‑1 (MQ=255)                                                            aacaaaaaCACCTTTTGGTGTGTGACGGTTTGTTTTTGGGCGATGCCCGAATCCCCCCAGAAAAAGT‑AATAGGTTTTCCCTTA‑AAAAA‑TT‑CCCG‑CAACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAt                                                                                                    <  1:890416/135‑1 (MQ=255)                                                            aaCTAAACCAACGTTTGTTGGGGGAGGGGTTTTTTTTTGGCGATGCCCGCATACCCCCAGGAAAAGA‑AAAAAGGTTT‑CCCTT‑AAAAA‑TTCCCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                          <  1:243630/150‑1 (MQ=255)                                                                         tttGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAt                                                                            >  1:35474/1‑151 (MQ=255)                                                                         tttGCTGTGGGGGGGGTTGTTTTTTGGCGATGCCCGCATACCCCCAGAAAAATA‑AAAAAGGTTTCCCCTA‑ATAAA‑TT‑CCCC‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                              <  1:364426/149‑1 (MQ=255)                                                                            gCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttg                                                                         >  2:24534/1‑151 (MQ=255)                                                                             cTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgt                                                                        >  2:180242/1‑151 (MQ=255)                                                                                tgtgCCGGGTTGTTTTTTGGCGATGTCCGCATACCCCAAGGAAAAGT‑AATAAGTTTTCCCTTG‑ATAAA‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgtt                                                                       <  2:3345/149‑1 (MQ=255)                                                                                 gtaaggggCTGTTTTTTGGCATTGTCCGCATA‑CCCCCGGAAAAAT‑TAAAAAGTTTTCCCTGAATAAT‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgtt                                                                       <  1:223864/143‑1 (MQ=255)                                                                                  tgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgg                                                                   >  2:410049/1‑151 (MQ=255)                                                                                    aCGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgggg                                                                 >  1:427427/1‑151 (MQ=255)                                                                                    aCGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgggg                                                                 >  2:699154/1‑151 (MQ=255)                                                                                          tgctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagca                                                            >  2:413308/1‑150 (MQ=255)                                                                                             tttttGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtg                                                        >  2:734590/1‑151 (MQ=255)                                                                                             tttttGGCGATGCCCGCATACCCCCAGGAAAATA‑AAAAAGTTTTCCCCTA‑ATAAA‑TT‑CCCC‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCaa                                                           <  2:931133/148‑1 (MQ=255)                                                                                                  ggCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                     >  1:81169/1‑149 (MQ=255)                                                                                                    cGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGc                                                  >  2:261698/1‑150 (MQ=255)                                                                                                     gATGCCCGCAAACCCCCAGGAAAAGA‑AATAAGTTTTCCCCTA‑AAAAT‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGc                                                  <  1:229397/149‑1 (MQ=255)                                                                                                         tCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAgg                                            >  1:619912/1‑151 (MQ=255)                                                                                                         cccccATCGCCCCAGAAAAATT‑AATAGGGTTTACCCTA‑AAAAA‑TT‑CCCC‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAaga                                              <  1:842236/145‑1 (MQ=255)                                                                                                           ccgaataccccAAGTAAAATAAAAAAGTTTTCCCCTA‑AAAAT‑TT‑CCCC‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAgg                                            <  2:502759/143‑1 (MQ=255)                                                                                                           cGCATACCCCGAGAAAAAGA‑AATAAGGTTTCCCCTA‑ATAAT‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCg                                          <  1:389156/151‑1 (MQ=255)                                                                                                            caaTAGCCCTGGGTAAATA‑AAGATGTTTTCCCTTA‑ATATA‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCgg                                         <  2:81817/149‑1 (MQ=255)                                                                                                                  cccggggTAAAGT‑AA‑AAGGTTTCCCCTT‑AATAA‑TT‑ACCGCCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatta                                   <  1:180222/146‑1 (MQ=255)                                                                                                                   cccAGAAAAAAT‑AATAAGTTTTCCCCTA‑ATAAA‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaa                                  <  2:1008375/151‑1 (MQ=255)                                                                                                                        gTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                             >  2:30440/1‑151 (MQ=255)                                                                                                                        gTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                             >  1:651975/1‑151 (MQ=255)                                                                                                                          aaaGT‑AAAAAGTTTTCCCCTG‑ATAAA‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg                           <  2:489542/151‑1 (MQ=255)                                                                                                                             gT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                          >  1:978878/1‑149 (MQ=255)                                                                                                                                aataaGTTTTCCCTTA‑AAAAA‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCa                         <  2:617784/148‑1 (MQ=255)                                                                                                                                 ataaGTTTTCCCTTA‑ATAAA‑TT‑CCCC‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGc                     <  2:195869/151‑1 (MQ=255)                                                                                                                                 ataaGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGc                     >  2:735115/1‑151 (MQ=255)                                                                                                                                 aaaaGGTTTCCCCTA‑AAAAA‑TC‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGc                     <  2:461932/149‑1 (MQ=255)                                                                                                                                       tttCCCCTA‑AAAAA‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtcg                 <  2:374214/149‑1 (MQ=255)                                                                                                                                         tACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtttt             >  2:953695/1‑151 (MQ=255)                                                                                                                                                   aaa‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt       <  2:595555/148‑1 (MQ=255)                                                                                                                                                   aaa‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAgg    <  2:944234/151‑1 (MQ=255)                                                                                                                                                   aaa‑TT‑CCCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAg     <  2:205005/150‑1 (MQ=255)                                                                                                                                                   aaa‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAgg    >  2:470573/1‑151 (MQ=255)                                                                                                                                                     a‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGt  >  2:431802/1‑151 (MQ=255)                                                                                                                                                          |                                                                                                                                                     CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCCTG‑ATAAA‑TT‑ACCG‑CCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGT  >  NC_000913/4378653‑4378949

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: ‑