Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,248,424	0	C	A	42.4%	79.2 / 32.7	66	V37F (GTT→TTT)	lysP	lysine transporter
Reads supporting (aligned to +/- strand):  ref base C (25/13);  new base A (0/28);  total (25/41)
Fisher's exact test for biased strand distribution p-value = 8.26e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGT  >  NC_000913/2248277‑2248571                                                                                                                                                       |                                                                                                                                                    cGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGcaacaa                                                                                                                                                   >  2:118734/1‑150 (MQ=255)  gTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACaaa                                                                                                                                                  >  1:941259/1‑150 (MQ=255)   tAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACaaaa                                                                                                                                                 >  2:775394/1‑150 (MQ=255)        ggCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCt                                                                                                                                           <  2:191082/151‑1 (MQ=255)           aaaCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCGTCAGGAAGTAAACC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGt                                                                                                                                         >  1:154301/1‑150 (MQ=255)           aaaCGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGt                                                                                                                                         >  1:581307/1‑150 (MQ=255)              cGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCg                                                                                                                                     >  1:709563/1‑151 (MQ=255)                aaCCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGAt                                                                                                                                   >  1:1008487/1‑151 (MQ=255)                    ggAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGaa                                                                                                                               <  2:785783/151‑1 (MQ=255)                     gAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCAACAGGAAGTAAACC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAAc                                                                                                                              >  1:453018/1‑151 (MQ=255)                      aaaCCGGCATATAGGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAACCC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAAc                                                                                                                              <  2:1147993/150‑1 (MQ=255)                             cATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATg                                                                                                                      >  1:262371/1‑151 (MQ=255)                                     cAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAatcat                                                                                                              >  1:506084/1‑151 (MQ=255)                                           gtccACCGAGACTGGTCATAAGGAAGTAACCC‑ATCAGGCCAATCAGAATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCa                                                                                                         <  1:894310/147‑1 (MQ=255)                                           gtccACCGAGACTGGTCATAAGGAAGTAACCC‑ACCAGGCCAATCAGCAAAAACGAGAGAAACCCCCCGCCCGGCCCTGCC‑TGAGAAAC‑CGTTGCGCCAGAGGCAAAAAAAAGCCCTGTACGGAGGAAACCGCCAATGGCAATCATCGTCAg                                                                                                        <  2:281919/148‑1 (MQ=255)                                           gTTCACCGAGACTGGTCATCAGGAAGTAAACC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTc                                                                                                          >  2:378625/1‑149 (MQ=255)                                               aCCGAGACTGGTCATCAGGAAGTAAACC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAgg                                                                                                       >  1:956149/1‑148 (MQ=255)                                                ccGGGACTGGTCTTCAGAAAGAAAACC‑ACCAGGCAAATAAACAAAGAGGAGAGAAACGCCCCCCCCGGCCCGGCC‑TGAGAAAT‑CGTGGCGCCAGAGCAAAAAAAAAGCCCTTTCCGGAGGAAACCGCCAATGGCAATCATCGTCAgg                                                                                                       <  2:807195/147‑1 (MQ=255)                                                     aCTGGTCATCAGAAAGTAACCC‑ATCAGGCCAACCAGCATATACGAGAGCAATCCCCCGCCCGGCCCTGCC‑GGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGCCTTGTACCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGAcgc                                                                                                 <  1:225991/148‑1 (MQ=255)                                                        gttcaccaGAAAGTAAACC‑ACCAGCCCAATCAGAATATACGAGACCAATCCCCCCCCCGGACCTGCC‑AAAAAAAT‑CGTGCCCCCAGAGCCAAAAAAAAGCCCTGTCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGAcgc                                                                                                 <  1:421437/143‑1 (MQ=255)                                                        ggcaaacaGGAAGTAACCC‑ATCAGGCCAACCAGAAAATACGAGAGCAATGCCCCCCCCGGACCTGCC‑GAAAAAAT‑CGTTCCGCAAGAGGAAAAAAAAAGCCCTGTCCCGTTGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCttt                                                                                           <  2:540691/145‑1 (MQ=255)                                                         gccaaggaggaatAAACC‑ATCAGGCCAATCAGAATATACGAGACAAACGCCGCACCGCCCCCAGCC‑AGAGAAAT‑CCTGGCGGAAGATGAAACAAAAAGCCTGGTGCCTTTAGGACCGCCAATGGCAATCATCGTCAGGTGACGCGCCttt                                                                                           <  2:523241/138‑1 (MQ=255)                                                           catcaGGAAGTAAACC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAAt                                                                                        >  1:1030441/1‑151 (MQ=255)                                                            atcaGGAAGTAAACC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAAtt                                                                                       >  2:977829/1‑151 (MQ=255)                                                             acaGAAAGTAACCC‑ATCAGGCCAACCAACATATACGAGAGAAACCCCCCGCCCGGACCGGCC‑TGAGAAAT‑CTTCGCGCGAGAGCCAACAAAAAGACCTGTGCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAAtt                                                                                       <  2:733828/149‑1 (MQ=255)                                                                  aaGTAACCC‑ACCAGCCCAATCAGAATATACGAGAGCAATCCCCCGCCCGGCCCTGCC‑TGAAAAAT‑CGTTGCGCCAGAGGCAACAAAAAGCCCTGTCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACg                                                                                   <  1:305258/149‑1 (MQ=255)                                                                  aaGTAAACC‑ATCAGGCCAATCAGCATATACGAGAGCAATCCCCCCCCCGGCCCTGCC‑TGAAAAAT‑CGTTGCGCCAGAGGAAAAAAAAAGCCCTGTCCCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACgg                                                                                  <  1:735264/150‑1 (MQ=255)                                                                   agaaaccc‑accaGCCCAACCAGAATATCCGAGGGAAAACCCCCGCCCGGCCCGGCG‑GGAAAAAT‑CTTGGCCCCAGAGCCAAAAAAAGGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCg                                                                                <  2:844390/145‑1 (MQ=255)                                                                    agtaacccctggggCAAATCAGAATATAGGAGAGAATAGCCCCCCCCGG‑CCTGCCGTGAGAATT‑CGTTGCGCCGAAGGCAACAAAAAGCCCGGGCCAGATGAACCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACgg                                                                                  <  1:217017/137‑1 (MQ=255)                                                                     aaaccc‑accaGCCAAACCAGCATATAGGAGAGCAACGCCCCGCCCGGCCCTCCC‑TGAAAAAC‑CGTGGCGCCAGAGGCAAAAAAAACACCTGTCCCGTTGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTa                                                                              <  2:936194/147‑1 (MQ=255)                                                                      aaaCCAATCAGGCAAAGAAGCATATACGAGAGCAACCCCCCGCCGGGACCGGCC‑GGAGAAAT‑GTTGGCGCGGGAGGCAAAAAAAAGCCGTGTACCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGt                                                                               <  2:275588/150‑1 (MQ=255)                                                                       aacccaccggcccAATCAGCATACAGGAGAGCAATCCCCCCCCCGGACCTGCC‑GGAGAAAT‑CGTGGCCCAGGAGGAAAAAAAAGGACCTGTACCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTaa                                                                             <  2:532087/140‑1 (MQ=255)                                                                        aCC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGc                                                                           >  2:656948/1‑151 (MQ=255)                                                                        aCC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGc                                                                           >  1:220608/1‑151 (MQ=255)                                                                                ggCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGgcgc                                                                    >  2:51318/1‑151 (MQ=255)                                                                                    aaacAGTATACAGGAGAACAATCCCCCGCCGGGACC‑GCC‑CGAAAAATCCGTGGCACCAGAGGCAACAAAAAGCCCTGACCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCt                                                                <  2:483316/148‑1 (MQ=255)                                                                                    aaTCGGCATATAGGAGAGAAATCCCCCCCCGGGCCCTGCC‑GGAAAATT‑CGTTGCGCCAGAGGCAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGgc                                                                      <  1:308122/145‑1 (MQ=255)                                                                                    aaTCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCt                                                                >  2:411831/1‑151 (MQ=255)                                                                                    aaTCAGCATATACGAGAGCAATCCCCCGCCCGGCCCTGCC‑TGAAAAAT‑CGTTGCGCCAGAGGAAAAAAAAAGCCCTGTCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCgg                                                                       <  2:994583/144‑1 (MQ=255)                                                                                           atatACGAGAGCACCGCCCCCCCCGACCCGGCC‑TGAAAATT‑CGTGGCGCCAGAGGAAACAAAAAGACTTGTACCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtt                                                          <  2:578381/150‑1 (MQ=255)                                                                                           atatACGAGACAATTCCCCCGCCGGGCCCTGCC‑GGAAAATT‑GGTTGCGCCAGAGCCAAAAAAAAGCCCTTCACCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGttt                                                         <  2:644657/151‑1 (MQ=255)                                                                                           aaataggaaaccAATCCCCCGCCCGACCCGGCC‑TGAAAATT‑CTTGGCGCAAGAGGAAAAAAAAAGCCCTGTACCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtgtg                                                             <  2:188871/136‑1 (MQ=255)                                                                                               aCGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTGGCGCCAGAGGCAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAg                                                      <  1:973941/150‑1 (MQ=255)                                                                                                       aagcccccGCCCGGCCCTGCG‑GAAAAAAC‑CTTGGCGCAGGAGAAAAAAAAAAGCCCTGTCCGGAGGAAACCGCCAATGGCAATCATCGTCAGGTCACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTc                                                  <  2:158020/142‑1 (MQ=255)                                                                                                       aagcccccGCCCGGACCGCCC‑TGAGAAAT‑CGTGCCGCCAGAGGCAAAAAAAAGCCCTGTCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGaaa                                             <  2:675199/147‑1 (MQ=255)                                                                                                       aaTGCCCCGCCCGGCCCTGCC‑TGAAAAAT‑CGTGGCGCCAGAGGAAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGaaa                                             <  1:587638/151‑1 (MQ=255)                                                                                                       aaTCCCCCGCCCGGCCCTGCC‑GAAAAAAT‑CGTTGCGCCAGAGGAAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGtt                                                    <  1:922083/144‑1 (MQ=255)                                                                                                        atcccccGCCGGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAAAAAAAAGCCCTGTACCGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAAc                                            <  2:306638/148‑1 (MQ=255)                                                                                                            cccgcccgGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCAt                                         >  1:315476/1‑150 (MQ=255)                                                                                                            cccgcccgGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCAt                                         >  2:553388/1‑150 (MQ=255)                                                                                                             ccgcccgGCCCTGCC‑TGAAAAAT‑CTTTGCGCCAGAGGCAACAAAAAGCCCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATAc                                       <  2:319967/151‑1 (MQ=255)                                                                                                              cgcccgGACCTGCC‑TGAAAAAT‑CGTTGCGCAGGAGGAAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACa                                      <  1:759790/151‑1 (MQ=255)                                                                                                                cccgGACCTGCC‑GGAAAAAT‑CTTTGCGCCAGAGCCAAAAAAAAGCCCTTTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACaaa                                    <  1:1172908/151‑1 (MQ=255)                                                                                                                         gcc‑tgAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATc                           >  1:718140/1‑151 (MQ=255)                                                                                                                          cc‑tgAAAAAT‑CGTTGCGCCAGAGCCAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTAt                            <  1:594534/149‑1 (MQ=255)                                                                                                                               agaAAT‑CGTGGCGCAAGAGACAACAAAAAGCCCTGTCCGGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCAt                        <  1:1075737/149‑1 (MQ=255)                                                                                                                                 aaaT‑CGTTGCGCCAGAGGCAACAAAAAGACCCGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCtt                     >  1:1048715/1‑150 (MQ=255)                                                                                                                                 aaaT‑CGTGGCGCCAGGGGAAAAAAAAAGCCCTGCCCCGAGGGACCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTa                    <  2:722229/151‑1 (MQ=255)                                                                                                                                      cGTTGCGCCAGAGGAAAAAAAAAGCCCTGTCCCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTaaaa                 <  1:82862/150‑1 (MQ=255)                                                                                                                                       gTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAAt               >  1:1061052/1‑151 (MQ=255)                                                                                                                                          gcgcCAGAGCAAAAAAAAAGCCCTATACAGAGGAACCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTaaaaa                <  2:577448/147‑1 (MQ=255)                                                                                                                                           cgcCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAgg           >  1:213770/1‑151 (MQ=255)                                                                                                                                           cgcCAGAGACAAAAAAAAGACCTGTACCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTTAAAAt               <  2:765609/147‑1 (MQ=255)                                                                                                                                             ccAGGGGCAAAAAAAAGCCCTGTCCCGAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGa             <  2:292746/147‑1 (MQ=255)                                                                                                                                               agagGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGcgcg       >  2:1192995/1‑151 (MQ=255)                                                                                                                                               agagGAAAAAAAAAGCCCTGTCCCGATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGgcgc        <  1:1135634/150‑1 (MQ=255)                                                                                                                                                   gCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTg   >  2:1124081/1‑151 (MQ=255)                                                                                                                                                      acaaAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGt  >  1:976354/1‑149 (MQ=255)                                                                                                                                                      acaaAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGt  >  1:978195/1‑149 (MQ=255)                                                                                                                                                       |                                                                                                                                                    CGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGT  >  NC_000913/2248277‑2248571

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