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.7%	102.9 / 17.2	60	V48G (GTA→GGA)	ampC	penicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (31/4);  new base C (0/25);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 4.57e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAA‑CCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGG  >  NC_000913/4378655‑4378948                                                                                                                                                          |                                                                                                                                                    ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAA‑CCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑a‑ccg                                                                                                                                                 >  1:1042557/1‑151 (MQ=255)  cTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAA‑CCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑a‑ccgc                                                                                                                                                >  2:1104453/1‑151 (MQ=255)    cGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAA‑CCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑a‑ccgcc                                                                                                                                               >  1:1031883/1‑150 (MQ=255)    cGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAA‑CCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑a‑ccgcc                                                                                                                                               <  2:137760/150‑1 (MQ=255)     gggaaaTCGCGTCCCCCCCACGCCCGCCAGTAATGTTTTTGCGGCCCGAC‑CATAACTCAAACACCTTTTGCTGTGGGACGGGTTGTTTTTGGGGGATGTCCGCATAGCCCCAGG‑TAAAGTAAAAAGGTTTA‑CCC‑TTATAAA‑TT‑g‑ccgccac                                                                                                                                             <  1:98225/147‑1 (MQ=255)         aaTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAA‑CCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑a‑ccgccaccgcc                                                                                                                                         >  1:416728/1‑151 (MQ=255)          aTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAA‑CCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑a‑ccgtcaccgcc                                                                                                                                         >  2:868502/1‑150 (MQ=255)          aTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAA‑CCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑ccgccaccgcca                                                                                                                                        <  2:430449/151‑1 (MQ=255)                    caccaAGCACGCCAGTAAATGTTTTGCTGACCGAA‑CCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGAt                                                                                                                              >  1:345341/1‑151 (MQ=255)                    caccaAGCACGCCAGTAAATGTTTTGCTGACCGAA‑CCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGAt                                                                                                                              >  1:287354/1‑151 (MQ=255)                       caAGCACGCCAGTAAATGTTTTGCTGACCGAA‑CCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCtt                                                                                                                           >  1:85780/1‑151 (MQ=255)                        aaGCACGCCAGTAAATGTTTTGCTGACCGAA‑CCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCttt                                                                                                                          >  1:16058/1‑151 (MQ=255)                           cgcGCAAGAAATTGTTTGGCGCACCGAC‑CCAAACACAAAAAAGTTTTGTTGTGTGAGGGGTTTTTTTTTGGCGATGCCCGCACACCCCCGGA‑AAAAGAAATAAGTTTTA‑CCC‑TGAAAAT‑TT‑C‑CCCCCACCGCCATACCCGGGATCtttt                                                                                                                         <  1:558708/147‑1 (MQ=255)                            aCGCCAGTAAATGTTTTGCTGACCGAA‑CCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTttgtt                                                                                                                      >  1:230864/1‑151 (MQ=255)                               ccAGTAAATGTTTTGCTGACCGAA‑CCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCTGGGATCTTTTGTTGct                                                                                                                   >  1:33170/1‑151 (MQ=255)                                cAGTAAATGTTTTGCTGACCGAA‑CCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGctc                                                                                                                  >  2:733046/1‑151 (MQ=255)                                cAGTAAATGGTTTGCTGACCGAA‑CCTAACTCAAACAACGTTTGCTGGGTGACGGGCTGCTTTTTGGCGATGACCGCATAGCCATAGG‑GAAAGGAATAAGGTTTA‑GCC‑TGATAAG‑TT‑A‑CCGCGACCGCCATACCGGGGATCTTTTGTTGc                                                                                                                    >  1:1011835/1‑149 (MQ=255)                                    aatttttttGCTAACCGAC‑CCTACCTAAAACACTGTTTGGGGTGTGAGGGGCTTTTTTTTGGCGATGCCCCCATA‑CCCCAAG‑AAAAGGAAATAAGTTTT‑CCC‑TTAAAAA‑TTCC‑CCCCCACCGCCATACCCGGGATCTTTTGTTGctc                                                                                                                  <  1:20760/142‑1 (MQ=255)                                        tttttGCTG‑CCGAACCCAAACTAAAACAAGGTTTGTTGTGGGGGGGGCTTTTTTTGGGCGACGCCCGAATACCCCCAGG‑AAAAGTAA‑AAGGTTTACCCC‑TAATAAA‑TC‑C‑CCGCCACCGCCATACCCGGGATCTTttgttg                                                                                                                     <  2:899965/139‑1 (MQ=255)                                          tttGCTGCCCGAC‑CCAAACAAAAAAAATTTTTGTTGTGTGGCGGGCTGTTTTTTGGCGATGCCCGCATAGCCCCAGG‑TAAAGTAATAGGTTTTC‑CCC‑TAATAAA‑TT‑C‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCgg                                                                                                         <  2:258531/150‑1 (MQ=255)                                                  ccaaa‑CCAAACTAAAACACCGTTTGTTGTGTGAGGGGTTTTTTTTTGGCGATGCCCCCATACCCCCAGT‑AAAAGTAATAAGTTTTC‑CCC‑TAAAAAT‑TG‑C‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAtt                                                                                                  <  1:232893/146‑1 (MQ=255)                                                   cGAA‑CCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgtg                                                                                               >  1:255752/1‑151 (MQ=255)                                                     aa‑CCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCg                                                                                             >  1:24159/1‑151 (MQ=255)                                                        ccaaaCCCAAACAACGTTTGTTGTGTGAGGGGTTGTTTTTTGGCGTTGCCCGCATACCCCCAGG‑AAAATAAATAAGTTTTC‑CCC‑TAATAAT‑TT‑C‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgtg                                                                                               <  1:253970/144‑1 (MQ=255)                                                         cTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATg                                                                                          >  2:43517/1‑151 (MQ=255)                                                           aaccaaaaCACCTTTTGTTGTGTGACGGGCTGTTTTTTGGGGATGCCCGCATACCCCCGGA‑TAAATAAATAAGTTTTC‑CCT‑GGATAAA‑TT‑C‑CCGCCACCGCCATACCTGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                         <  2:33170/145‑1 (MQ=255)                                                             cTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac                                                                                       >  1:1092297/1‑150 (MQ=255)                                                                aaaCAACGTTTGCTGTGTGAGGGGCTGTTTTTGGGCGATGCCCCAATACCCCCGGG‑AAAAGAAATAAGTTTCC‑CCC‑TAATAAA‑TT‑C‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                    <  2:617293/150‑1 (MQ=255)                                                                  acaacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATc                                                                                 >  2:287832/1‑151 (MQ=255)                                                                    atcttttGCTGTGGGAGGGGCTGTTTTTGGGCGATGTCCGCATACCCCCAGA‑AAAAGAAATAAGTTTTC‑CCT‑TAAAAAT‑TT‑C‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCg                                                                                <  1:509592/146‑1 (MQ=255)                                                                    accttttGCTGTGGGACGGTCTTTTTTTTGGCGTTGCCCCCATAGCCCAAGA‑AAAAGTAATAGGTTTTC‑CCT‑TAATAAA‑TC‑C‑CCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGt                                                                               <  1:1000073/147‑1 (MQ=255)                                                                    aacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGt                                                                               >  2:374613/1‑151 (MQ=255)                                                                     acGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                              >  1:335464/1‑151 (MQ=255)                                                                      cGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                             >  1:570587/1‑151 (MQ=255)                                                                         ttGGTGTGTGAGGGGCTGTTTTTTGGGGAGGCCCGCATGCCCCAAGG‑AAAGGAAAAAAGTTTTA‑CCT‑TGAAAAA‑TT‑C‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAt                                                                           <  2:889008/150‑1 (MQ=255)                                                                              gtggggcgggttttttttgGGCGATGCCCGCAAACCCCGAGA‑TAAAGAAATAGGTTTTA‑CCT‑TAAAAAA‑TT‑C‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgt                                                                       <  1:912878/135‑1 (MQ=255)                                                                                 tgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgg                                                                  >  1:740801/1‑151 (MQ=255)                                                                                     gggtttttttttGGCGATGCCCGCATGGCCCCAGG‑AAAATAAATATGTTTTA‑CCC‑TAAAAAT‑TT‑C‑CCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgca                                                              <  1:174949/144‑1 (MQ=255)                                                                                          gctTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCaa                                                          >  2:984176/1‑150 (MQ=255)                                                                                           ctTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAt                                                        >  1:57523/1‑151 (MQ=255)                                                                                           ctTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAt                                                        >  2:296734/1‑151 (MQ=255)                                                                                            tttttGTCGATGCCCGCATA‑CCCCAGGATAAAGAAATATGTTTTA‑CCT‑TTAAAAA‑TT‑C‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCaaa                                                         <  1:47294/149‑1 (MQ=255)                                                                                                    gATGCCCCCATACCCCAAGA‑AAAAGAAAAAAGTTTTA‑CCT‑TAATAAA‑TT‑C‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCaa                                               <  1:887234/151‑1 (MQ=255)                                                                                                     aTGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                >  2:176215/1‑149 (MQ=255)                                                                                                     aTGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                >  1:258281/1‑149 (MQ=255)                                                                                                       gcccGCATCGCCCGAGG‑TAAAG‑AATTAGGTTTT‑CCCTTAATAAT‑TT‑C‑CCGCCGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAaga                                             <  2:341030/148‑1 (MQ=255)                                                                                                             acacccccAGG‑‑AAAATAAATAGGTTTT‑CCC‑TTAAAAATTT‑C‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGta                                      <  1:479365/147‑1 (MQ=255)                                                                                                                 ccccAGG‑AAAAGTAATAAGTTTTC‑CCC‑TGATAAT‑TT‑C‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaat                                    <  2:198924/149‑1 (MQ=255)                                                                                                                      gg‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaata                               >  1:433625/1‑149 (MQ=255)                                                                                                                       ggaaaataaataaGGTTTC‑CCT‑TAATAAT‑TC‑C‑CCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaata                               <  2:345400/141‑1 (MQ=255)                                                                                                                       gaaaaggTAATAAGGTTTC‑CCC‑TGATAAT‑TT‑C‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaataa                              <  1:37897/144‑1 (MQ=255)                                                                                                                       g‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                            >  1:361689/1‑151 (MQ=255)                                                                                                                          aaagaaaaaggttttC‑CCC‑TAAAAAT‑TT‑C‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaat                                <  1:117515/133‑1 (MQ=255)                                                                                                                          aaaGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg                          >  1:101729/1‑151 (MQ=255)                                                                                                                          aaaGGAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATttataggtcg                          >  1:957449/1‑146 (MQ=255)                                                                                                                           aataaatagGTTTTC‑CCC‑TAAAAAA‑TC‑C‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                         <  1:714161/147‑1 (MQ=255)                                                                                                                               aataa‑GTTTT‑CCC‑TTAAAAA‑TT‑ACCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCaga                      <  1:541906/150‑1 (MQ=255)                                                                                                                                ataaGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGc                    >  2:93923/1‑151 (MQ=255)                                                                                                                                   aGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtc                 >  2:739802/1‑151 (MQ=255)                                                                                                                                     ttttA‑CCC‑TGAAAAA‑TT‑C‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtc                 <  2:342490/148‑1 (MQ=255)                                                                                                                                       ttA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGttt             >  1:492947/1‑151 (MQ=255)                                                                                                                                           ccc‑TAATAAA‑TT‑C‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGttt             <  2:1080976/148‑1 (MQ=255)                                                                                                                                             cctgtaaaa‑TC‑A‑CCTCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGaa         <  2:108314/145‑1 (MQ=255)                                                                                                                                               tGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGaa         >  1:953570/1‑149 (MQ=255)                                                                                                                                                 aTAAA‑TT‑C‑CCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt      <  2:882622/150‑1 (MQ=255)                                                                                                                                                   aaa‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGaa         <  2:457623/145‑1 (MQ=255)                                                                                                                                                       tt‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg  >  1:1005474/1‑149 (MQ=255)                                                                                                                                                          |                                                                                                                                                    CCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAA‑CCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGG‑TAAAGTAATAAGGTTTA‑CCC‑TGATAAA‑TT‑A‑CCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGG  >  NC_000913/4378655‑4378948

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