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

CCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAT  >  NC_000913/4378655‑4378944                                                                                                                                                         |                                                                                                                                               ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑accg                                                                                                                                             >  2:821589/1‑151 (MQ=255) ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AACAAGG‑TTTACCC‑TGATAAA‑TT‑accg                                                                                                                                             >  2:642036/1‑151 (MQ=255)       gCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑accgccaccg                                                                                                                                       >  1:853218/1‑151 (MQ=255)       gCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑accgccaccg                                                                                                                                       >  2:711687/1‑151 (MQ=255)        cAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑accgccaccgc                                                                                                                                      >  1:1184898/1‑151 (MQ=255)         aaTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑accgccaccgcc                                                                                                                                     >  1:359876/1‑151 (MQ=255)         aaTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGGAAAGT‑AATACGG‑TTTACCC‑TGATAAA‑TT‑accgccaccgcc                                                                                                                                     >  2:337364/1‑151 (MQ=255)          aTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑accgccaccgcc                                                                                                                                     >  2:289728/1‑150 (MQ=255)                 cGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCgg                                                                                                                             >  1:234417/1‑151 (MQ=255)                     accaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATc                                                                                                                         >  2:1187994/1‑151 (MQ=255)                       caAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCtt                                                                                                                       >  1:409419/1‑151 (MQ=255)                        aaGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCttt                                                                                                                      >  1:1141973/1‑151 (MQ=255)                             cGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGAGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTttgttg                                                                                                                 >  2:86449/1‑151 (MQ=255)                                 aGTAAATTTTTCGCTGCCCGAACCAAACCCAACCATGTTTTCGTGTGTGCCGGGCTTTTTTTGGGCGAGTCCCGCA‑GCCCCCAGGTAAAATAAAAAAGG‑TTTCCCT‑TAAAAAA‑TC‑CCCCCAACCGCCATACCCGGGATCTTTTGTTGctct                                                                                                             <  1:325860/151‑1 (MQ=255)                                  gTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGctct                                                                                                             >  1:302725/1‑150 (MQ=255)                                  gTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGctct                                                                                                             >  2:76502/1‑150 (MQ=255)                                       tttttggcggcccaaaccaaaccaAAACATCTTTTGTTGTGGGAGGGGCTGTTTTTGGGGGTTGCCCAAATACCCCCAGAAAAATT‑AATAGGT‑TTTCCCC‑TAAAAAT‑TC‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATaa                                                                                                         <  1:510157/135‑1 (MQ=255)                                             gCTGACCAAACCTACCAAAACAAAGGTTGGCGGTGTGGGGGCCTTTTTTTTGGCTATGCCCGAATAGCCCGAGAAAAAGA‑AAT‑AGG‑TTTTCCA‑CAAAAAA‑TTCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCggg                                                                                                    <  2:493147/148‑1 (MQ=255)                                              ctacccaaaCCAAACAAAAAGAACTTGTGGTGTGTGTGGGTCTGTTTTTTGGCGATCCCCCCATGGCCCCGGAAAAAAA‑AATAAGT‑TTTCCCT‑TAAAAAA‑TC‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGt                                                                                                  <  1:382196/142‑1 (MQ=255)                                              cTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTa                                                                                                 >  2:1000287/1‑150 (MQ=255)                                              cTGACCGAACCTAACTCAAACAACGTTTGCCGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGGAAAGT‑AATAAGG‑CTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTa                                                                                                 >  2:839713/1‑150 (MQ=255)                                                 accaaaccaaCCACAAACACGTTTTGGTGTGTGGGGGTCTTTTTTGGGGGGATGCCCTCATGCCCCCAGAAAAAAT‑AATAGGG‑TTTCCCT‑TGATAAT‑TT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAtt                                                                                              <  1:156487/150‑1 (MQ=255)                                                   cGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgtg                                                                                           >  2:498071/1‑151 (MQ=255)                                                      aCCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCg                                                                                         >  1:883598/1‑150 (MQ=255)                                                        cTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGA‑AATAAGG‑CTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGGTCTTTTGTTGCTCTATAAGCGGTGTAATTGTGCGAt                                                                                       >  2:208117/1‑150 (MQ=255)                                                         tAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                     >  1:471817/1‑151 (MQ=255)                                                          acccaaaccaccttttGCTGTGTGAGGGGTTTTTTTTTGGCGATCCCCGAATACCCCCAGAAAAAGA‑AATAAGT‑TTTCCCT‑TGAAAAT‑TC‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGa                                                                                        <  1:1031916/134‑1 (MQ=255)                                                          acccaaaaCATCGTTTGTTGGGTGGCGGGTTGTTTTTGGGGGTTGCCCGCATGCCCCCAGAAAAATA‑AAAAAGT‑TTTCCCC‑TAATAAA‑TT‑CCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGca                                                                                    <  1:1188071/146‑1 (MQ=255)                                                               aaaCATCGTTGGCTGTTGGTGGGGCTTTTTTTGGGCGTTCTCCGCATACCCCCTGGTAAAAT‑AAAAGGG‑TTTACCA‑TAATAAA‑TT‑CCCCCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                <  1:326944/150‑1 (MQ=255)                                                                aacatcTTTTGTTGTGTGAGGGGTTTTTTTTGGGGGATTCCCGCATACCCCAAGAAAAAGA‑AAAAAGT‑TTTCCCT‑TGAAAAT‑TT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAatat                                                                              <  1:421834/151‑1 (MQ=255)                                                                    acGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCTCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                          >  2:19086/1‑151 (MQ=255)                                                                       tttGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAt                                                                       >  1:206894/1‑151 (MQ=255)                                                                        ttgttgtggtgggggCTTTTTTTGGGCGAGGTCCACATCGCCCGGGGAAAGGT‑AA‑AAGG‑TTTTCCCTTAATAAT‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAt                                                                       <  1:659636/138‑1 (MQ=255)                                                                         tGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAttt                                                                     >  1:1183133/1‑151 (MQ=255)                                                                           cTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgt                                                                   >  1:431344/1‑151 (MQ=255)                                                                                 gaggggTTGTTTTTTGGCGTTGCCCCCAT‑CCCCCAAGTAAAAT‑TAAAAGG‑TTTTCCC‑TTAAAAA‑TTCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAggg                                                             <  2:147911/148‑1 (MQ=255)                                                                                  aCGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgggg                                                            >  1:26864/1‑151 (MQ=255)                                                                                    gggtttttttttGGCTATGCCCACACCCCCCCAGTAAAGGA‑AATATGG‑TTTCCCT‑TAAAAAT‑TT‑CCCGCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgggg                                                            <  2:1190706/142‑1 (MQ=255)                                                                                     ggttttttttgGGAGATGCCCGCATACCCCCGGAAAAATA‑AAAAATT‑TTTCCCC‑TAAAAAT‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcag                                                         <  1:51088/145‑1 (MQ=255)                                                                                        tgctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagca                                                       >  1:830744/1‑150 (MQ=255)                                                                                           tttttGGCGTTCTCCGCATACCCCCAGGTAAAGA‑AAAAGGT‑TTTACCC‑TAAAAAA‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAt                                                    <  2:371620/150‑1 (MQ=255)                                                                                           tttttGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtg                                                   >  2:1049704/1‑151 (MQ=255)                                                                                                    aTGCCCGCATCGCCCCAGGAAAAGA‑AAAAAGT‑TTTACCT‑TAAAAAA‑TA‑CCCGCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                            <  1:645061/149‑1 (MQ=255)                                                                                                       ccccaaTACCCCCAGAAAAAGA‑AATAAGT‑TTTACCT‑TAAAAAA‑TC‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAaga                                         <  2:120747/144‑1 (MQ=255)                                                                                                             tAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCggga                                  >  1:999865/1‑148 (MQ=255)                                                                                                             tAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCgg                                    >  2:931739/1‑148 (MQ=255)                                                                                                             accccccAAGAAAAAA‑AAAAAGT‑TTTCCCC‑TAAAAAATTT‑CCCGCACCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                   <  2:588732/147‑1 (MQ=255)                                                                                                              taccccAAGAAAAAG‑AAAAAGGTTTTCCCC‑TAAAAAA‑TC‑CCCCCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                   <  1:602354/147‑1 (MQ=255)                                                                                                              aGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaa                                 >  2:718744/1‑150 (MQ=255)                                                                                                              aGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCACCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaat                                >  1:216880/1‑151 (MQ=255)                                                                                                                ccccAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatt                               >  2:1147872/1‑150 (MQ=255)                                                                                                                ccccAAGTAAAAT‑AAAAAGT‑TTTTCCT‑GTAAAAA‑TTCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatt                               <  2:1166778/151‑1 (MQ=255)                                                                                                                 cccAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatt                               >  1:1146826/1‑149 (MQ=255)                                                                                                                 cccAAGTAAAATAAATATGT‑TTTCCCC‑TAAAAAT‑TT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatta                              <  2:760690/151‑1 (MQ=255)                                                                                                                  ccAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaa                             >  2:339081/1‑150 (MQ=255)                                                                                                                   aaGGAAAAGA‑AATAAGG‑TTTCCCC‑TAAAAAT‑TC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaata                           <  1:454300/150‑1 (MQ=255)                                                                                                                    aGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaata                           >  1:49372/1‑150 (MQ=255)                                                                                                                      gTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                        >  1:94716/1‑151 (MQ=255)                                                                                                                          aaagaaaaaatttttCCCC‑TAATAAA‑TC‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaata                           <  2:471794/134‑1 (MQ=255)                                                                                                                            t‑aataaGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                     >  2:191045/1‑148 (MQ=255)                                                                                                                                     gtttctccaggATAAA‑TG‑CCCCCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgt              <  2:1198031/139‑1 (MQ=255)                                                                                                                                      tttACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtt          >  2:489719/1‑151 (MQ=255)                                                                                                                                          ccc‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtttt        >  2:52032/1‑149 (MQ=255)                                                                                                                                           cc‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGCAATTAATAAGGCGCAGAGCGGCGTTTTGa      >  2:217906/1‑150 (MQ=255)                                                                                                                                            c‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGaa     >  2:1190624/1‑150 (MQ=255)                                                                                                                                               gATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACa   >  2:560854/1‑150 (MQ=255)                                                                                                                                                  aaa‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt  <  1:317362/148‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                               CCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGG‑TTTACCC‑TGATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAT  >  NC_000913/4378655‑4378944

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