Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,248,431	0	A	C	43.8%	97.8 / 23.4	65	G34G (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

AGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGC  >  NC_000913/2248281‑2248574                                                                                                                                                         |                                                                                                                                                  aGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGa                                                                                                                                                   <  1:430875/151‑1 (MQ=255)    ggCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGAcc                                                                                                                                                 <  2:1088514/150‑1 (MQ=255)     gCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑Tg                                                                                                                                              <  1:1039560/151‑1 (MQ=255)     gCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑Tg                                                                                                                                              >  1:792925/1‑151 (MQ=255)      cAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGt                                                                                                                                             <  1:56149/151‑1 (MQ=255)           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            >  1:970646/1‑151 (MQ=255)                                                                                                                                        ccAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGgcg            >  2:455553/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)                                                                                                                                                         |                                                                                                                                                  AGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGTAGC  >  NC_000913/2248281‑2248574

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