Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,248,424	0	C	A	41.7%	90.4 / 34.5	61	V37F (GTT→TTT)	lysP	lysine transporter
Reads supporting (aligned to +/- strand):  ref base C (23/12);  new base A (0/25);  total (23/38)
Fisher's exact test for biased strand distribution p-value = 5.09e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGT  >  NC_000913/2248275‑2248571                                                                                                                                                        |                                                                                                                                                      aCCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGcaac                                                                                                                                                       >  2:1142994/1‑150 (MQ=255)       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         >  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)                                                                                                                                                        |                                                                                                                                                      ACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCACCGAGACTGGTCATCAGGAAGTAAACCATCAGGCCAATCAGCATATACGAGAGCAATGCCCCGCCCGGACCTGCCT‑GAGAAATCGTT‑GCGCCAGAGGCAACAAAAAGACC‑TGTACC‑GATGGAACCGCCAATGGCAATCATCGTCAGGTGACGCGCCTTTAATTCACGGCGTAAGCCCGGCGCTTCTGTGGTTTTAGTTTCGGAAACCATACAAAAATGCTATCCATCTTAAAAATGAGGCGCGATTGT  >  NC_000913/2248275‑2248571

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