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

CCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑CGTTTGCTGTGTGACGGGCTGCTTTTTGGCG‑ATGTCCGC‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTC  >  NC_000913/4378654‑4378950                                                                                                                                                         |                                                                                                                                                       cccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑CGTTTGCTGTGTGACGGGCTGCTTTTTGGCG‑ATGTCCGC‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTa‑cc                                                                                                                                                     <  2:139220/151‑1 (MQ=255)  ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑CGTTTGCTGTGTGACGGGCTGCTTTTTGGCG‑ATGTCCGC‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTa‑ccg                                                                                                                                                    <  2:1081601/151‑1 (MQ=255)  ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑CGTTTGCTGTGTGACGGGCTGCTTTTTGGCG‑ATGTCCGC‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTa‑ccg                                                                                                                                                    <  2:308830/151‑1 (MQ=255)    tCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑CGTTTGCTGTGTGACGGGCTGCTTTTTGGCG‑ATGTCCGC‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTa‑ccgcc                                                                                                                                                  >  1:768432/1‑151 (MQ=255)      gAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑CGTTTGCTGTGTGACGGGCTGCTTTTTGGCG‑ATGTCCGC‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTa‑ccgcca                                                                                                                                                 <  2:627337/150‑1 (MQ=255)                gTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑CGTTTGCTGTGTGACGGGCTGCTTTTTGGCG‑ATGTCCGC‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATAccc                                                                                                                                     >  1:184679/1‑151 (MQ=255)                gTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑CGTTTGCTGTGTGACGGGCTGCTTTTTGGCG‑ATGTCCGC‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATAccc                                                                                                                                     >  1:786395/1‑151 (MQ=255)                  cGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑CGTTTGCTGTGTGACGGGCTGCTTTTTGGCG‑ATGTCCGC‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCgg                                                                                                                                   >  1:525043/1‑151 (MQ=255)                            cccgcaattAAATGTTTTGCTGACCAAACCTAACAAAAACAC‑CTTTGGCTGTGAGACGGGTTTTTTTTTGGCT‑ATGCCCGC‑ATCCCCCAAGAAAAAAAAAAATGT‑TTTACCCTAAAAATTTA‑CCGCCG‑CCGCCATACCCGGGATCtttt                                                                                                                           <  1:996647/141‑1 (MQ=255)                            cACGCCAGAAATTTTTTGGCTGACCGACCCAAACTAAAACAG‑CTTTTTTTGTGTGAGGGGTTTTTTTTTGGCG‑ATCCCCGC‑ATACCCCCAGGAAAAGAAATAAGT‑TTTCCCCTAAAAATTTC‑CCCCCA‑CCGCCATACCCGGGATc                                                                                                                               <  1:602468/145‑1 (MQ=255)                             cccccATTAATTTTTTTGCTGCCCAACCCTAACACAAACAC‑CGTTTGCTGTGTGACGGCCTGTTTTTTGGCG‑ATGTCCGC‑ATCCCCCCAGGTAAAAAAATAATT‑TTTCCCCTGAAAAATCC‑CCCCCA‑CCGCCATACCCGGGATCTTttgtt                                                                                                                        <  1:313697/148‑1 (MQ=255)                                 cAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑CGTTTGCTGTGTGACGGGCTGCTTTTTGGCG‑ATGTCCGC‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGctc                                                                                                                    >  1:186507/1‑151 (MQ=255)                                 aattAAATTTTTTGCGGCCCGCACCTAAC‑CAAACAATCTTTTCTTGTGTGTGGGGCTG‑TTTTTGGTG‑GTTATCGCAATAGCCCCGGGAAAAATAATAAGG‑TTTGCCTTGATAAATTC‑CCGCCA‑CCGCCATACCCGGGATCTTttgttg                                                                                                                       <  2:689365/145‑1 (MQ=255)                                     aattttttgtcttcccGACCCAAACTCAAACAA‑CTTTGGTTGGGGGA‑GGGTCGTTTTTTGGGGAATGCCCCC‑TTACCCCAAGGAAAAGAAATAAGT‑TTTCCCCGAAAAATTTC‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCtata                                                                                                                <  2:716005/137‑1 (MQ=255)                                     aaaTGTTTTGCTGACCGAACCTAACTCAAACAA‑CGTTTGCTGTGTGACGGGCTGCTTTTTGGCG‑ATGTCCGC‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                                 >  1:276042/1‑150 (MQ=255)                                     aaaTGTTTTCCTGACGAAACTTAACTCAAACAC‑CTTTTGTTGTGTGAGGGGTTGTTTTTGGGCG‑ATCCCCGC‑ATCCCCCCAGAAAAAAAAAAAGGT‑TTTCCCCTGATAAATAC‑CCCCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                                 <  2:174766/150‑1 (MQ=255)                                        ttttttGCTGCCCAACCCAACCAAAAACAA‑CGTTTGTTGTGTGGCGGGCTTTTTTTGGGCG‑ATGCCCGC‑ATCCCCCCAGATAAAGAAAAAAGT‑TTTCCCCTAAAAAATTC‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAg                                                                                                              <  1:580803/148‑1 (MQ=255)                                        tGTTTTGCTGACCGAACCTAACTAAAACAC‑CGTTTGCTGTGAGACGGGTTGTTTTTTGGCG‑ATGTCCGC‑ATACCCCCAGTAAAAGAAATAAGT‑TTTCCCCTAATAAATTC‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAg                                                                                                              <  2:795848/150‑1 (MQ=255)                                          ttttGCGGCCCGAACCAAACTAAAACAA‑CTTTTGCTGTGGGACGGGCTGTTTTTTGGGG‑ATGTCCGC‑ATAGCCCCAGGAAAAATAATAAGG‑TTTCCCCTGAAAAATTC‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCg                                                                                                            <  2:1030464/150‑1 (MQ=255)                                            ttgcgaacctacccTAACAAAAACAA‑CTTTCGCTGTGGGACGGGCTTCTTTTTGGCA‑AGTTCCAC‑TAACCCCCAGGTAAA‑TAAAAAGGTTTTCCCCTAAAAAATTC‑CCGCCC‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCggg                                                                                                          <  2:639927/138‑1 (MQ=255)                                            ttGCTGACCGAACCTAACTCAAACAA‑CGTTTGCTGTGTGACGGGCTGCTTTTTGGCG‑ATGTCCGC‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCgggg                                                                                                         >  2:554279/1‑151 (MQ=255)                                                tGACCGAACCTAACTCAAACAA‑CGTTTGCTGTGTGACGGGCTGCTTTTTGGCG‑ATGTCCGC‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAt                                                                                                     >  1:891477/1‑151 (MQ=255)                                                 cacaaaaCCTAACTAAAACAC‑GTTTTGTTGTGGGAGGGGTTGTTTTTGGGCG‑TTGCCCGA‑ATCCCCCCAGGTAAAGTAATAAGT‑TTTCCCTTGA‑AAATTACCCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTaa                                                                                                      <  2:123532/144‑1 (MQ=255)                                                   ccaacccaaaCTAAAACAC‑CGTTTGGTGTGGGAGGGGCTGTTTTTGGGGG‑ATGCCCGC‑ATAGCCCAAGTTAAATTAATAGGT‑TTTCCCCTGATAAATCC‑CCGCCC‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTa                                                                                                       <  2:689197/138‑1 (MQ=255)                                                   ccGACCCAAACAAAAACAC‑CTTTTGTTGTGTGACGGGCTTTTTTTTGGCG‑ATGCCCGC‑ATACCCCCAGGAAAAGAAAAAAGT‑TTTCCCCTGAAAAATTC‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAt                                                                                                     <  1:380052/148‑1 (MQ=255)                                                      acccTACCTCAAACAC‑CTTTTCCTGTGTGACGGTCTTCTTTTTGGCG‑TTGCCCCC‑ATCCCCCCTGATAAAGTAATAAGT‑TTTCCCCTAATAATTTC‑CCCCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGc                                                                                                <  2:520315/148‑1 (MQ=255)                                                      aaCCTAACTCAAACAA‑CGTTTGCTGTGTGACGGGCTGCTTTTTGGCG‑ATGTCCGC‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCg                                                                                               >  2:747334/1‑151 (MQ=255)                                                           aaCTCAAACAA‑CGTTTGCTGTGTGACGGGCTGCTTTTTGGCG‑ATGTCCGC‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                           >  1:81929/1‑150 (MQ=255)                                                                aaaCAA‑CGTTTGTTGTGTGACGGGCTGTTTTTTGGCG‑ATGTCCGC‑ATACCCCAAGGTAAAGTAATAAGT‑TTTCCCTTAATAATTTC‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAata                                                                                     <  2:291157/151‑1 (MQ=255)                                                                  acaa‑cTTTGGCTGTGTGAGGGGTTGTTTTTGGGCG‑ATGCCCGC‑ATGCCCCAAGGAAAGGAAATAAGT‑TTTCCCCAAATAATTCC‑CCCCCC‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAatat                                                                                    <  2:787002/150‑1 (MQ=255)                                                                   caa‑cGTTTGCTGTGTGACGGGCTGCTTTTTGGCG‑ATGTCCGC‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCg                                                                                  >  2:688030/1‑151 (MQ=255)                                                                      acttttGTTGTGTGCGGGGTTTTTTTTGGGCG‑TTGCCCGC‑ATACCCCCAGTTAAATTATTAAGT‑TTTCCCCTAATAATTTC‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCg                                                                                  <  2:840891/146‑1 (MQ=255)                                                                              tgtgtgACGGGCTGCTTTTTGGCG‑ATGTCCGC‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgtt                                                                        >  1:1066745/1‑151 (MQ=255)                                                                              tgtgtgACGGGCTGCTTTTTGGCG‑ATGTCCGC‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgtt                                                                        >  1:163394/1‑151 (MQ=255)                                                                                        gctgctTTTTGGCG‑ATGTCCGC‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagc                                                              >  1:349141/1‑151 (MQ=255)                                                                                         ctgctTTTTGGCG‑ATGTCCGC‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagc                                                              >  1:611293/1‑150 (MQ=255)                                                                                              ttttGGCG‑ATGTCCGC‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgt                                                        >  2:777869/1‑151 (MQ=255)                                                                                                   gggtttgcccGC‑ATCCCCCGAGAAAAAGTAATAAGT‑TTTACCCTAAAAAATTC‑CCCCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGa                                                     <  1:1118268/142‑1 (MQ=255)                                                                                                          cccGC‑ATAGCCCCAGTAAAAGTAATAAGT‑TTTCCCCTAATAAATTC‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAgg                                             <  1:882389/150‑1 (MQ=255)                                                                                                            cGC‑ATAGCCCCAGGAAAAGTAAAAAGT‑TTCACCTTGAAAAATTC‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAagag                                              <  2:851035/148‑1 (MQ=255)                                                                                                              c‑ATCGCCCAAGAAAAAGAAATAAGT‑TTTCCCTTAAAAAATCC‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCgg                                          <  1:1099934/150‑1 (MQ=255)                                                                                                              c‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                         >  2:595023/1‑151 (MQ=255)                                                                                                                aTAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                         >  1:363404/1‑150 (MQ=255)                                                                                                                    cccaagaaaaaGAAAAAGGG‑TTTCCCTTAAAAATTTC‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCg                                           <  1:712822/136‑1 (MQ=255)                                                                                                                     cccAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatta                                    >  2:1085537/1‑150 (MQ=255)                                                                                                                     cccAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatta                                    >  1:557631/1‑150 (MQ=255)                                                                                                                     ccaagataaataaataaGT‑TTTCCCTTAATAAATCC‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaa                                   <  1:415030/139‑1 (MQ=255)                                                                                                                      ccAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaa                                   >  1:78254/1‑150 (MQ=255)                                                                                                                       cAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaat                                  >  2:1064702/1‑150 (MQ=255)                                                                                                                          gTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                              >  1:247859/1‑151 (MQ=255)                                                                                                                           tAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgc                             >  1:749203/1‑151 (MQ=255)                                                                                                                            aaaGAAATAGGT‑TTTCCCTTAATAATTTC‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg                            <  2:95459/151‑1 (MQ=255)                                                                                                                             aaGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg                            >  2:134463/1‑150 (MQ=255)                                                                                                                             aaGAAATAGGT‑TTTCCCATGATAAATTC‑CCCCCG‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                           <  2:104835/151‑1 (MQ=255)                                                                                                                              aGTAATAAGT‑TTTCCCTTGATAATTTC‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCa                          <  2:280966/151‑1 (MQ=255)                                                                                                                                 aataa‑G‑TTTCCCCCGAAAAAATACCCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCaga                        <  2:961796/150‑1 (MQ=255)                                                                                                                                    aaGG‑TTTCCCCTGAAAAATTC‑CCCCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgt                    <  2:313437/151‑1 (MQ=255)                                                                                                                                    aaGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGAAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGGGGAGCAAGAGGCGGAAATTAATAAGGGGCagag                       >  1:95941/1‑148 (MQ=255)                                                                                                                                      gg‑TTTGCCCCGAAAAATTG‑CCCCCAGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCagag                       <  2:661183/147‑1 (MQ=255)                                                                                                                                      gg‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtcg                  >  1:952387/1‑151 (MQ=255)                                                                                                                                      aagttttccacaaaaaaTTA‑CCGCCAGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgt                    <  1:56359/137‑1 (MQ=255)                                                                                                                                       gttttCCCCTAATAAATTC‑CCGCCT‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGc                      <  1:815434/145‑1 (MQ=255)                                                                                                                                         tttACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtt                >  2:74823/1‑151 (MQ=255)                                                                                                                                          tttcccTGAAAAATTC‑CCCCCACCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtc                   <  1:630888/145‑1 (MQ=255)                                                                                                                                               cTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAAc          >  2:121429/1‑151 (MQ=255)                                                                                                                                                tGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAAc          >  2:926615/1‑150 (MQ=255)                                                                                                                                                    aaaTTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAgg     >  2:500745/1‑151 (MQ=255)                                                                                                                                                      aTTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg    >  1:304867/1‑150 (MQ=255)                                                                                                                                                       ttA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTc  >  2:950527/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                       CCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAA‑CGTTTGCTGTGTGACGGGCTGCTTTTTGGCG‑ATGTCCGC‑ATAGCCCCAGGTAAAGTAATAAGG‑TTTACCCTGATAAATTA‑CCGCCA‑CCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTC  >  NC_000913/4378654‑4378950

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