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

CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACG‑TTT‑GCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTC  >  NC_000913/4378653‑4378950                                                                                                                                                                   |                                                                                                                                                       ccccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACG‑TTT‑GCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑a                                                                                                                                                        <  1:390032/150‑1 (MQ=255)   ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACG‑TTT‑GCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑acc‑g                                                                                                                                                    <  2:628884/151‑1 (MQ=255)       gaacAATGGGGTCGCCACCAAGCAAGCCAGTAAATGGTTTGCTGACCGAACCTAACTCAAACACCG‑TTT‑GCTGTGTGACGGGCTACTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑acc‑gcc                                                                                                                                                  <  2:1098336/146‑1 (MQ=255)       gAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACG‑TTT‑GCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑acc‑gcca‑c                                                                                                                                               >  2:691459/1‑151 (MQ=255)         gCAATAGGGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAGCG‑TTT‑GCTGTGTGACGGGCTGTTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTCCCC‑T‑GATAAA‑TT‑gcc‑gcca‑cc                                                                                                                                              <  2:832488/150‑1 (MQ=255)         gCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACG‑TTT‑GCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑acc‑gcca‑ccg                                                                                                                                             >  1:1247681/1‑151 (MQ=255)           aaTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACG‑TTT‑GCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑acc‑gcca‑ccgcc                                                                                                                                           >  1:605693/1‑151 (MQ=255)              ggCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACG‑TTT‑GCTGTGTGACGGGATGCTTTTTGGCGATGTACGCATAGCCCC‑A‑TG‑T‑AAAG‑T‑AATAG‑GGTTTACCC‑T‑GATAAA‑TT‑ACA‑GCCA‑CCGCCAt                                                                                                                                         >  2:1146363/2‑150 (MQ=255)              aGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACG‑TTT‑GCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCAt                                                                                                                                         >  1:377487/1‑150 (MQ=255)               gCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACG‑TTT‑GCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATAc                                                                                                                                       <  2:855086/151‑1 (MQ=255)                       accaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACG‑TTT‑GCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATc                                                                                                                               >  1:866452/1‑151 (MQ=255)                          aaGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACG‑TTT‑GCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCttt                                                                                                                            >  1:279000/1‑151 (MQ=255)                              aCGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACG‑TTT‑GCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTttgtt                                                                                                                        >  1:195440/1‑151 (MQ=255)                               cGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACG‑TTT‑GCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTttgttg                                                                                                                       >  1:650854/1‑151 (MQ=255)                                      aatttttttGCTGCCCGACCCTACCAAAAACAACT‑TTT‑GTTGTGTGCGGGTTTTTTTTTTGGCGATCCCCGCATACCCCC‑A‑GG‑A‑AAAG‑A‑AAAAA‑GTTTTCCCT‑T‑AAAAAA‑TT‑CCC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCtata                                                                                                                <  1:829094/146‑1 (MQ=255)                                      aatttttttGCTGACCAACCCAAACTAAAACCCGT‑TTT‑GCTGTGTGAGGGTCTGTTTTTTGGCAATCCCCGCATCCCCCC‑A‑GA‑T‑AAA‑‑T‑AATAAGGGTTTCCCC‑T‑AATAAT‑TC‑CCC‑CCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCta                                                                                                                  <  1:1260072/144‑1 (MQ=255)                                        aTGTTTTGCTGACCGAACCTAACTCAAACAACG‑TTT‑GCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCtata                                                                                                                >  1:426963/1‑149 (MQ=255)                                         ttttttCCTGCCCAAACCAACCACAAACATCT‑TTT‑TTTGTGGGGGGGGTTTTTTTTGGGCGATCCCCCCATCCCCCA‑A‑GT‑A‑AAAG‑T‑AATAT‑GTTTTCCCT‑T‑AAAAAA‑TT‑CCC‑GCAA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAg                                                                                                              <  2:846090/148‑1 (MQ=255)                                         tGTTTTGCTGACCGAACCTAACTCAAACAACG‑TTT‑GCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGc                                                                                                             >  2:107298/1‑151 (MQ=255)                                           ttttGCTGACCGACCCAAACTAAAAAACCG‑TTT‑GTTGGGGGCGGGGTTGTTTTTGGGGGATGCCCGAATACCCCC‑A‑‑G‑A‑AAAA‑T‑AAATA‑GGTTTTCCC‑T‑TAAAAATTT‑CCC‑CCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAg                                                                                                              <  1:74494/148‑1 (MQ=255)                                           ttttGCTGACCGAACCTAACTCAAACAACG‑TTT‑GCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCgg                                                                                                           >  2:460426/1‑151 (MQ=255)                                           ttttGCGGCCCAACCCTACCTCAAAC‑ACG‑TTTGGTTGTGTGGGGGGTTTTTTTTGGGGGATGCCCGCATACCCCC‑A‑GA‑A‑AAAG‑A‑AAAAG‑GTTTTCCCC‑T‑AAAAAA‑TT‑CCC‑CCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCg                                                                                                            <  1:409439/150‑1 (MQ=255)                                             ttgcgacccGAACATACCTAAAACACCT‑TTT‑GTTGTGGGAGGGGTTTTTTTTGGGCGATGCCCGCAACCCCCC‑A‑GA‑A‑AAAG‑T‑AATAA‑GTTTTCCCT‑T‑AAAAAA‑TC‑CCC‑CCCG‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCgg                                                                                                           <  2:600324/142‑1 (MQ=255)                                                   aCCGACCCTACCACAAACATCT‑TTT‑GCGGGGTGCCGGGCTTTTTTTTGGCGTTGCCCGCTCCCCCCG‑A‑GG‑A‑AAAG‑T‑AATAG‑GTTTTCCCC‑T‑AATAAA‑TC‑CCC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGt                                                                                                        <  1:495902/146‑1 (MQ=255)                                                   aCCGAACCTAACTCAAACAACG‑TTT‑GCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtg                                                                                                   >  1:184687/1‑151 (MQ=255)                                                      gAACCTAACTCAAACAACG‑TTT‑GCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑A‑AAAG‑T‑AATAA‑GGATTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgtg                                                                                                 >  2:966428/1‑150 (MQ=255)                                                       aaCGTAACTAAAACAGCT‑TTT‑GTGGTGTGACGGGTTG‑TTTTTTGCGAATCGCCCAAACCCCCGA‑GT‑T‑AAAG‑A‑AAAAA‑GTTTCCCCT‑T‑GAAAAA‑TC‑GCC‑CCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgtg                                                                                                 <  1:20310/149‑1 (MQ=255)                                                        accctaacaaaaCAACGTTTT‑GTTGTGTGAGGGGCTTTTTTTTGGCGATGCCCAAATACCCCC‑A‑‑G‑A‑AAAG‑TAAATAG‑GGTTTCCCT‑A‑AAAAAA‑TT‑CCC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgtg                                                                                                 <  1:501271/140‑1 (MQ=255)                                                         ctaaatcaaaaCAGCG‑TTT‑GTTGTGTGAGGGGCTGTTTTTGGGCGATGCCCG‑ATTACCCC‑A‑AG‑A‑AAAGTT‑AATAA‑GTTTTCCCT‑T‑AAAAAA‑TT‑CCC‑CCCG‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGa                                                                                              <  2:1127044/142‑1 (MQ=255)                                                          caaccaaaaccaacT‑TTT‑GCTGTGGGACGGGCTTTTTTTTGGCGATGCCCCCATCCCCCA‑A‑GA‑A‑AAAT‑T‑AATAA‑TTTTTCCCT‑A‑AAAAAA‑TT‑CCC‑GCAA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgt                                                                                                  <  2:141494/135‑1 (MQ=255)                                                               ccAAACAACG‑TTT‑GCTGGGTGATGGGCTTTTTTTTGGCGATGCCCGCAT‑CCCCC‑A‑GGTT‑AAAG‑A‑AATAA‑GGTTTCCCT‑T‑AATAAA‑TT‑CCC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACaa                                                                                       <  1:95977/150‑1 (MQ=255)                                                                 aaaCACTG‑TTG‑GTTGTGTGGCGGGTTTTTTTTTGGCGATGCCCGCATACCCCC‑A‑GA‑A‑AAAG‑T‑AA‑AA‑GGTTTTCCCTT‑AATAAA‑TT‑CCC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAata                                                                                     <  1:838075/151‑1 (MQ=255)                                                                 aaaCACCG‑TTT‑GCTGTGTGACGGGCTGTTTTTTGGCGTTGCCCGCATACCCCA‑A‑GG‑A‑AAAG‑T‑AATAA‑GGTTTCCCT‑T‑GAAAAA‑TT‑CCC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACaa                                                                                       <  1:1040482/149‑1 (MQ=255)                                                                 aaaCAATG‑TTG‑GCTGTGGGAGGGGTTGTTTTTGGGCGTTGCCCGCATCCCCCC‑A‑GG‑A‑AAAT‑T‑AATAG‑GTTTTCCCC‑T‑AAAAAA‑TT‑CCC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcaca                                                                                        <  1:751977/148‑1 (MQ=255)                                                                 aaaCAACG‑TTT‑GCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                      >  2:321459/1‑150 (MQ=255)                                                                   acaagG‑TTT‑GTTGTGGGACGGGCTGTTTTTGGGCGATGCCCGCATACCCCC‑A‑GG‑A‑AAAT‑A‑AATAG‑GTTTTACCC‑T‑AAAAAA‑TT‑CCC‑GCAA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATc                                                                                   <  1:1357142/151‑1 (MQ=255)                                                                     aacG‑TTT‑GCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGt                                                                                 >  1:479335/1‑151 (MQ=255)                                                                           tttgggtgggggcGGGCTGTTTTTTGGCTATGCCCGCAT‑CCCCC‑A‑GG‑TAAAAT‑T‑AATAA‑GTTTTCCCT‑T‑AATAAT‑TT‑ACC‑GCCT‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                               <  2:928368/137‑1 (MQ=255)                                                                                  tgtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                       >  2:754833/1‑150 (MQ=255)                                                                                  ggtgAGGGGCTGTTTTATGGGGAAGCCCGC‑CCGCCCC‑C‑AG‑T‑AAAT‑G‑AATTA‑GGTTTACCC‑A‑GTAAAA‑AC‑ACCGGACA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgt                                                                         <  2:481942/147‑1 (MQ=255)                                                                                    tgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgg                                                                    >  2:27910/1‑151 (MQ=255)                                                                                       cGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                                 >  1:213360/1‑151 (MQ=255)                                                                                               tttttGCCGTTGTCCCCAGCCCCCA‑A‑GT‑A‑AAAG‑A‑AATAA‑GTTTTACCT‑T‑AAAAAT‑TC‑CCC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCaaa                                                           <  2:697608/149‑1 (MQ=255)                                                                                                tttgGGCGATGCCCGCATACCCCC‑G‑GG‑A‑AAAG‑T‑AAAAA‑GTTTTCCCC‑T‑GAAAAT‑TT‑CCC‑CCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtg                                                         <  1:736864/150‑1 (MQ=255)                                                                                                          gccccaaTAGCCCG‑A‑GG‑A‑AAAG‑T‑AA‑AA‑GGTTTTCCC‑TAAAAAAA‑TA‑ACC‑GCCG‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                  <  1:946890/141‑1 (MQ=255)                                                                                                          gTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAag                                                >  2:324483/1‑149 (MQ=255)                                                                                                           tCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGCAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAgg                                             >  2:794310/1‑151 (MQ=255)                                                                                                            ccGCATACCCCC‑A‑GG‑A‑AAAG‑T‑AAAAA‑GTTTTCCCC‑T‑AATAAA‑TT‑CCC‑CCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                  <  1:577885/145‑1 (MQ=255)                                                                                                             ccgcaaccccc‑AGGG‑T‑AAAG‑A‑AATAA‑GTTTTCCCT‑T‑GATAAA‑TT‑CCC‑CCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGc                                            <  2:284199/144‑1 (MQ=255)                                                                                                             cGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCg                                           >  1:610177/1‑151 (MQ=255)                                                                                                              cgcttacccc‑A‑GG‑GAAAAT‑T‑AATAA‑GTTTTCCCT‑T‑AAAAAT‑TT‑CCC‑CCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCg                                           <  1:1150597/145‑1 (MQ=255)                                                                                                                catccccc‑C‑GG‑A‑AAAA‑T‑TAAAA‑GGTTTTCCC‑T‑TAAAAA‑TTCCCC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCgg                                          <  2:1357981/146‑1 (MQ=255)                                                                                                                 tAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGta                                        >  1:854652/1‑150 (MQ=255)                                                                                                                     ccc‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatta                                    >  2:139526/1‑150 (MQ=255)                                                                                                                     ccc‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatta                                    >  1:754134/1‑150 (MQ=255)                                                                                                                     ccc‑A‑GG‑A‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatt                                     >  2:1011793/1‑149 (MQ=255)                                                                                                                        cccatgaaaaaG‑T‑AATAT‑TGTTCCCCC‑T‑GATAAA‑TT‑CCC‑CCCT‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatta                                    <  1:288999/143‑1 (MQ=255)                                                                                                                         a‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACACTATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatt                                     >  2:24251/1‑146 (MQ=255)                                                                                                                         a‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaat                                  >  2:71302/1‑149 (MQ=255)                                                                                                                                aaaG‑T‑AATAA‑GTTTTCCCC‑T‑AAAAAA‑TT‑CCC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg                            <  1:779741/151‑1 (MQ=255)                                                                                                                                  aG‑T‑AATAA‑GGTTTCCCT‑T‑AAAAAT‑TT‑GCC‑CCCC‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                           <  2:437999/150‑1 (MQ=255)                                                                                                                                   aggtaatgattttttCCCC‑T‑AATAAA‑TT‑‑TC‑GCCACCCGCCATACCCGGGATCTTTCGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg                            <  2:5068/139‑1 (MQ=255)                                                                                                                                        ataa‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGGGCag                         >  2:107608/1‑148 (MQ=255)                                                                                                                                        ataa‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCagag                       >  1:1061037/1‑150 (MQ=255)                                                                                                                                        aaaaggttttCCCC‑T‑AAAAAT‑TT‑CCC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCag                         <  1:346011/142‑1 (MQ=255)                                                                                                                                              gTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtcgt                 >  1:589312/1‑151 (MQ=255)                                                                                                                                               tttCCCC‑T‑GAAAAA‑TT‑CCC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgt                    <  1:1139349/147‑1 (MQ=255)                                                                                                                                               tttACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtt                >  2:1143374/1‑151 (MQ=255)                                                                                                                                                   ccc‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg             >  1:1008158/1‑150 (MQ=255)                                                                                                                                                     c‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGaa           >  2:882693/1‑150 (MQ=255)                                                                                                                                                      ctagaaaaa‑TT‑CCC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg             <  1:1059953/143‑1 (MQ=255)                                                                                                                                                       t‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGaa           >  2:1314406/1‑149 (MQ=255)                                                                                                                                                                 t‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGGGCGGCGTTTTGAACATAGGGTc  >  2:461195/1‑150 (MQ=255)                                                                                                                                                                   |                                                                                                                                                       CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACG‑TTT‑GCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCC‑A‑GG‑T‑AAAG‑T‑AATAA‑GGTTTACCC‑T‑GATAAA‑TT‑ACC‑GCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTC  >  NC_000913/4378653‑4378950

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