Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,133,302	0	T	G	55.7%	42.9 / 84.5	61	H219P (CAC→CCC)	yghS	putative ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (7/20);  new base G (34/0);  total (41/20)
Fisher's exact test for biased strand distribution p-value = 1.42e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-06

ATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAG‑TTCAAGGATCTTTGCGCCATTG‑AATGTCAGTTGC‑GGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGG  >  NC_000913/3133154‑3133450                                                                                                                                                     |                                                                                                                                                         aTGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGGGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgtgt                                                                                                                                                        <  2:594383/151‑1 (MQ=255)   gTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgggtgg                                                                                                                                                      >  2:445089/1‑149 (MQ=255)    tAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCg                                                                                                                                                      >  2:398207/1‑150 (MQ=255)    tAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCg                                                                                                                                                      >  1:308244/1‑150 (MQ=255)      gTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGgggttggc                                                                                                                                                    >  1:848457/1‑146 (MQ=255)       ttGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGccc                                                                                                                                                   <  2:994604/150‑1 (MQ=255)       ttGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGcc                                                                                                                                                    <  2:94193/149‑1 (MQ=255)         gggCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGTGCCCAt                                                                                                                                                 >  1:85245/1‑150 (MQ=255)           gCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATcc                                                                                                                                               <  1:776794/150‑1 (MQ=255)               cggcgCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAG‑TTc                                                                                                                                         <  2:935523/151‑1 (MQ=255)                  cgcgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAG‑TTCaa                                                                                                                                       <  2:22679/150‑1 (MQ=255)                    cgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAggttcaggg                                                                                                                                     >  1:693384/1‑143 (MQ=255)                      gTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGGCCCCCATCCAG‑TTTAAgggt                                                                                                                                   >  1:338865/1‑148 (MQ=255)                      gTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTGGCCCATCCAC‑TTCAAgggta                                                                                                                                  >  2:863428/1‑148 (MQ=255)                      gTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCCCCAG‑TTCaatggtt                                                                                                                                  >  2:731604/1‑146 (MQ=255)                           aCTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCCG‑TTCCAGGATCTTTgg                                                                                                                             >  1:501372/1‑150 (MQ=255)                              aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTGGCCCCTCCAA‑TTTAAGGATATTTgtgc                                                                                                                           >  1:391031/1‑150 (MQ=255)                              aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAG‑GTCAAGGATCTTTgcgc                                                                                                                           >  2:459793/1‑150 (MQ=255)                                   cAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAG‑TTTAAGGATCTTTTCGCCCTTG‑a                                                                                                                    >  1:399164/1‑151 (MQ=255)                                     gtttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCCCCCCG‑GTCAAGGGTTTTTTTGCCAtttgat                                                                                                                   >  2:4453/1‑147 (MQ=255)                                         gagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCATCCAG‑GTTAAGGGTGTTTTTGCCATTT‑aaaggc                                                                                                               >  1:1044172/1‑146 (MQ=255)                                         gagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCCCCCG‑TTCAAAGGTCTTTGTGCCCTTT‑AATGtta                                                                                                              >  1:828221/1‑149 (MQ=255)                                         gagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCAAG‑GTCAAGGGTCTTTTCGCCATTG‑GATGGCa                                                                                                              >  2:89355/1‑151 (MQ=255)                                          agtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCAG‑TTCAAGGGTCTTTTCGCCCTTT‑AATGTCAg                                                                                                             >  2:61074/1‑151 (MQ=255)                                            tttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAG‑TTCAAGGATCTTTGCGCCATTG‑AATGTCAg                                                                                                             <  2:596976/149‑1 (MQ=255)                                              tgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAG‑TTCAAGGATCTTTGCGCCATTG‑AATGTCAGTTgc                                                                                                         <  2:487357/151‑1 (MQ=255)                                               gTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAC‑TTCAAAGGACCTTGCGCCATTG‑AATGTCAGGTgc‑g                                                                                                       >  1:663324/1‑151 (MQ=255)                                               gTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCT‑TTTCAAGATGTTTTTGGCATTT‑AATGGCAGTTgc‑g                                                                                                       >  1:230690/1‑151 (MQ=255)                                               gTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCG‑TTCAAGGAGTTTTTCGCCCTTG‑AATGTTCGTTgc‑g                                                                                                       >  1:651571/1‑151 (MQ=255)                                               gTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCCCCAG‑TTCAAGGGTCTTTTTGCCCTTT‑AATGTCAGTTgc‑g                                                                                                       >  2:630380/1‑151 (MQ=255)                                                  gTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGGCGGCTCTCCAGCTTTAAGGGCTTTTTGTCCATT‑‑AATGTCCGTTgc‑ggcg                                                                                                    >  1:1047587/1‑151 (MQ=255)                                                    gCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAG‑TTCAAGGAACTTTTCGCCATTT‑AATGTCAGTTGC‑CGCGTa                                                                                                  >  1:618134/1‑151 (MQ=255)                                                        ccATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAG‑TTCAAGGATCTTTGCGCCATTG‑AATGTCAGTTGC‑GGCGTAACgg                                                                                              <  2:200141/151‑1 (MQ=255)                                                           ttCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAG‑TTCCAGGATCTTTGCGCCATTG‑GATGTCAGTTGT‑GGCGTTAAGGCgg                                                                                           >  1:298783/1‑150 (MQ=255)                                                                     catcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAG‑TTCAAGGATCTTTGCGCCATTG‑AATGTCAGTTGC‑GGCGTAACGGCGATTTTTTCCTg                                                                                 <  1:261479/151‑1 (MQ=255)                                                                        ccccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCG‑GTCAAGGAAATTTGCGCCCATG‑AATGTGAGGTGC‑GGGGGAACGGCGGTTTTTTTTTGTaa                                                                              >  2:1023078/1‑151 (MQ=255)                                                                        ccccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCA‑GTCCAGGATCTTTTCGCCCTTG‑AATGTCAGGTGC‑GGCGGAACCGCGAATTTTTTCTGTaa                                                                              >  1:723456/1‑151 (MQ=255)                                                                          ccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCG‑TTCAAGGATCTTTGCGCCATTG‑GATGTCAGTTGCGGGGGAAAGGGGGATTTTTTCTTGTaag                                                                             >  2:582533/1‑150 (MQ=255)                                                                          ccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCAG‑TTCAAGGATCTTTTCGCCCTTT‑AATGTCAGTTGC‑GGCGTGACGGCGGTTTTTTTCTGTAACg                                                                            >  1:545761/1‑151 (MQ=255)                                                                               ccATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGGGCCCCTCCCG‑TTCTAGGGTCTTTGCGCCATTG‑AATGTTAGTTGG‑GGCGCAACGGCGGTTTTTTTCTGttaagcggc                                                                        >  2:73583/1‑142 (MQ=255)                                                                                 aTTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAG‑TTCAAGGATCTTTGCGCCATTG‑AATGTCAGTTGC‑GGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAg                                                                    <  1:444922/151‑1 (MQ=255)                                                                                     tGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAG‑TTCAAGGATCTTTGCGCCATTG‑AATGTCAGTTGC‑GGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTtgg                                                                <  1:854837/151‑1 (MQ=255)                                                                                         gATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCCCCAG‑TTCAAGGGATTTTTCGCCCTTG‑AATGTCCGTTGG‑GGGGGAACGGCGGTTTTTTTCTGTTAAGCCG‑CCCGtttggggg                                                            >  2:68763/1‑145 (MQ=255)                                                                                                    cgcgGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAG‑TTCAAGGATCTTTGCGCCATTG‑AATGTCAGTTGC‑GGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTAcg                                                  <  1:200171/150‑1 (MQ=255)                                                                                                      cggcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCATCCCG‑TTCAAGGGTCTTTTTGCCCTTG‑AATTTCAGTTGG‑GGGGGAACGGCGGTTTTTTTCTGTAAAGCTG‑CCAGGtggtggggggggtttcgcg                                                >  2:242654/1‑137 (MQ=255)                                                                                                        gcgTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCATCCAG‑TTCCAGGATCTTTGCGCCATTG‑GATGTCAGTTGC‑GGGGTAAAGGCGATTTTTTTTTGTAAAGCTG‑GCCGTTGGTGGTTAGGGTTACGGGCaa                                             >  2:833234/1‑151 (MQ=255)                                                                                                            gAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCTCCCCTCCAG‑TTTAAGGATTTTTGGGGCCTTT‑TATGTCAGTTGC‑GGCGGAACGGCGAGTTTTTCCTGTTACGCTG‑CCCGTTGGGGGTCACGGTTTCGCGCAAAcg                                          >  2:736670/1‑150 (MQ=255)                                                                                                              aTTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTTTTGCCCCACCCG‑TTCAAGGGTCTTTTGGCCAATT‑AAAGGCTGTTTT‑GGCGTTACGGCGAATTTTTTCTGTAACGCTG‑CCAGCTGGTGGTCAGGTTTTAGCGCCAAAGCgg                                       >  2:922326/1‑151 (MQ=255)                                                                                                                    cgTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAG‑TTCAAGGATCTTTGCGCCATTG‑AATGTCAGTTGC‑GGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGt                                  <  2:374600/150‑1 (MQ=255)                                                                                                                     gTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCCG‑TTCAAGGGTCTTTTCGCCATTT‑AATGTCAGTTGC‑GGGGTAACGGCGAGTTTTTTCTGTCAAGCTGTCCGGTTGG‑GGGCCGGTTTTCGCGCCAACGCGGgttgttt                                >  1:119426/1‑150 (MQ=255)                                                                                                                     gTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCATCCAG‑TTCAAGGATCTTTGGGCCCTTGTATTTTCGGTTG‑‑GGGGGAACGGCGGTTTTTTTCTGTAAAGCTG‑TCCGCTTGTGGGTAGGGTTTCCGGCCAACGCGGGTTggtt                                >  1:934245/1‑148 (MQ=255)                                                                                                                      tAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCAACCCG‑GGCAAGGAGCTTTGCGGCATTT‑AATGTGAAGTGC‑GGGGCAACGGGGAGTTTTTTCTGTAAAGCTG‑CCAGTTGGGGGGTGGGTGTACGGGCAAACTCGGGTttttt                                >  1:1001508/1‑149 (MQ=255)                                                                                                                         tGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAG‑TTCAAGGATCTTTGCGCCATTG‑AATGTCAGTTGC‑GGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTtcatc                             <  1:281410/150‑1 (MQ=255)                                                                                                                          taCGCTTGCAGAATTTAATCGGCGGGGTGTCGCCCATCCAG‑TTCAAGGATCTTTGCGCCATTG‑AATGTCAGTTGC‑GGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCaa                           <  1:564970/150‑1 (MQ=255)                                                                                                                           aCGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAG‑TTCAAGGATCTTTGCGCCATTG‑AATGTCAGTTGC‑GGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCaa                           <  1:249871/150‑1 (MQ=255)                                                                                                                           aCGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAG‑TTCAAGGATCTTTGCGCCATTG‑AATGTCAGTTGC‑GGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCaa                           <  2:306786/150‑1 (MQ=255)                                                                                                                            cGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAG‑TTCAAGGATCTTTTCGCCATTG‑TATTTCCGTTGC‑GGCGTAACGGCGATTTTTTTCTGTAAAGCTG‑CCCGTTGGGGGTCAGGTTTTCGCGCAAACGCGGGTTTTTTATCaaat                         >  2:13533/1‑149 (MQ=255)                                                                                                                                    gAATTTCATCGGCGGGGGGGCGCCCATCCAG‑TTCAAGGAACTTTGGGCCATTG‑AATGTCAGTTGC‑GGCGGAACGGGGGGTTTTTTCTGTTACGCTG‑CCAGTTGGTGGTCAGGGTTACGCGCGAACGCGGTTTGTTCATCCATGCCAAGaag                 >  1:634830/1‑149 (MQ=255)                                                                                                                                     aaTTTCATCGGCGGGGTGTCGCCCATCCAG‑TTCAAGGATCTTTGCGCCATTG‑AATGTCAGTTGC‑GGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGAc                  <  1:902494/149‑1 (MQ=255)                                                                                                                                           aTCGGCGGGGTGTCGCCCATCCAG‑TTCAAGGATCTTTGCGCCATTG‑AATGTCAGTTGC‑GGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATcaac          <  2:385377/151‑1 (MQ=255)                                                                                                                                                cggGGTGTCGCCCATCCAG‑TTCAAGGATCTTTGCGCCATTG‑AATGTCAGTTGC‑GGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATAc     <  2:1044552/151‑1 (MQ=255)                                                                                                                                                  gggTGTCGCCCATCCAG‑TTCAAGGATCTTTGCGCCATTG‑AATGTCAGTTGC‑GGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACgg   >  2:400545/1‑151 (MQ=255)                                                                                                                                                   ggTGTCGCCCATCCAG‑TTCAAGGATCTTTGCGCCATTG‑AATGTCAGTTGC‑GGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACggg  >  2:298757/1‑151 (MQ=255)                                                                                                                                                   ggTGTCGCCCATCCAG‑TTCAAGGATCTTTGCGCCATTG‑AATGTCAGTTGC‑GGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACggg  >  2:171326/1‑151 (MQ=255)                                                                                                                                                    gtgtCGCCCATCCAG‑TTCAAGGATCTTTGCGCCATTG‑AATGTCAGTTGC‑GGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACggg  >  1:369726/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                         ATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAG‑TTCAAGGATCTTTGCGCCATTG‑AATGTCAGTTGC‑GGCGTAACGGCGATTTTTTCCTGTAACGCTG‑CCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATCAATGCCAAGACGAATCAACAATACGGG  >  NC_000913/3133154‑3133450

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