Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,480,371	0	T	G	54.0%	68.3 / 94.8	63	G390G (GGT→GGG)	dsdA	D‑serine dehydratase
Reads supporting (aligned to +/- strand):  ref base T (6/23);  new base G (34/0);  total (40/23)
Fisher's exact test for biased strand distribution p-value = 5.05e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-09

GTTGGTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGC‑ACGGT‑TTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGATGAATCAATATCTGGC  >  NC_000913/2480225‑2480517                                                                                                                                                    |                                                                                                                                                     gTTGGTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTATgg                                                                                                                                                  <  1:160182/151‑1 (MQ=255)        gcgcATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGGTGGGCgg                                                                                                                                              >  2:496657/1‑148 (MQ=255)              cAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGATTGGCGGGGCACAgcg                                                                                                                                     >  2:159717/1‑151 (MQ=255)               aGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGATGGGCgggac                                                                                                                                           >  2:972908/1‑141 (MQ=255)                 gCTTTGTCGGGCGGGCAATGGAGCGTCTTCTGGATGGCTTCTATACCCTTAGCGATCAAGCCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGCAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTATGGCCTGACCTCatc                                                                                                                                      <  2:26577/147‑3 (MQ=255)                   tttGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGATGGGCGGGCCTCAGcgcg                                                                                                                                   >  2:732231/1‑148 (MQ=255)                          ggcgggcAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGATGGGCGGGGCTCCGCGCGggtgtg                                                                                                                             >  1:347967/1‑147 (MQ=255)                                  aTGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTATGGGCGGGCCTCAGCGCGGGTGTGCATCAGTAAggt                                                                                                                 >  1:496273/1‑149 (MQ=255)                                  aTGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGtt                                                                                                                 <  1:27373/151‑1 (MQ=255)                                   tGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGGACCTTCGGCACTGGC‑GGGGGTTGGCGGGCCTCCAGGCGGGTGTGGCTCAGTaa                                                                                                                    >  2:739249/1‑147 (MQ=255)                                      aGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCa                                                                                                             <  1:189712/151‑1 (MQ=255)                                      aGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGATGGGCGGGCCCCCACGCGGGGGTGCATCAGTAAGTTatc                                                                                                              >  1:751878/1‑148 (MQ=255)                                       gCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGGGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCaa                                                                                                            <  1:1071414/151‑1 (MQ=255)                                         gTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACa                                                                                                          <  1:115766/151‑1 (MQ=255)                                         gTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGATGGGCGGGCCTCAGCGCGGGTGGGCATCAGTAAGGTACCaa                                                                                                            >  1:244813/1‑149 (MQ=255)                                          tCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGATGGGCGGGCCTCCGCGCGGGTGTGCATCAGTAAGTTAacaaca                                                                                                          >  1:219267/1‑150 (MQ=255)                                            tgctgGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGATGGTCGGGCCTCCGCCCGGGTGGGGCGCAGGAAGTTaac                                                                                                              >  2:1074456/1‑142 (MQ=255)                                              ctgGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGATGGGCGGGGCTCAGCGCGGGTGGGGATCAGTAAGTTAACAACAGaag                                                                                                      >  1:351753/1‑148 (MQ=255)                                                 gATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGC‑ACg                                                                                                 <  2:586709/151‑1 (MQ=255)                                                  aTGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGATGGCCGGGACTCAGCGCCGGTGGGCCTCAGTAAAGTACCCAAAGAAGc                                                                                                     >  2:238236/1‑147 (MQ=255)                                                   tGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGGTGGGCGGGCCTCCGCGCGGGAGTGCTTCACTAAGTTTACAAAAGATgg                                                                                                     >  2:873653/1‑145 (MQ=255)                                                        tCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGGTGGGCGGGCCCCCGAGCGGGTGTGCATCAGGAAGTTACCAACAAAAGC‑CCGGg‑tt                                                                                            >  2:63179/1‑148 (MQ=255)                                                        tCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGATGGCCGGACCTCCACGCGTGTGTGCATCAGTAAATTACCAACAAAAGC‑ACGGg‑tt                                                                                            >  1:92315/1‑148 (MQ=255)                                                           ataCCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGC‑ACGGT‑TTCAgc                                                                                        <  2:1155956/149‑1 (MQ=255)                                                           ataCCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGATGGGGGGGCCTCAGCGCGGGTGTGCATTCGGAAGGTAACAACAGATGC‑ACGGg‑ttttgc                                                                                        >  2:876466/1‑145 (MQ=255)                                                              cccTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGCTGGCGGGGATTCACCGCGGGCGTGGATCCGCTAGGTACCAACACATGCAACGGt‑tctgtg                                                                                        >  2:1172999/1‑142 (MQ=255)                                                                  tAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGATGGGCGGGCCTCCGGGCGGGGGGGTCTCCAGAAGTTACCAACAAAAGG‑ACGGG‑TTCAGCGCCGaa                                                                                  >  2:40337/1‑148 (MQ=255)                                                                   aGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGC‑ACGGT‑TTCAGCGCAGaacaa                                                                               <  1:718254/150‑1 (MQ=255)                                                                   aGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGATGGGCGGACCTCCGCGCGGGTGTGCATCAGGAAAGTAACAACAGAAGC‑ACGGG‑TTTAGCGCCGaacaa                                                                               >  1:314922/1‑150 (MQ=255)                                                                    gCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTATGGCCGGGCCTCCCCGCGTGTGTGCCTCAGTAAGTTAACAACAGATGC‑ACGGG‑TTCAGCGCAGaacaa                                                                               >  1:853148/1‑149 (MQ=255)                                                                    gCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGATGGGCGGGCCCCAGCGCGGGGGGGCATCAGTAAAGTACCAACAGAAGC‑ACGGG‑TTTAGCGCAGaacaa                                                                               >  1:993328/1‑149 (MQ=255)                                                                      gATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGATGGGCGGGCCCCAGCGCGGGTGGGCCTCAGTAAGTTACCAACAAATGG‑ACGGG‑TTTAGCGCAGaaaaa                                                                               >  1:117114/1‑147 (MQ=255)                                                                          aaaCCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGC‑ACGGT‑TTCAGCGCAGAACAACTGCGTaa                                                                       <  2:1105537/151‑1 (MQ=255)                                                                             ccATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGATGGGCGGGGCCCCGCGCGTGTGGGCGACAGCAAGGTACCAACAGATGC‑ACGGG‑TTCTACGCAGAACAAATGCGTGATa                                                                     >  2:456521/1‑150 (MQ=255)                                                                              cATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGGTGGGCGGGCCTCCGCGCGCGGGGGCATCCGTAAGGTACCAACAGATGC‑ACGGG‑TTCAGCGCAGAACAACTGCGTAATAc                                                                    >  2:2473/1‑150 (MQ=255)                                                                                tGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGGTAGGCGGGCCTCCGCGCGCGTGTGCATTCGTAAGGTAGCTACAAAAGC‑CCGGG‑TTGCGCGCGGAAAAACTGCGGAATAc                                                                    >  1:815780/1‑148 (MQ=255)                                                                                tGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GAGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGC‑ACGGT‑TTCAGCGCAGAACAACTGCGTAATa                                                                     <  2:314405/147‑1 (MQ=255)                                                                                     gACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGC‑ACGGT‑TTCAGCGCAGAACAACTGCGTAATACCACTCATc                                                            <  1:151579/151‑1 (MQ=255)                                                                                      aCATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTGTGGCCGGGCCTCCGCGCGGGGGTGGCTCCGTAAATTAACCACAGCAGC‑ACGGG‑TTTCGCGCAGAAAAACAGCGTaaaaacccccc                                                              >  2:855306/1‑139 (MQ=255)                                                                                         tGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGGTGGGCCGGCCTCAGGGCGGGTGTGCATAAGGAAGGTACCCACAGAAGC‑ACGGG‑TTCAGCGCCGAACAACTGCGTAATAc                                                                    >  2:150144/1‑139 (MQ=255)                                                                                          gCTTGGATGGCTGGCGCAGGAAGAAGGGATTCGTCTTGAACCTTCGGCACTGGC‑GGGGGGTTGCGGGGCACCGCGAGGGTGTGGCTGTGCAAATTGCCCACACATGC‑ACGGG‑TTCACCGCAGAACAACACCGGAATACAAcaatcatagtg                                                       >  2:1100100/1‑141 (MQ=255)                                                                                             tggctggctggcGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGATTGGCGGGCCCCCGCGCGGGTGGGGATCAGTAAGTTAACAACAAATGC‑ACGGG‑TTCCGCGCCGAAAAAACGCGGAAAAACACCCATCTGGGGTgg                                                    >  1:39404/1‑151 (MQ=255)                                                                                              ggctggctggcGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGC‑ACGGT‑TTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTggg                                                   <  1:370008/151‑1 (MQ=255)                                                                                              ggctggctggcGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGC‑ACGGT‑TTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTgg                                                    <  1:1126854/150‑1 (MQ=255)                                                                                               gctggctggcGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGCTGGGCGGGCCTCCCAGGGGGTGTGTATCCATAAGTTACCCACAGATTC‑ACGGG‑TTCCGCGCAGAACAACCCCGCAATAACAACCATCTGGTgggg                                                    >  2:553102/1‑149 (MQ=255)                                                                                                  ggctggcGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTGTGGGCGGGCCTCAGCCCGGGTGGGGATCAGACAGTTACCACAACAAGG‑CCGGTGTTCAG‑GCAGAAAAACTGAGGCAAACCACCCAACTGGTGGGGGCg                                                 >  2:609654/1‑149 (MQ=255)                                                                                                  ggctggcGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGC‑ACGGT‑TTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCg                                                 <  1:567844/149‑1 (MQ=255)                                                                                                  ggctggcGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGC‑ACGGT‑TTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCg                                                 <  1:882712/149‑1 (MQ=255)                                                                                                  ggctggcGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGC‑ACGGT‑TTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGa                                                <  1:407685/150‑1 (MQ=255)                                                                                                     tggcGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGATGGGCGGACCCCAGCGCGGGTGTGGATCAGTAAGTTAACAACAGATGC‑ACGGG‑TTTAGCGCAGAAAAAATGCGCAATAACACTCATCTGGGGTGGGCGACggg                                            >  1:893248/1‑151 (MQ=255)                                                                                                      ggcGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGC‑ACGGT‑TTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGa                                           <  1:42828/151‑1 (MQ=255)                                                                                                       gcgcAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGATTGGCGGGCCTCCGCGCGGGTGGGTGTTACGAAGTTAACCACAGATGG‑ACGGG‑TTCAGCGCAGAACAACTTCGGAATACCACCCATCTTGGGTGGGCGGCggggg                                          >  2:365847/1‑149 (MQ=255)                                                                                                          cAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGCGGGGGGTGGCC‑GACCTCAATCGGGGTGTGCAGCAGGAAGATGTCACCAGCAGC‑CCGGG‑GTCACACGGGAACAACTGCCGAATACCACTCCACTGGTGGGGGGGGCGGGAgggg                                       >  2:553288/1‑149 (MQ=255)                                                                                                              aagaagGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTATGGCCGGGCCTCCGCGCGTGTGTGCATCAGTAAGATAACAACAAATGC‑ACGGG‑TTCAGCGCCGAACAAATGCGGAATAACACTCATCTTGGGGGGGGGGCGggggggggga                                    >  1:443945/1‑148 (MQ=255)                                                                                                                 aagGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGATGGGCGGGCCTCAGCGCGGGTGTGGATCAGTAAGTTAACAACAGATGC‑ACGGG‑TTTAGCGCCGAAAAACAGCGGAATACCACTCCTCTGGGGGGGGGGGCGGGGGGGGGGATgg                                 >  1:455100/1‑150 (MQ=255)                                                                                                                    gTATTCGTCTTGAACCTTCGGCACTGGC‑GGGGATGGGCGGGCCTCCGCGCGGGTGTGCCTCCGTAAGTTACCAACAGATGC‑ACGGG‑TTCAGCGCAGAACAACAGCGGAATACCACACATCTGGGGGGGGGGACGGGGGGGGGGATgg                                 >  2:460338/1‑147 (MQ=255)                                                                                                                       ttCGTCTTGAACCTTCGGCACTGGC‑GGGGATGGCCGGGCCTCCGCGCGTGTGGGCATCAGTAAGATACCAACAGACGG‑ACGGT‑TTTAGCGCAGAAAAAACGCGTAATACCACCCATCTTGTGTGGGGGAGGGGGGGGGGGATGGTGCCggg                          >  1:221089/1‑150 (MQ=255)                                                                                                                         cGTCTTGAACCTTCGGCACTGGC‑GGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGC‑ACGGT‑TTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGgaag                        <  1:926310/151‑1 (MQ=255)                                                                                                                               gAACCTTCGGCACTGGC‑GGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGC‑ACGGT‑TTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAgaga                  <  1:695275/151‑1 (MQ=255)                                                                                                                                      cGGCACTGGC‑GGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGC‑ACGGT‑TTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGATGAATc            <  2:73267/150‑1 (MQ=255)                                                                                                                                      cGGCACTGGC‑GGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGC‑ACGGT‑TTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGATGAATCa           <  2:720181/151‑1 (MQ=255)                                                                                                                                              gC‑GGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGC‑ACGGT‑TTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGATGAATCAATATCTgg   <  2:1165050/151‑1 (MQ=255)                                                                                                                                                 gggTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGC‑ACGGT‑TTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGATGAATCAATATCTGGc  >  1:85814/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GTTGGTCGCGCATCAGGCTTTGTCGGGCGGGCAATGGAGCGTCTGCTGGATGGCTTCTATACCCTTAGCGATCAAACCATGTATGACATGCTTGGCTGGCTGGCGCAGGAAGAAGGTATTCGTCTTGAACCTTCGGCACTGGC‑GGGTATGGCCGGACCTCAGCGCGTGTGTGCATCAGTAAGTTACCAACAGATGC‑ACGGT‑TTCAGCGCAGAACAACTGCGTAATACCACTCATCTGGTGTGGGCGACGGGAGGTGGAATGGTGCCGGAAGAAGAGATGAATCAATATCTGGC  >  NC_000913/2480225‑2480517

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