Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,826,379	0	T	C	51.8%	58.5 / 108.9	81	E182G (GAG→GGG)	astE	succinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base T (29/10);  new base C (0/42);  total (29/52)
Fisher's exact test for biased strand distribution p-value = 7.82e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.12e-10

TGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGG‑CAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATA  >  NC_000913/1826231‑1826527                                                                                                                                                           |                                                                                                                                                     tGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGG‑CAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTcc                                                                                                                                                    <  2:512042/151‑1 (MQ=255)      tGGCAGTTACTGCAAACCAGCGCAGATCGTTTTGCCCGAAGGG‑CAACGATTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAA‑CGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGcac                                                                                                                                               <  1:799136/151‑3 (MQ=255)       ggCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGG‑CAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGccc                                                                                                                                               >  1:1042391/1‑150 (MQ=255)       ggCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGG‑CAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGccc                                                                                                                                               >  1:666800/1‑150 (MQ=255)         gaGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGG‑CAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCccg                                                                                                                                             >  1:1308976/2‑150 (MQ=255)                        ggCGAAGATCGTTTTGCCCAAAGGG‑CAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGt                                                                                                                             >  1:539491/1‑151 (MQ=255)                              gATCGTTTTGCCCAAAGGG‑CAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGaaa                                                                                                                       >  2:416753/1‑151 (MQ=255)                                 tgTTTTGGCCAAAGGGCCAACGCTTTGCCAAGTCCCAGCGTACAGCCCAGCGCGCAAAAAGGTCTGGCGCTGAAA‑TGGGAAAC‑CGCCCCCCCAGG‑TTCCTGATGG‑AACC‑CCACCGCCTCCACCCCCGCCGCACCCAGCCACGTCAGAAAt                                                                                                                      <  1:708714/149‑1 (MQ=255)                                  gTTTTGCCCAAAGGG‑CAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAAttt                                                                                                                    >  1:811270/1‑150 (MQ=255)                                   ttttCCCAAAAGGG‑CAACGTTTTGCCAATTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTtct                                                                                                                  <  2:1162459/151‑1 (MQ=255)                                       gCCCAAAGGG‑CAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGTCGCTGAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGt                                                                                                               >  1:88356/1‑150 (MQ=255)                                        ccaaaaGGC‑CACCTCTTTGCAAAGTCCCAGCGTACAGGCCAGCGCCCAAAAAGGTCTGCCGCTAAAA‑TGGGAAAC‑CGTCCCCCCAGT‑TCCCTAAGGA‑AACC‑CCAGCCCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTcc                                                                                                             <  1:64238/148‑1 (MQ=255)                                               ggccagcGCTTTGCAAAGTCCAGGCGAACGGGCCAGCGCGCAAAAATGTCTGGCGTAAAAA‑GGGTTAAG‑CGTCCCCCCAGT‑TTCCTGAGGA‑ACCC‑CCGGCCCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAggg                                                                                                        <  1:573069/144‑1 (MQ=255)                                               gg‑CACCGTTTTGCCAAGTCCCGGCGTACAGCCCCGCGCCCCAAAATGTCGGGCGCTAAAA‑TGGGTAAA‑CGTACCCCCAGG‑TTCCTGATGG‑AACC‑CCAGCCCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGa                                                                                                       <  2:147012/150‑1 (MQ=255)                                               gg‑CAACGCTTTGCCAATTTCCAGCGTACAGGCCAGCGCGCAAAAATGTCTGGCGCTGAAA‑TGGGTAAA‑CGTCCCCCCAGG‑TTCCTGATGA‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGa                                                                                                       <  2:1084850/150‑1 (MQ=255)                                               gg‑CAACCCTTTGCCCAGTTCCAGCGGACAGGCCCGGGCGCCAAAATGGCTGGGCCTGAAA‑TGGGTAAC‑CGTCCCCCCAGG‑TTCCGGAGGA‑AACC‑CCAGCCCCCCCACCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGa                                                                                                       <  2:672556/150‑1 (MQ=255)                                                ggaaaCGTTTTGCCAATTTCCGGCGTACAGGCCAGCGCGAAAAAATGTCTGGCGCTGAAA‑TGGGTAAA‑CGTCCCCCCAGG‑TTCCTGATGG‑AACC‑CCAGCGCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGa                                                                                                       <  2:181078/147‑1 (MQ=255)                                                gccaccGCTTTGCCAAGTCCCAGCGTACGGCCCAGCGCCCAAAAATGTCGGGCGCTGAAA‑TGGTTAAC‑CGTCCCCCCAGG‑TCCCTGAGGA‑AACC‑CCAGCGCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGa                                                                                                       <  1:844660/145‑1 (MQ=255)                                                   aaCCCTTTGCCAAGTTCCAGCGTACGGGCCAGCGCCCAAAAATGTCGGGCGCTGAAA‑GGGGTAAA‑CGTCCCCCCAGG‑TTCCTGAGGG‑AACC‑CCAGCGCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGt                                                                                                   <  2:301915/151‑1 (MQ=255)                                                                   cgggCGTACGGCCCAGCGCCCAAAAATGTCTGGCGCTAAAA‑GGGGAAAC‑CGTCCCACCAGT‑TTCCTGATGG‑AACC‑CCAGCGCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATac                                                                                   <  2:497235/148‑1 (MQ=255)                                                                   ccAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATa                                                                                    >  1:616528/1‑150 (MQ=255)                                                                   ccAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATa                                                                                    >  1:708944/1‑150 (MQ=255)                                                                   ccAGCGTACAGGCCAGCGCCCAAAAATGTCTGCCCCTAAAA‑TGGGTAAC‑CTTCCCCCCATG‑TTCCTGATGG‑AACC‑CCAGCGCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATac                                                                                   <  2:162061/151‑1 (MQ=255)                                                                     cgagtcccggCCCGCGCCCAAAAAGCTCCGCCGCAGAAG‑GGGGCACA‑CGTCCCACCAGG‑TCCCTGAAGG‑ACCC‑CCCGCCCCCCAACCCCCGCCGCCCCCAGCCACGTCAGGAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATacgc                                                                                 <  1:27422/143‑3 (MQ=255)                                                                     accgtccggccccgcgcgcAAAAATGTCGGGCGCTAAAA‑GGGGTAAC‑CTTCCCCCCGGG‑TTCCTGAGGA‑AACA‑CCAGCCCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATacac                                                                                 <  1:588009/138‑1 (MQ=255)                                                                      gcgtccggcccAGCGCCCAAAAGTTTCGCGCGCTAAAG‑TGATTAAC‑AGTCCCCCCAGT‑TTCTTGTTGA‑AACC‑CCCCCCCCCCCAACCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATacac                                                                                 <  1:412587/141‑1 (MQ=255)                                                                      gCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATacac                                                                                 >  1:1077884/1‑150 (MQ=255)                                                                      gCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATac                                                                                   >  1:497043/1‑148 (MQ=255)                                                                       cgaaCAGGCCGGCGCGCAAAAATGTCGGGCGCTAAAA‑TGGGAAAC‑CGCCCCCCCAGG‑TTCCTGATGG‑AACC‑CCAGCGCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATacac                                                                                 <  2:728944/146‑1 (MQ=255)                                                                       cGTACGGGCCGGCGGGCAAAAATGTCTGGCGCTGAAA‑GGGGAAAA‑CGTCCCCCCGGG‑TTCCTGATGG‑AACC‑CCAGCGCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATacac                                                                                 <  2:278496/149‑1 (MQ=255)                                                                       cGTACGGGCCAGCGCGCCAAAGTGTCGGGCGGTAAAA‑TGGTTAAC‑CGTCCCCCCAGT‑TTCCGGATGA‑AACA‑CCGGCCCCCCCGGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATacac                                                                                 <  1:520575/149‑1 (MQ=255)                                                                       cGTACGGGCCAGCGCGCCAAAATGTCTGGCGCTGAAA‑TGGGTAAA‑CGTCCCCCCGGG‑TTCCGGATGG‑AACA‑CCAGCGCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATacac                                                                                 <  2:118000/149‑1 (MQ=255)                                                                       cGTACAGGCCAGCGCGCCAAAATGTCGGGCGCTAAAA‑TGGGTAAA‑CGTCCCCCCAGG‑TTCCTGATGG‑AACC‑CCAGCGCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATacac                                                                                 <  2:600504/149‑1 (MQ=255)                                                                        gTACAGGCCGGCGCGCCAAAATGTCGGGCGCTAAAA‑GGGGTAAC‑CCTCCCCCCAGG‑TTCCTGAGGA‑AACA‑CCAGCGCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATac                                                                                   <  1:1275646/146‑1 (MQ=255)                                                                        gTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATac                                                                                   >  2:1276255/1‑146 (MQ=255)                                                                        gTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGa                                                                              >  2:868432/1‑151 (MQ=255)                                                                        gTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGa                                                                              >  2:1117551/1‑151 (MQ=255)                                                                          acgggCCCGGGCGCAAAAATGTCTGCCGCTGAAT‑GGGGTAAA‑CGTCCCCCCAGT‑TCCTTGAGGG‑AACA‑CCAGCTCTTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAAc                                                                            <  2:691225/148‑1 (MQ=255)                                                                          acgggCCAGCGCGCAAAAATGTCTGGCGCTGAAA‑TGGTTAAA‑CGTCCCCCCAGG‑TTCCTGAGGG‑AACC‑CCAGCGCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAAc                                                                            <  2:66182/148‑1 (MQ=255)                                                                             ggCCAGCGCGCCAAAATGTCTGGCGCTGAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTg                                                                          >  1:1347066/1‑150 (MQ=255)                                                                             ggCCAGCGCGCCAAAATGTCTGGCGCTGAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTg                                                                          >  1:1086759/1‑150 (MQ=255)                                                                              gCCAGCGCGCAAAAAGGTCGGGCGCTGAAA‑TGGGTAAA‑CGTCCCCCCAGG‑TTCCGGAGGG‑AACC‑CCGGCGCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCAGTAATACACCGAAc                                                                            <  2:1411911/147‑1 (MQ=255)                                                                              cccggcgcgcAAAAAGGTCTGGC‑CCGAAA‑GGGGGAAACCGTCCCACCAGT‑TTCCTGAGAG‑AACC‑CCAGCGCCCCCACCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAAc                                                                            <  1:934732/143‑1 (MQ=255)                                                                              cccAGCGCGCAAAAATGTCTGGCGCTAAAA‑GGGGTAAA‑CGTACCACCAGG‑TTCCTGATGA‑AACC‑CCAGCCCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAAc                                                                            <  2:1279903/146‑1 (MQ=255)                                                                               ccggggggcaaaaaTGTTGGGCGTTGAAA‑GGGGAAAC‑CGGCCCACAGGG‑TTCCTGATGA‑AACA‑CCGGCCCCCCCGGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCgg                                                                       <  2:224141/141‑1 (MQ=255)                                                                               ccAGCGCGCAAAAAGGTCTGGCGCTGAAA‑TGGGTAAA‑CGTCCCACCGGT‑TTCCTGATGG‑AACA‑CCGGCGCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCgg                                                                       <  2:660689/151‑1 (MQ=255)                                                                               ccAGCGCCCAAAAATGTC‑GGCCGTGAAAGTGGTAAAA‑CGTCCCCCAGGG‑TTCCTGATGG‑AACC‑CCAGCCCCCCCACCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCgg                                                                       <  2:900374/151‑1 (MQ=255)                                                                                       caaaaattgcgggggcTAAAA‑GGGGAAAA‑CGCACCACCGGG‑TTCGTGATGG‑AACC‑CCGGCCCCCCCGGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAgg                                                               <  2:1441597/137‑1 (MQ=255)                                                                                           attttgtggctataaaa‑GGGGAAAA‑CATCCCACCAGT‑TTCCTGAAGA‑AACC‑CCAGCCCCTCCTGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGcc                                                           <  2:335156/137‑1 (MQ=255)                                                                                              gTCTGGCGCTGAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGCTGCGGGAATACACCGAACTGCGGATGCAAGGAGCCACg                                                        >  1:1117776/1‑151 (MQ=255)                                                                                               tcggccccaaaaa‑TGGATAAA‑CGTCCCCCCAGG‑TCCCTGAGGG‑AACC‑CCGGCGCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACg                                                        <  2:678806/140‑1 (MQ=255)                                                                                                      cTGAAA‑TGGGAAAA‑CGTCCCCCCGGG‑TTCCGGAGGA‑AACC‑CCGGCGCCCCCACCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCgg                                                <  2:539723/151‑1 (MQ=255)                                                                                                        gAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGTATTGAGCCaattactgaat                                                >  2:1385344/1‑139 (MQ=255)                                                                                                        gAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGt                                               >  2:1010728/1‑150 (MQ=255)                                                                                                           a‑TGGATAAA‑CGACCCCCCGGG‑TTCCTGATGA‑AACC‑CCGGCCCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGt                                           <  1:1100491/151‑1 (MQ=255)                                                                                                            atgagaaaa‑CGTACCACC‑GG‑ATACTAATGGCAAAC‑CCAGCGCCCCTTGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGt                                           <  2:944582/145‑1 (MQ=255)                                                                                                              gggTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAg                                         >  2:691817/1‑151 (MQ=255)                                                                                                              gggTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAg                                         >  2:293516/1‑151 (MQ=255)                                                                                                              gggTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAg                                         >  2:502450/1‑151 (MQ=255)                                                                                                              gggTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGAGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAg                                         >  2:1033854/1‑151 (MQ=255)                                                                                                              gggAAAC‑CGTCCCACCAGG‑TCCCTGAGGG‑AACC‑CCAGCCCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGt                                           <  2:311629/149‑1 (MQ=255)                                                                                                              gggAAAA‑CCTCCCCCCGGG‑TTCCTGAGGG‑AACC‑CCGGCCCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGt                                           <  2:706167/149‑1 (MQ=255)                                                                                                                 tAAC‑CGTCCCCCCGGG‑TCCCTGAGGA‑ACCC‑CCAGCGCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATc                                      <  2:1192577/151‑1 (MQ=255)                                                                                                                 tAAA‑CCTACCGCCAGT‑GCCCTTAGGA‑AACA‑CGGCCGCCCCCGCCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGAt                                       <  1:913527/150‑1 (MQ=255)                                                                                                                 aaaa‑CTTCCCCCCGGG‑TTCTTGAGGG‑AACA‑CCCCCGCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATc                                      <  2:1090663/150‑1 (MQ=255)                                                                                                                   aaaccgtcccccaGGTTTCTTGAGGG‑AACC‑CAAGCGCCCCCGGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATc                                      <  1:238671/141‑1 (MQ=255)                                                                                                                       gTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGt                                 >  2:1483659/1‑151 (MQ=255)                                                                                                                        ccccaccgGT‑TTCCTGA‑GG‑AACACCCGGCCCCCCCCGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTg                                <  1:70529/149‑1 (MQ=255)                                                                                                                          cccccggg‑TTCCTGAGGA‑AACC‑CCAGCCCCTCCACCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGcc                              <  2:1309628/145‑1 (MQ=255)                                                                                                                            cccggg‑TTCCGGAGGG‑ACCC‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAg                            <  2:757100/147‑1 (MQ=255)                                                                                                                                 g‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGCATGTAGATCAAGGTGCCAGCGca                        >  1:710974/1‑150 (MQ=255)                                                                                                                                  tttCGTGAGGC‑AACA‑CCGGCCCCCCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGcac                       <  1:1498195/150‑1 (MQ=255)                                                                                                                                   ttCCTAAGGA‑AACC‑CCGGCGCCCCCGCCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGcac                       <  1:986449/150‑1 (MQ=255)                                                                                                                                      ctaaTGA‑AACA‑CCACCCCCCCAAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGAt                   <  1:101184/148‑1 (MQ=255)                                                                                                                                         aTGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGCAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATc                                      >  2:1208985/1‑129 (MQ=255)                                                                                                                                          tGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGCATGTAGATCAAGGTGCCAGCGCACAGATTCtt               >  1:749080/1‑151 (MQ=255)                                                                                                                                           tggaacccccggcgcCCCCGGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTc                 <  1:1249805/139‑1 (MQ=255)                                                                                                                                              aaacccaggcgcCCCCCGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTg             <  1:116774/143‑1 (MQ=255)                                                                                                                                              aaCA‑CCGGCCCCCCCGCCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGc            <  2:141758/151‑1 (MQ=255)                                                                                                                                               aca‑cCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGc            >  1:827901/1‑150 (MQ=255)                                                                                                                                                         ccTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATa  >  2:200132/1‑151 (MQ=255)                                                                                                                                                           |                                                                                                                                                     TGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGG‑CAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAA‑TGGGTAAA‑CGTACCACCAGG‑TTCCTGATGG‑AACA‑CCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCGAACTGCGGATGCAAGGAGCCACGAATTGCGGTATGTAGATCAAGGTGCCAGCGCACAGATTCTTTGCCCTGGTCATA  >  NC_000913/1826231‑1826527

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