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

CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTG‑CTGACCGAACCT‑AACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTC  >  NC_000913/4378653‑4378950                                                                                                                                                              |                                                                                                                                                     ccccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTG‑CTGACCGAACCT‑AACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ac                                                                                                                                                     <  1:697692/151‑1 (MQ=255)   ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTG‑CTGACCGAACCT‑AACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑accg                                                                                                                                                   >  2:1038883/1‑151 (MQ=255)      cGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTG‑CTGACCGAACCT‑AACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑accgcc                                                                                                                                                 >  2:849102/1‑150 (MQ=255)       gAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTG‑CTGACCGAACCT‑AACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑accgccac                                                                                                                                               >  1:225673/1‑151 (MQ=255)           aaTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTG‑CTGACCGAACCT‑AACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑accgccaccgcc                                                                                                                                           >  1:1021311/1‑151 (MQ=255)           aaTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTG‑CTGACCGAACCT‑AACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑accgccaccgcc                                                                                                                                           >  1:810019/1‑151 (MQ=255)           aaTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTG‑CTGACCGAACCT‑AACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑accgccaccgcc                                                                                                                                           >  1:1161116/1‑151 (MQ=255)              aGCGTCGCGACCAAGCACGCCAGTAAATGTTTTG‑CTGACCGAACCT‑AACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCAt                                                                                                                                         >  1:911024/1‑150 (MQ=255)                   cGCCACCAAGCACGCCAGTAAATGTTTTG‑CTGACCGAACCT‑AACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCgg                                                                                                                                   >  2:681796/1‑151 (MQ=255)                       accaAGCACGCCAGTAAATGTTTTG‑CTGACCGTACCT‑AACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGAt                                                                                                                                >  1:1108254/1‑150 (MQ=255)                             cACGCCGTTAATTTTTTTG‑CTGCCCAACCAT‑AACACAAACAGCTTTTGCT‑GTGTGACGGGCTGTTTTTG‑GGTGATGCCCGAATACCCCGAGAT‑AATAAAAAAAGTTTTCCCTA‑T‑AAAAATA‑CCCGCCACCGCCATACCCGGGATCtt                                                                                                                             <  1:269462/147‑1 (MQ=255)                              aCGCCAGTAAATGTTTTG‑CTGACCGAACCT‑AACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTttgtt                                                                                                                        >  1:120671/1‑151 (MQ=255)                               cGCAAGTAATTGTTTGG‑CTGCCCGACCCA‑ACCCAAAACATCTTTTGCT‑GTGGGAGGGGTTTTTTTTT‑GGCGTTTCCCGCATACCCCCAGGA‑AAATTAATAAGTTTTCCCCT‑A‑AAAATTT‑CCCCCCACCGCCATACCCGGGATCTTttgtt                                                                                                                        <  1:501911/150‑1 (MQ=255)                                      aaaTTTTTTG‑CTGCCCAACCCA‑AACAAAAACAAGTTTTGTT‑GGGGGCGGGCCTTTTTTTT‑GGCGATGTCCGCATCCCCCCAGG‑‑AAAAAAAAAAGGTTTTCCCC‑GAATAAATT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCtata                                                                                                                <  2:868568/151‑1 (MQ=255)                                       aaTTTTTTT‑CTGACCCGACCTAAACCAAAACATCTTTGGGT‑GTGGGGCGGGTTTTTTTTG‑GGCGTTGCCCGAATACCCCCAGGA‑AAATTAAAAAGTTTTCCCTT‑A‑AAAATTT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCta                                                                                                                  <  2:48815/149‑1 (MQ=255)                                       aaTTTTTTGTCTGACCAAACTA‑ACCTAAAACAATTTTTGTT‑GTGTG‑CGGGCTTTTTTTTGGGGGTTCCCCGAATGCCCCAAGGA‑AAAATAATAGTTTTTCCCTT‑A‑AAAAATC‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGctct                                                                                                                   <  2:1204621/148‑1 (MQ=255)                                                     ccGAACCT‑AACAAAAACAACGTTGGCG‑GTGTGCGGGGCTGTTTTTT‑GGCGATGCCCGCATACCCCGAGTA‑AAAGAAATAAGTTTTCCCCT‑A‑AAAAATC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtg                                                                                                   <  2:114691/150‑1 (MQ=255)                                                      ccaaccttaacaaaaaCATCTTTTGTT‑GTGGGAGGGGTTGTTTTTT‑GGCGATGCCCGCATACCCCGAGGA‑AAAGAAATAGTTTTTCCCCT‑A‑AAAATTT‑CCCGCCACCGCCATACCCGAGATCTTTTGTTGCTCTATAAGCGGGGTAAtt                                                                                                    <  2:325055/136‑1 (MQ=255)                                                       gAACCT‑AACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtggg                                                                                                 >  2:570045/1‑148 (MQ=255)                                                       gAACCT‑AACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGc                                                                                                >  1:269171/1‑151 (MQ=255)                                                          ccT‑AACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGAt                                                                                             >  2:226251/1‑151 (MQ=255)                                                           cT‑AACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATg                                                                                            >  2:152690/1‑151 (MQ=255)                                                            t‑AACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                           >  2:1116097/1‑151 (MQ=255)                                                                   aaaCAACGTTTGCT‑GTGTGACGGGCTGCTTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                      >  2:1252089/1‑150 (MQ=255)                                                                   aaaCAA‑GTTTGGTGGTGTGAGGGGCTGTTTTTG‑GGCGATGTCCGCATACCCCCAGGA‑AAAGTAAAAGGGTTTCCCTT‑A‑ATAAATT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                      <  2:815108/150‑1 (MQ=255)                                                                         cGTTTGCT‑GTGTGACGGGCTGCTTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                               >  1:919219/1‑151 (MQ=255)                                                                           ttggtt‑gtgggAGGGGTTGTTTTTG‑GGCGATGCCCGCATAGCCCAAGGA‑AAAGTAATAAGGTTTCCCCT‑A‑AAAATTC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                                <  2:913562/143‑1 (MQ=255)                                                                               gctggatgggggtctgttTTTG‑GTCGATGTCCACAGCCCCCAAGGA‑AAGATAATATTGTTCCCCCT‑A‑ATAAATT‑TCCCTCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttg                                                                          <  2:970540/138‑1 (MQ=255)                                                                               cT‑GTGTGACGGGCTGCTTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgt                                                                         >  1:194844/1‑151 (MQ=255)                                                                                     tgACGGGCTGCTTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgg                                                                    >  2:1029397/1‑151 (MQ=255)                                                                                      gaggggTTGTTTTTT‑GGCGATGCCCGCATCCCCCCAGAA‑AAAG‑AAAAAAGTTTTCCCT‑GAAAAAATC‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgg                                                                    <  2:433771/147‑1 (MQ=255)                                                                                         ggtttttttttg‑GGCGTTGCCGGCATAGCCCCAGGA‑AAAAAAATAAGTTTTCCCTT‑A‑AAAAATT‑CCCGCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAttt                                                                           <  2:892048/133‑1 (MQ=255)                                                                                         ggtctttttttt‑GGCGATGCCCGCATCCCCCAAGAT‑AAATCAATAAGTTTTCCCCT‑A‑ATAAATT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAggg                                                                   <  2:44581/141‑1 (MQ=255)                                                                                         gggCTGTTTTTG‑GGAGATGTCCGCATAGCCCCAGAA‑AAGGGAAAAAGGTTTCCTTT‑A‑AAAATTT‑ACCGCAACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgca                                                                <  1:1081145/151‑1 (MQ=255)                                                                                           gctgctTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagc                                                              >  1:199944/1‑151 (MQ=255)                                                                                                  tttcggtgcggccgccATAGCCCGAGGT‑AAGTTAAAAAGTTTTCCCTT‑G‑ATAAATC‑ACCCCTTCCGCCATTCCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGTCAGCAAAtg                                                         <  1:589379/135‑1 (MQ=255)                                                                                                   tt‑GGCGGTGCCCGCAGCCCCCAAGGA‑AAAGAAATAGGTTTTCCCTT‑A‑AAAAATC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgtg                                                       <  1:961552/150‑1 (MQ=255)                                                                                                   tt‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                      >  2:86372/1‑151 (MQ=255)                                                                                                    t‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                      >  2:1095382/1‑150 (MQ=255)                                                                                                          aTGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                  >  1:379318/1‑149 (MQ=255)                                                                                                           tGTCCGCATAGCCCAAGAA‑AAAGAAATAAGTTTTCCCCT‑A‑ATAAATC‑CCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCaa                                                 <  2:403311/149‑1 (MQ=255)                                                                                                           atccccgcgcaccccGGGA‑AAAATAAGGAGGTTTTCCCTAA‑AAAAATT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCaa                                                 <  2:259073/139‑1 (MQ=255)                                                                                                              cggaacaccccaAGGA‑AAAATAAGAGGTTTTCCCAT‑G‑AAAATTT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAagag                                              <  1:979801/141‑1 (MQ=255)                                                                                                               cGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCg                                           >  2:943872/1‑151 (MQ=255)                                                                                                               cGCATACCCCCAGAA‑AAATAAATAGGTTTTCCCTT‑A‑AAAAATT‑CCCCCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGc                                            <  2:267344/150‑1 (MQ=255)                                                                                                                 cATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCgg                                          >  1:104494/1‑150 (MQ=255)                                                                                                                  aTACCCCCAGGT‑AAATTAATAAGTTTTCCCCT‑A‑ATAATTT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                         <  2:278907/150‑1 (MQ=255)                                                                                                                  aTACCCCCAGGA‑AAAGAAATAGTTTTTCCCCT‑A‑AAAAATT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                         <  2:621781/150‑1 (MQ=255)                                                                                                                    acccccAGGAAAAAGTAATAAGTTTCCCCCA‑A‑AAAAATT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                         <  1:217999/147‑1 (MQ=255)                                                                                                                       ccgagat‑AAAGTAATATGTTTTCCCTT‑G‑ATAATTT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatt                                     <  1:973298/143‑1 (MQ=255)                                                                                                                            gttaaataaatatGTTTTCCCCT‑G‑AATAATT‑CCCACCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaata                                 <  2:1021349/141‑1 (MQ=255)                                                                                                                             t‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                              >  1:711405/1‑150 (MQ=255)                                                                                                                               aaaTTAATAGGTTTTCCCCT‑A‑ATAAATC‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg                            <  1:495504/151‑1 (MQ=255)                                                                                                                               aaaGTAATAAGTTTTCCCTT‑A‑AAAATTT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg                            <  2:101919/151‑1 (MQ=255)                                                                                                                               aaaGTAATAA‑GTTTTCCCT‑T‑AAAAATTCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaataa                                <  1:849064/147‑1 (MQ=255)                                                                                                                                aaGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                           >  1:627185/1‑151 (MQ=255)                                                                                                                                 aGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgc                             >  1:319413/1‑148 (MQ=255)                                                                                                                                    aatagGTTTTCCCCT‑A‑AAAATTT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCaga                        <  2:1161150/150‑1 (MQ=255)                                                                                                                                    aataaGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCag                         >  1:494500/1‑149 (MQ=255)                                                                                                                                         ggTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtcg                  >  2:774060/1‑151 (MQ=255)                                                                                                                                              aCCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtttt              >  2:1122287/1‑150 (MQ=255)                                                                                                                                                 ccttaaTAAATC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGttt               <  1:186824/143‑1 (MQ=255)                                                                                                                                                  t‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGGTTTGAAc          >  1:1216522/1‑150 (MQ=255)                                                                                                                                                      aTAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGGGCAGAGCGTCGTTTTGaa           >  1:428826/1‑147 (MQ=255)                                                                                                                                                       aaaaTT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt        <  2:979916/148‑1 (MQ=255)                                                                                                                                                         aaTT‑GCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg    <  2:67416/151‑1 (MQ=255)                                                                                                                                                         aaTT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg    >  2:795469/1‑151 (MQ=255)                                                                                                                                                            t‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGGGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCCGAGCGTCGTTTTGAACATAggg    >  2:877433/1‑148 (MQ=255)                                                                                                                                                            t‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAACGTCGTTTTGAACATAGGGTc  >  1:473302/1‑150 (MQ=255)                                                                                                                                                              |                                                                                                                                                     CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTG‑CTGACCGAACCT‑AACTCAAACAACGTTTGCT‑GTGTGACGGGCTGCTTTTT‑GGCGATGTCCGCATAGCCCCAGGT‑AAAGTAATAAGGTTTACCCT‑G‑ATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTC  >  NC_000913/4378653‑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: ‑