Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	650,464	0	A	T	41.4%	80.3 / 25.6	58	R9R (CGT→CGA)	citE	citrate lyase, citryl‑ACP lyase (beta) subunit
Reads supporting (aligned to +/- strand):  ref base A (8/26);  new base T (24/0);  total (32/26)
Fisher's exact test for biased strand distribution p-value = 1.40e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGC  >  NC_000913/650314‑650610                                                                                                                                                        |                                                                                                                                                    aaCCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTa                                                                                                                                                     <  2:779494/151‑1 (MQ=255) aaCCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTa                                                                                                                                                     <  1:206182/151‑1 (MQ=255)         ggcgggcgGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCCCCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGttgttg                                                                                                                                             <  1:161450/151‑1 (MQ=255)          gcgggcggTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGc                                                                                                                                            <  1:181898/151‑1 (MQ=255)               cggTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGA‑GGtttttgtttttgga                                                                                                                                           >  2:963208/1‑145 (MQ=255)                ggTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGaa                                                                                                                                      <  1:214841/151‑1 (MQ=255)                  tgtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTTTTGCAGCGaa                                                                                                                                      >  1:347884/1‑149 (MQ=255)                  tgtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGG‑GGTTTTCgttttt                                                                                                                                              >  2:438134/1‑141 (MQ=255)                   gtCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCg                                                                                                                                   <  2:30366/151‑1 (MQ=255)                       ttttCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGaa                                                                                                                                <  1:366069/150‑1 (MQ=255)                        tttCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGaaa                                                                                                                               <  1:195191/150‑1 (MQ=255)                          tCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATc                                                                                                                             <  1:697027/150‑1 (MQ=255)                           cACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑t                                                                                                                          <  2:444487/151‑1 (MQ=255)                           cACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑t                                                                                                                          <  2:220375/151‑1 (MQ=255)                            aCGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGG‑GGTTTTTGTTTTTGGAGCGAAAAGGAAATCaat                                                                                                                          >  2:144694/1‑149 (MQ=255)                               cAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGGGG‑GGTTTTCGTTGTTGCAGCGGAACGGGAAACA‑tt                                                                                                                         >  2:985723/1‑148 (MQ=255)                               cAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGG‑GGTTTTTGTTGTTTGAGCGAAGCGGGAAACA‑tttg                                                                                                                       >  1:458449/1‑148 (MQ=255)                                 aTGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGTTGCGG‑GGTTTACGTTGTTTGAGCGAAACGGAAATCA‑tttg                                                                                                                       >  1:209594/1‑146 (MQ=255)                                  tGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGGaa                                                                                                                    >  1:494497/1‑150 (MQ=255)                                     tACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGGGG‑GTTTTTCGTTGTTGCCACGCAGCGGGAAACC‑TTGGCGAtcc                                                                                                                 >  1:794005/1‑150 (MQ=255)                                     tACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGG‑GTTTTTTGGTTTTGCAGCGAAGCGGGAATCA‑TTTGCAAtcc                                                                                                                 >  2:647538/1‑150 (MQ=255)                                       cGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAAtcctcc                                                                                                              <  1:508858/151‑1 (MQ=255)                                          aGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCAt                                                                                                           <  2:1383/151‑1 (MQ=255)                                             cTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGG‑GTTTTTCGTTGTTGCAGCGAAGCGGAAAACA‑TTGGCAATCCTccc                                                                                                             >  1:362126/1‑146 (MQ=255)                                                cGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCaga                                                                                                     <  2:911256/151‑1 (MQ=255)                                                cGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGA‑GGTTTTCGTTGTTGCAGCGAAGCGGGAATCA‑TTGGGAAACCTCCCCTg                                                                                                          >  1:98760/1‑146 (MQ=255)                                                  aGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGc                                                                                                   <  1:738583/151‑1 (MQ=255)                                                      cAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCgg                                                                                                <  2:638845/150‑1 (MQ=255)                                                       aaaCATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGa                                                                                              <  1:910399/151‑1 (MQ=255)                                                       aaaCATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGG‑GGTTTTTGTTGTTTGAGCGGAGCGGGAAACA‑TTGGGAATCCTcccctagc                                                                                                        >  1:944053/1‑136 (MQ=255)                                                                ggCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTgg                                                                                     <  2:165568/151‑1 (MQ=255)                                                                  cATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGG‑GGTTTTCGTTGTTGCAGCGAAGCGGAAATCATTTGGCAA‑CCTCCCATGGCAGAGCCGGGGTACCGCTgg                                                                                     >  1:112725/1‑149 (MQ=255)                                                                     cAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGcgcg                                                                                 <  1:704749/150‑1 (MQ=255)                                                                     cAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGG‑GGTTTTTGTTTTTGGAGGGGAACGGGAATAA‑ATTGCAAACCTCCCCTGGCAGAGGCGGGGTACCGCCgg                                                                                     >  2:263674/1‑146 (MQ=255)                                                                     cAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGG‑GGTTTTCGTTGTTGCCGCGAAACGGAAATCA‑TTGGCAATCCTCCCCTGGCAGAGCCGGGGTACCCCTGgcg                                                                                   >  1:697742/1‑148 (MQ=255)                                                                      aGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGG‑GATTTATGCTGTTGCAGCGAAACGGGAAACA‑TTGGCAATCCTCCCCTGGCAGAGCCGGGGGACCGCTGgcg                                                                                   >  1:787808/1‑147 (MQ=255)                                                                         cGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGTGGACC‑GTTTTTCGGTTATGCAGCGAGGGGGGAATCC‑GACGCCATGCTCCCATGCCTGAGGCggggctcccccg                                                                                      >  1:885035/1‑132 (MQ=255)                                                                           ggTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGG‑GGTTTTTGTTGTTGGAGGGGAGCGGAAAACA‑ATGGGCATACTCCCCCGGCAGAGCCGGGGGACCGCCGGCGCGTggc                                                                             >  2:94465/1‑146 (MQ=255)                                                                            gTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcagca                                                                         <  1:279520/151‑1 (MQ=255)                                                                              aGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGG‑GGTTTTCGTTGTTGCAGCGAAGCGGGAAACA‑ATTGCAATCCTCCCATGGCAGAGCCGGGGTAACGCTGGCGCGGGGCAGCAgg                                                                       >  1:667351/1‑151 (MQ=255)                                                                                     ggAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGATGCGG‑GGTTTACGTTTTTGCAGGGAAGCGGGAATCA‑TTGGCCATCCTCCCCTGGCAGAGCCGGGGGACCGCCGGCGCGGGCCAGCAGGGCtt                                                                   >  1:132346/1‑148 (MQ=255)                                                                                        gTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTc                                                             <  1:968947/151‑1 (MQ=255)                                                                                        gTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTc                                                             <  1:391399/151‑1 (MQ=255)                                                                                          tGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGGGG‑GTTTTTCGTTTTTTCAGCGAAACGGGAAACA‑TTGGCAATCCTCCCACGGCACGGCCGGGGTACCGCCGGCGCGGGCCAGCCGGGGTTccc                                                                >  2:662605/1‑145 (MQ=255)                                                                                          tGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGG‑GTTTTTCGTTGTTTCAGCGAAGCGGGAAACA‑TTTGCAAACCTCCCCTTGCAGAGCCGGGGTACCGCTGGCGCGCGCCCGCAGGGCTTcccctc                                                             >  1:817417/1‑145 (MQ=255)                                                                                           gCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTg                                                          <  1:849277/151‑1 (MQ=255)                                                                                              gACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCAc                                                       <  1:1033259/151‑1 (MQ=255)                                                                                                 cATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGA‑GGTTTTCGTTGTTTCAGCGAAGCGGAAATCA‑TTGGCAAACCTCCCATGGCAGAGCCGGGGTAACGCCGGCGCGGGGCCGCAGGGCTTCCAGTCGTGCAc                                                       >  1:914083/1‑148 (MQ=255)                                                                                                  aTCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGGGAGGTTTTTCG‑TGTTTCGGCGAAACGGGAAACA‑TTGGGAATCCTCCCATGGCCGGGCCGGGGTACCGCCGGCGCGGGCCGCCAGGGGCTCCCCTCGTGcgcgga                                                    >  2:510931/1‑148 (MQ=255)                                                                                                    cgcgGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGG‑GGTTTTCGTTGTTGCAGAGAAGCGGGAAACA‑TTGGCAAACCTCCCCTGGCAGAGCCCGGGGACCCCTGGGGCGGGCCAGCAGGGGTTCCCGTCGTGCCCGTaaa                                                  >  1:766289/1‑150 (MQ=255)                                                                                                      cgGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGG‑GGTTTTTGTTTTTGTAGCGAAAAGGGAAACA‑TTTGCAAACCCCCCCCGGCAGAGCCGGGGTAACCCTGGGGCGGGGCCGCAGGGCGTCCCGGCGTGCACGGaaaaa                                                >  2:155319/1‑149 (MQ=255)                                                                                                            tGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGG‑GTTTTTTGTTGTTGCAGCGAAAAGGGAATCA‑TTTGCAATCCTCCCATGGCAGAGCCGGGGGACCGCTGGCGCGTGCCCGCAGGGCTTCCCGTCCTGCACGCAAAAAGCAGGcc                                         >  2:129780/1‑151 (MQ=255)                                                                                                            tGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGG‑GGTTTTCGTTGTTTCCGCGGAGCGGAAATCA‑TTGGGAATCCTCCCATGGCAGAGCCGGGGGACCGCTGGCGCGTGGCAGCCGGGCGTCCCGTCGTGCCCGGAAAAAGCAGGcc                                         >  1:787175/1‑151 (MQ=255)                                                                                                             ggCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTcc                                         <  1:478170/150‑1 (MQ=255)                                                                                                               cACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGt                                      <  1:279244/151‑1 (MQ=255)                                                                                                                   aGGCACAAACAACATGCTGCGGCGGGGGCGG‑GTTTTTCGTTTTTGCAGCGAAACGGAAAACA‑TTGGCAATCCTCCCCTGGCAGAGCCGGGGTACCGCTGGCGCGGGCCAGCAGGGGCTCCCGTCGTGCACGTAAAAAGCAGTCCAGGgcg                                   >  1:288802/1‑150 (MQ=255)                                                                                                                                aTGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGAc                     <  2:112698/151‑1 (MQ=255)                                                                                                                                      cggcggGGGCGA‑GTTTTACGTTTTTGCAGCGAAGCGGGAATCA‑TTGGCAATCCTCCCCTGGCAGAGCCGGGATACCGCTGGCGCGGGCCAGCAGGGCTTCCAGTCGGGCACGTAAAAAGCAGTCCAGTGCGCCTTTGTCATCGACCtt                  >  2:308412/1‑148 (MQ=255)                                                                                                                                       ggcggGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAg               >  2:26007/1‑150 (MQ=255)                                                                                                                                                 a‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTAcgccgc    <  1:1034103/151‑1 (MQ=255)                                                                                                                                                     tttACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCgcgc  >  1:1022336/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                    AACCATGCGGCGGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGA‑GTTTTACGTTGTTGCAGCGAAGCGGAAATCA‑TTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTAAAATGCAGTCCAGTGCGCCTTTGTCATCGACATTCAGCTGTACGCCGCGC  >  NC_000913/650314‑650610

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