Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,248,437	0	A	C	43.8%	96.6 / 20.1	65	G32G (GGT→GGG)	lysP	lysine transporter
Reads supporting (aligned to +/- strand):  ref base A (26/10);  new base C (0/28);  total (26/39)
Fisher's exact test for biased strand distribution p-value = 4.44e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCA  >  NC_000913/2248291‑2248586                                                                                                                                                      |                                                                                                                                                       gAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑g                                                                                                                                                    >  2:1182779/1‑150 (MQ=255)   aCCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATg                                                                                                                                                 <  2:1065673/151‑1 (MQ=255)    ccGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATgg                                                                                                                                                <  1:1012264/151‑1 (MQ=255)     cGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATgg                                                                                                                                                >  2:110463/1‑150 (MQ=255)       gAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAAc                                                                                                                                             <  1:1036376/151‑1 (MQ=255)        aaaCCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAAcc                                                                                                                                            <  1:763430/151‑1 (MQ=255)                atatatGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATgg                                                                                                                                    >  1:984882/1‑151 (MQ=255)                   tatGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCa                                                                                                                                  >  1:69127/1‑150 (MQ=255)                             gtccacggggctgggTCACAAGAAAGTAACCCACCAGGCAAACCAGCAAATACGAGAGCAATCCCCCGCCCGGCCCTGCCG‑GAAAAATCGTG‑GCGCAAGGGGAAACAAAAAGCCC‑TGTCCC‑GATGGAACCGCCAATGGCAatca                                                                                                                              <  1:1118139/131‑1 (MQ=255)                             gTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTc                                                                                                                         >  2:905154/1‑149 (MQ=255)                                   cGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGAcg                                                                                                                 >  1:86653/1‑151 (MQ=255)                                     agaCTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGAcgcg                                                                                                               >  1:143575/1‑151 (MQ=255)                                          ggTCATCAGAAAGAAAACCATCAGGCAAATCAGCATATACGAGAGCAACCCCCCGCCCGGCCCTGCCT‑GAAAATTCGTG‑GCGCAAGAGGCAAAAAAAAGCCC‑TGTCCC‑GATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCttt                                                                                                          <  2:436299/151‑1 (MQ=255)                                           gTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCttt                                                                                                          >  2:189512/1‑150 (MQ=255)                                              atcaGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAAtt                                                                                                      >  1:1031161/1‑151 (MQ=255)                                              acaagaaaGAAACCCACCAGCCCAACCAGCATATACGAGACAAATCCCCCCCCCGGCCCTGCCT‑GAAAAATCGTG‑GCGCAAGAGGCAAAAAAAAGCCC‑TGTACC‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTa                                                                                                         <  1:65591/142‑1 (MQ=255)                                                  ggAAGAAAACCATCAGCCAAATCAGCATATACGAGAGAAACCCCCCGCCCGGCCCTGCCT‑GAAAAATCGTT‑GCGCAAGAGGAAAAAAAAAGCCC‑TGTCCC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCa                                                                                                    <  2:629611/149‑1 (MQ=255)                                                    aaacaacccaccaGGCAAATCACAATAAACGAGACAAACCCCCCGCCCGCCCCTGCGG‑AAAAAATCGTG‑GCGCGAGAGGCAAAAAAAAGCCC‑GGTACC‑AATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACgg                                                                                                 <  1:824694/143‑1 (MQ=255)                                                        aaacctgacgccaaaTCAGAATAACCGAGAGCAACCCCCCGCCCGGCCCTGCCT‑GAAAAATCTTT‑CCGCGAGAGCCAACAAAAAGCCC‑TATCCC‑GATGGATCCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTa                                                                                             <  2:570819/137‑1 (MQ=255)                                                        aaaCCACCAGCCAAATCAGCATATACGAGAGAAACCCCCCCCCGGGCCTTGCCT‑GAAAAACCGTG‑GCGCCA‑AGGAAAACAAAAGCCCTTGTCCC‑GAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGt                                                                                              <  2:238888/149‑1 (MQ=255)                                                                 gccaaaTCAGCATATAGGAGAGCATCCCCCCGCCGGGCCCTGCCT‑GAGAATTGTTG‑GGGCCGGAGGCAAAAAAAAGCCC‑GGCCCT‑GATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGcc                                                                                         <  2:1129546/141‑1 (MQ=255)                                                                 gaccaaaaagaatatACGAGAGCAATCCCCCGCCGGGACC‑GCCTGGAGAATTGTTG‑GCCCAAGGGGCAAAAAAAAGCCC‑GTCCCC‑GAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGgc                                                                                     <  2:1089610/138‑1 (MQ=255)                                                                  gCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCt                                                                                  >  2:124073/1‑151 (MQ=255)                                                                   ccAATCAGAATATACGAGAGCAATCCCCCCCCCGGCCCGGCCT‑GAAAAATCGTT‑GCGCAAGAGCCAAAAAAAAGACC‑TGTACC‑GAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCt                                                                                  <  1:1142987/150‑1 (MQ=255)                                                                     aaccagataataCGAGAGCAATCCCCCCCCCGGCCCTGCCT‑GAAAATTGGTT‑GCCCCAGAGGAAAAAAAAGGCCC‑TGTCCC‑GATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCt                                                                               <  2:1164288/142‑1 (MQ=255)                                                                     aaccAGCATATGGGAGAGCAACCCCCCCCCCGGCCCTGCCG‑GAGAAATCGTG‑GCGCAAGGGGAAAAAAAAAGCCC‑TTTCCC‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCtt                                                                                 <  2:451509/146‑1 (MQ=255)                                                                     aaccAGAATATACGAGAGAAATCCCCCCCCGGGCCCTCCTT‑GAAAAATCGTT‑GCGCGAGAGCCAAAAAAAAGCCC‑TGTCCC‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCt                                                                               <  2:361873/148‑1 (MQ=255)                                                                     aaTAAGCATATAGGAAAGCAGCCCCCCGCCCGCACCAGCGT‑GAAAATTGGTG‑GGGCAAGAGAAAAAAAAAAGACC‑TGCCCA‑GAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTc                                                                                <  1:521892/150‑1 (MQ=255)                                                                        cAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtgtg                                                                            >  2:331125/1‑151 (MQ=255)                                                                             atagcaaagagcacgCCCCCGCCGGGCCTGGCC‑GAAAAATTGTTGGCGCCAGAGAAAAAAAAAACCCT‑TCTCCC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCtgt                                                                             <  2:1158819/132‑1 (MQ=255)                                                                                aCGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGt                                                                    >  2:818207/1‑151 (MQ=255)                                                                                 cGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGtt                                                                   >  2:785100/1‑151 (MQ=255)                                                                                       cacccccccGCCCGGACCTGCCG‑GAGAAATGTCG‑GCGCAAGGGGCAA‑AAAAAGGCC‑CGTTCCGGAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGtt                                                                         <  1:31713/134‑1 (MQ=255)                                                                                       caaccccccGCCGGGCCCTGCCA‑AAAAAATCGTT‑GCCCAAGAGCAAAAAAAAAGCCC‑TGTCCA‑GAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCg                                                                <  1:233274/143‑1 (MQ=255)                                                                                        aaTGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGaa                                                             >  1:963138/1‑150 (MQ=255)                                                                                         accccccGCCCGCCCCTGCCT‑GAAAAATCGTC‑GCGCCAGAGGAAAAAAAAAGCCC‑TTTCCC‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGttt                                                                  <  1:547175/141‑1 (MQ=255)                                                                                           gCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAAcc                                                          >  1:1107713/1‑150 (MQ=255)                                                                                             cccccccgCCCCTGCCG‑GAAAAATCGTT‑GCGCCAGAGCCAAAAAAAAGCCC‑TGTCCC‑GATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATa                                                       <  2:416463/151‑1 (MQ=255)                                                                                                 cccgcccctgcct‑aAAAAATGGTG‑CCGCCGAAGCCAAAAAAAAGCCC‑TGCCCC‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACa                                                     <  1:1057308/143‑1 (MQ=255)                                                                                                 cccgGCCCTGCCT‑GAAAAATCGTT‑GCGCCAGAGGAAAAAAAAAGCCC‑TGTCCC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACaa                                                    <  2:219148/150‑1 (MQ=255)                                                                                                           cct‑gAAAAACCTTG‑GCGCCAGAGGAAAAAAAAAGCCC‑TGTCCC‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACaaaaa                                                 <  1:992253/143‑1 (MQ=255)                                                                                                            ct‑gAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATcc                                         >  2:145797/1‑150 (MQ=255)                                                                                                               gagaAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATc                                      >  1:408047/1‑151 (MQ=255)                                                                                                               gaaaaaTCTTT‑GCACCAGAGCCAAAAAAAAGCCC‑TGTACC‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATc                                          <  2:567710/144‑1 (MQ=255)                                                                                                                aaaaaTCGTT‑GCGCCAGAGGCAAAAAAAAGCCC‑TGTCCC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCAt                                       <  2:736901/147‑1 (MQ=255)                                                                                                                aaaaaTCGTT‑GCGCCAGAGAAAAAAAAAAGCCC‑TGTACC‑GAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATcc                                         <  1:1016826/145‑1 (MQ=255)                                                                                                                aaaaaCCGTG‑GCGCCAGAGGAAAAAAAAAGCCC‑GGTCCC‑GAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCa                                        <  1:238714/146‑1 (MQ=255)                                                                                                                 gaAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCtt                                    >  1:787269/1‑151 (MQ=255)                                                                                                                  aaaGCGTT‑GCGCCAGAGCAAAAAAAAAGACC‑TGTCCC‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCAtat                                     <  1:1141827/149‑3 (MQ=255)                                                                                                                  aaaGCGTT‑GCGCCAGAGAAAACAAAAAGCCC‑TGTCCA‑GATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCtt                                    <  2:193385/150‑1 (MQ=255)                                                                                                                   aaTCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTa                                   >  1:953537/1‑150 (MQ=255)                                                                                                                        tt‑GCGCAAGAGGCAAAAAAAAGCCC‑TGTCCC‑GAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTaaaa                                <  2:63661/148‑1 (MQ=255)                                                                                                                         tggcgcCAGAGGAAACAAAAAGCCC‑TGTCCC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAAt                              <  1:230280/148‑1 (MQ=255)                                                                                                                            cgcCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAg                           >  2:207296/1‑150 (MQ=255)                                                                                                                              ccAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGgcg                        >  2:455553/1‑151 (MQ=255)                                                                                                                              ccAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGgcg                        >  1:970646/1‑151 (MQ=255)                                                                                                                                     caacaaAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGt                 >  1:551697/1‑151 (MQ=255)                                                                                                                                      aacaaAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGt                 >  1:18317/1‑150 (MQ=255)                                                                                                                                      aaaaaaaaGCCC‑TGTCCC‑GATGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGAtt                   <  2:692306/145‑1 (MQ=255)                                                                                                                                       aaaaaaaGCCC‑TGCCCC‑GAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGgc                         <  1:1037363/139‑1 (MQ=255)                                                                                                                                          aaaaGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGc              >  1:1177791/1‑149 (MQ=255)                                                                                                                                                cc‑TGTCCG‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACg         <  1:1045051/148‑1 (MQ=255)                                                                                                                                                cc‑TGTCCC‑GAGGAAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGAtt      <  1:909311/151‑1 (MQ=255)                                                                                                                                                    gTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTg     >  1:341170/1‑149 (MQ=255)                                                                                                                                                     tACC‑GAGGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTg     <  1:89345/148‑1 (MQ=255)                                                                                                                                                      aCC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCa  >  2:829239/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                       GAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGCAGACGATTGGCA  >  NC_000913/2248291‑2248586

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