Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,248,431	0	A	C	49.2%	73.4 / 37.0	67	G34G (GGT→GGG)	lysP	lysine transporter
Reads supporting (aligned to +/- strand):  ref base A (26/8);  new base C (0/33);  total (26/41)
Fisher's exact test for biased strand distribution p-value = 8.30e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGAC  >  NC_000913/2248284‑2248578                                                                                                                                                       |                                                                                                                                                    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                                                >  2:553388/1‑150 (MQ=255)                                                                                                     cccgcccgGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCAt                                                >  1:315476/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:978195/1‑149 (MQ=255)                                                                                                                                               acaaAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGt         >  1:976354/1‑149 (MQ=255)                                                                                                                                                    aaGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGAc  >  2:573045/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                    GGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACC‑ATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCC‑TGAGAAAT‑CGTTGCGCCAGAGGCAACAAAAAGACCTGTACCGATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGAC  >  NC_000913/2248284‑2248578

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