Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	838,239	0	G	T	54.4%	6.7 / 70.4	46	A73D (GCT→GAT)	ybiI	DksA‑type zinc finger protein
Reads supporting (aligned to +/- strand):  ref base G (8/13);  new base T (11/14);  total (19/27)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAAGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAACCGGATGCCATATCGACCTCCCCATATCAATACTTGTACAGTTAAGTGTAGCTAATCCAGGGACGAACTCGGGCAGTTCAAGCATCAGATCTC  >  NC_000913/837953‑838538                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                            ataAGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTaaa                                                                                                                                                                                                                                                                                               <  2:323823/301‑1 (MQ=255)    aGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAAt                                                                                                                                                                                                                                                                                              >  2:173929/1‑299 (MQ=255)      tGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCtt                                                                                                                                                                                                                                                                                           <  2:42272/300‑1 (MQ=255)        aCGTTGTTTTGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCTTtc                                                                                                                                                                                                                                                                                         <  1:280066/300‑1 (MQ=255)                      aTGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATg                                                                                                                                                                                                                                                                           <  1:165512/300‑1 (MQ=255)                           atgaGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATAc                                                                                                                                                                                                                                                                      <  2:41591/300‑1 (MQ=255)                            tgaGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACAt                                                                                                                                                                                                                                                                    <  2:104978/301‑1 (MQ=255)                                tgctcttccgatctACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAgc                                                                                                                                                                                                                                                                <  2:106985/288‑1 (MQ=255)                                 tgctcttccgatctcGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAgcg                                                                                                                                                                                                                                                               <  2:325556/287‑1 (MQ=255)                                         tttttACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAg                                                                                                                                                                                                                                                       <  1:173763/301‑1 (MQ=255)                                          ttttACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAgg                                                                                                                                                                                                                                                      <  2:236319/301‑1 (MQ=255)                                              aCGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAgcagatcggaagagcg                                                                                                                                                                                                                                                  >  1:106878/1‑287 (MQ=255)                                              aCGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATACGCGCACGCCAGGAATg                                                                                                                                                                                                                                                  >  1:40162/1‑301 (MQ=255)                                               cGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAgcgagatcggaagagcg                                                                                                                                                                                                                                                 >  1:325192/1‑287 (MQ=255)                                                  aagaTAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCtt                                                                                                                                                                                                                                              <  2:74622/300‑1 (MQ=255)                                                           tCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACggg                                                                                                                                                                                                                                     <  2:31691/301‑1 (MQ=255)                                                                      tGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATg                                                                                                                                                                                                                          >  2:348227/1‑301 (MQ=255)                                                                         ccGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCAAGGAATGGCTTCCCGACGGGCCTGCGGGATgggg                                                                                                                                                                                                                       >  1:307357/1‑301 (MQ=255)                                                                         ccGCCTTGTTCAAATAAATTGATGATTATCTAACCTACTGATCAGTCATGTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGTCAGTCACAGCAGCGTTAAAGGCACTGAATTACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCATTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAGTCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATgggg                                                                                                                                                                                                                       <  1:53490/301‑1 (MQ=255)                                                                                                    atCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAgg                                                                                                                                                                                            >  1:270914/1‑301 (MQ=255)                                                                                                                        atttatttGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTc                                                                                                                                                                        <  1:51259/301‑1 (MQ=255)                                                                                                                           tatttGTGGTACTTATGAACTTGTTGTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCAcc                                                                                                                                                                     <  1:355927/301‑1 (MQ=255)                                                                                                                                gttgtACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCCCCGCCCTCTTCACATTCCTCCAGGCTTTCGCCGCGCGGAATTTCACCCCGa                                                                                                                                                                 >  1:186798/1‑300 (MQ=255)                                                                                                                                        tATGAAATTGTTTTATACCGCGTGGCAGTCACAGCATCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCa                                                                                                                                                        <  2:266024/301‑1 (MQ=255)                                                                                                                                           gAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAAt                                                                                                                                                      >  1:148223/1‑300 (MQ=255)                                                                                                                                                  ttttATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTAACCCCGAGCGCGGGCAATCGCATCtt                                                                                                                                              >  1:415059/1‑301 (MQ=255)                                                                                                                                                       taCCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAAtt                                                                                                                                         <  2:40201/301‑1 (MQ=255)                                                                                                                                                        aCCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTg                                                                                                                                        <  2:307698/301‑1 (MQ=255)                                                                                                                                                                    tCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCt                                                                                                                            <  2:271203/301‑1 (MQ=255)                                                                                                                                                                              gTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGACTTTCACCCCGAGCCCGGGCATTCGCATCTTCAATTCTACTGTTGACTTGTTCGTTgg                                                                                                                   >  1:281901/1‑299 (MQ=255)                                                                                                                                                                                              aacacGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGcc                                                                                                  <  2:186976/300‑1 (MQ=255)                                                                                                                                                                                                            aaaTTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTTACCTGTTCGTTGACGGCCTCATCGTTAGCCCACCCGGCTGCCa                                                                                     >  1:366057/1‑300 (MQ=255)                                                                                                                                                                                                               tgtgctcttccgatcttctcGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAACCGGATGCCATATc                                                                                 <  2:384781/285‑1 (MQ=255)                                                                                                                                                                                                                             cATCTAGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAACCGGATGCCATATCGACCTCCCCATATc                                                                   <  2:185551/301‑1 (MQ=255)                                                                                                                                                                                                                               tctcGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAACCGGATGCCATATcagatcggaagagcgt                                                                  >  1:384354/1‑285 (MQ=255)                                                                                                                                                                                                                                ccgatatgTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGGATAAGCTGGTTTTTGTAAATCCTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAACCGGATGCCATATCGACCTCCCCATATCaa                                                                 <  1:32142/293‑1 (MQ=255)                                                                                                                                                                                                                                 tcGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAACCGGATGCCATATCGACCTCCCCATATCAATa                                                               <  2:148364/301‑1 (MQ=255)                                                                                                                                                                                                                                  cGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAACCGGATGCCATATCGACCTCCCCATATCAATa                                                               >  2:403224/1‑300 (MQ=255)                                                                                                                                                                                                                                     cAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTCCGTTGACGGCGTCATCGTTAGCCCAACCGGATGCCATATCGACCTCCCCATACCAATACttt                                                           >  1:251535/1‑300 (MQ=255)                                                                                                                                                                                                                                       gTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAACCGGATGCCATATCGACCTCCCCATATCaaagatcgga                                                         >  2:32171/1‑293 (MQ=255)                                                                                                                                                                                                                                        ttcagacgtgtgctcttccgatcttttCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAACCGGATGCCATATCGACCTCCCCATATCAATACTTGt                                                          <  2:232032/275‑1 (MQ=255)                                                                                                                                                                                                                                                          gCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAATCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAACCGGATGCCATATCGACCTCCCCATATCAATACTTGTACAGTTAAGTGTAGCTAAt                                       <  2:251796/300‑1 (MQ=255)                                                                                                                                                                                                                                                               cTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAACCGGATGCCATATCGACCTCCCCATATCAATACTTGTACAGTTAAGTGTAGCTAATCCAggg                                 >  2:93584/1‑301 (MQ=255)                                                                                                                                                                                                                                                                tttCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCACAACCGGATGCCATATCGACCTCCCCATCTCACTACTTGTagatcggaagagcgtcgtgtcgtgaa                                >  1:231810/1‑276 (MQ=255)                                                                                                                                                                                                                                                                                    atCCTGTATAATCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGCCAATCGCATCTTCAATTGTACTGTTGACCTGTTCGTTGACGGCGTCATCGTTATCCCAACCGGCTGCCATATCGACCTCCCCATATCAATACTTTTACAGTTAAGTGTAGCTAATCCATGGACGAACTCCGGCATTTCAAg             >  1:337345/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                            taAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAACCGGATGCCATATCGACCTCCCCATATCAATACTTGTACAGTTAAGTGTAGCTAATCCAGGGACGAACTCGGGCAGTTCAAGCATCAGAtc    <  1:93485/301‑1 (MQ=255)                                                                                                                                                                                                                                                                                               gCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAACCGGATGCCATATCGACCTCCCCATATCAATACTTGTACAGTTAAGTGTAGCTAATCCAGGGACGAACTCGGGCAGTTCACGCATCAGAtctc  >  1:174048/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                                            ATAAGTGACGTTGTTTCGTCGATGTGATGAGTATCACGTTTTTTTACGAGAGATAAAATCGAAGAGAATTGACCGCCTTGTTCAAATAAATTGATTGATATCTAACCTACTGATAAACCATTTATTTGTTGTACTTATGAACTTGTTTTATACCGCGTGGCAGTCACAGCAGCGTTAAAGGCACGGAATGACACGCAAAATACAAATTATACTCGCTATGCATCTCGCCAGTTGTAGTTAACGTAACTGGCTGTCTTTCGAACCTCTGCGATTATATCCTGTATAAGCTGGTTTTTGTAAATCTTTCTCCTGCTGACAATGAATACATAAGCGCACGCCAGGAATGGCTTCCCGACGGGCCTGCGGGATGGGGGCACCGCACTCTTCACATTCATCCAGGCTTTCGCCGCGCGGAATTTCACCCCGAGCGCGGGCAATCGCATCTTCAATTGTACTGTTGATCTGTTCGTTGACGGCGTCATCGTTAGCCCAACCGGATGCCATATCGACCTCCCCATATCAATACTTGTACAGTTAAGTGTAGCTAATCCAGGGACGAACTCGGGCAGTTCAAGCATCAGATCTC  >  NC_000913/837953‑838538

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑