Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,812,792	0	C	A	38.2%	119.9 / 28.0	68	T21N (ACC→AAC)	dfp	fused 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase, FMN‑binding
Reads supporting (aligned to +/- strand):  ref base C (10/32);  new base A (26/0);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 8.60e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

CGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCAGCCGAAGC  >  NC_000913/3812642‑3812942                                                                                                                                                         |                                                                                                                                                           cGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAAc                                                                                                                                                            <  1:1017267/151‑1 (MQ=255)    aaaGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAAcccc                                                                                                                                                         <  1:807070/151‑1 (MQ=255)       gTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGa                                                                                                                                                      <  2:82618/151‑1 (MQ=255)          aaTCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTTCCTTTTAAAAccccgtacc                                                                                                                                                    >  2:978982/1‑146 (MQ=255)             ctcAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGA‑ATTGCTGCC‑CATAAAAACCCCTAACttg                                                                                                                                                 >  2:625240/1‑148 (MQ=255)                 ggTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCCCCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGcgt                                                                                                                                            <  1:42850/151‑1 (MQ=255)                 ggTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGcgt                                                                                                                                            <  1:391958/151‑1 (MQ=255)                   tGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑AATGGTGCC‑TTTAAAAACCCCGAACTGGGgcggcg                                                                                                                                          >  1:771899/1‑151 (MQ=255)                     aaGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTTCTGCC‑TATAAAACCCCCGAAATGGGGCGTCGtt                                                                                                                                        >  1:218038/1‑151 (MQ=255)                       ggCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTg                                                                                                                                      <  2:92024/151‑1 (MQ=255)                         cACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTgc                                                                                                                                     <  1:291826/150‑1 (MQ=255)                                acaaAAAGAATGTGAAAAAATTGCCAACTTCTGGTTGCCCCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATc‑g‑c                                                                                                                           <  1:962786/151‑1 (MQ=255)                                    aaGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAAAACCTGCACTCGGGCGTCGTTTGCGCGATC‑G‑Cgggg                                                                                                                       >  2:433031/1‑151 (MQ=255)                                     aGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCCGAACTTGGGCGGCGCTTGCGCGATC‑G‑Cgggg                                                                                                                       >  1:270093/1‑150 (MQ=255)                                     aGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAAACCCCGAAATGGTGGGTCGTTTTCGCGGTC‑G‑Cgggg                                                                                                                       >  1:103934/1‑150 (MQ=255)                                      gTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGcc                                                                                                                     <  2:806318/151‑1 (MQ=255)                                       tATTGTGTTAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAAACCCTGAACTGGTGCGTCGTTTTCGCGATC‑G‑CGGGGcc                                                                                                                     >  1:447201/1‑150 (MQ=255)                                       tATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAAACCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑Ccggggc                                                                                                                     >  1:472032/1‑148 (MQ=255)                                        aTTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGa                                                                                                                   <  2:1033866/151‑1 (MQ=255)                                             gATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TTTAAAAACCCCGAACTGGGGCGGCGTTTGCGCGATC‑G‑CGGGGCCCACGGcc                                                                                                              >  1:196776/1‑151 (MQ=255)                                             gATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAAACCCTGAACTGGTGCGTCGTTTTCGCGATC‑G‑CGGGGCCGACGTcc                                                                                                              >  1:258858/1‑151 (MQ=255)                                              aTAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAAACCCCGAACTGGTGGGGCGGTTGCGCGATC‑G‑CGGGGGCGCCggc                                                                                                               >  1:153898/1‑147 (MQ=255)                                               tAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTTCTGCC‑TCTAAAAACCCCCAACCTGGGTGTCGGTTTCGGGATT‑G‑CGGGGgcggcgcc                                                                                                               >  2:26139/1‑146 (MQ=255)                                                aaaaTCGCCAACTTCTGGTGTCTCCCTACAGGAAAAATCATCATGGGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAGCCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGCCCGACGTCcgcg                                                                                                           <  2:1088827/151‑1 (MQ=255)                                                aaaaTCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCcgc                                                                                                            <  1:65360/150‑1 (MQ=255)                                                   aTCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTTCTGCC‑TATAAAAACCCCGAAATGGGGCGTCGTTTGCGCGAGC‑GCCGGGGccccgcccg                                                                                                              >  1:1202264/1‑139 (MQ=255)                                                     cGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGc                                                                                                       <  1:274065/150‑1 (MQ=255)                                                       ccAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAACCCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCa                                                                                                     <  1:1152563/150‑1 (MQ=255)                                                       ccAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCAt                                                                                                    <  1:617115/151‑1 (MQ=255)                                                        cAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TCTAAAACCCCTGCACTGGTGCGTCGCTTTGGCGATC‑G‑CGGGGGCCACGTCCGCGTAGCCATg                                                                                                   >  1:193070/1‑151 (MQ=255)                                                          aCTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTGATTG‑TGCC‑CATAAAAACCCTGCAGTGGGGTTTTGTTTTCTCGATT‑G‑CGGGGc                                                                                                                      >  2:437270/1‑130 (MQ=255)                                                          aCTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGGC‑CATAAAAACCCCGAAATGGGGCGCCGTTTTTGCGATC‑G‑CGGGGGCGCCGTCCGCGTAGCCATGa                                                                                                  >  1:1255617/1‑150 (MQ=255)                                                          aCTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAAACCCCGAAATGGGGCGGCGTTTTCGCGATC‑G‑CGGGGCCGGCGTCCGCGGAGcc                                                                                                      >  1:338899/1‑146 (MQ=255)                                                                ggTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAACCCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGaa                                                                                             <  2:223202/149‑1 (MQ=255)                                                                ggTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGaa                                                                                             <  1:707842/149‑1 (MQ=255)                                                                ggTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAg                                                                                            <  1:124098/150‑1 (MQ=255)                                                                ggTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑CATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGaa                                                                                             <  1:503544/149‑1 (MQ=255)                                                                    tCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTTCTGCC‑CATAAAAACCCTGGCCCTGGGCGGCGGTTTCGCGATCGG‑CGGGGCCACCGGCGGC‑TAGCCAAGACCGAAggcgg                                                                                        >  2:832989/1‑146 (MQ=255)                                                                    tCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGT‑AAAAAATCGTTCTCGGCGTTAGCGGCGGG‑TTTTCTTGC‑TATTAAACCCCTTAACCTGTGCGGCGGTTTCGCGAAT‑G‑CGGGGGCCGCGTCCGCGTAGCCCTGCACCAAGC‑ggg                                                                                       >  1:156255/1‑148 (MQ=255)                                                                     cACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTTCTGGC‑TATAAAAACCCCGTACTTGGGCGGCGGTTTCGCGGTC‑G‑CGGGGGCGACGTCCGCGGAGCCATGAACGCAAC‑gg                                                                                        >  1:231873/1‑148 (MQ=255)                                                                     cACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAAACCCCGAAATGGGGCGTCGTTTTCGCGATC‑G‑CGGGGGCGACGGCCCCGGAGCCCTGACCGAAGC‑gg                                                                                        >  2:646220/1‑148 (MQ=255)                                                                       cccTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCaaaa                                                                                   <  2:752730/151‑1 (MQ=255)                                                                             aGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCAAAAGCCttt                                                                             <  2:894929/151‑1 (MQ=255)                                                                                aaaaaTCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTTATGGC‑CATTAAACCCCCCTACTGGGGCGTCGTTTTCGCGAGC‑G‑CGGGGCCGCCGGCCCCGGCGGCATGACCGAAGC‑GGCGAAAGCCctt                                                                             >  1:559306/1‑148 (MQ=255)                                                                                  aaaTCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTTCTGCC‑CATAAAAACCCCTAACTGGTGCGTCGTTTTCGCGGGC‑G‑CGGGGGCGACGTCCGCGCAGCCATGACCGAAGC‑GGCAAAAAACttttt                                                                           >  2:458050/1‑146 (MQ=255)                                                                                  aaaTCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TCTAAAAACCCCGAAACTGTGGGGTGGTTGTGCGCTT‑G‑CGGGGGCGGGGGTCGCCGAGCCATGACCGGAGCGGGCAAAAGCttttt                                                                            >  1:1205144/1‑147 (MQ=255)                                                                                   aaTCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCAAAAGCCTTTATCAcc                                                                       <  1:524783/151‑1 (MQ=255)                                                                                   aaTCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCAAAAGCCTTTATCAcc                                                                       <  1:393926/151‑1 (MQ=255)                                                                                     tcatcaTGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑CTTAAAAACCCCGAACTGGGGCGGCGTTTGGGCGGTC‑G‑CGGGGGCGGCGTCCGCCGAGCCCTGACCGCAGC‑GGCCAAAACCCTTTTCAcccc                                                                     >  2:444044/1‑151 (MQ=255)                                                                                      catcatGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAAACCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCCCGTAGCCATGACCGAAGC‑GGCAAAAGCCTTTATCAccccc                                                                    >  1:447150/1‑150 (MQ=255)                                                                                              gCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCAAAAGCCTTTATCACCCCACTTAGCt                                                             <  2:1156687/150‑1 (MQ=255)                                                                                                 tGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAAAACCCCAACTGGTGCGTCGTTTGCGCGAGC‑G‑CGGGGGCGACGTCCGCGCAGCCATGACCGCAGC‑GGGCAAAACCTTTTTTACCCCCCCTAGCtt                                                            >  2:732077/1‑148 (MQ=255)                                                                                                          aaaaaaTCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTc                                               <  1:20874/151‑1 (MQ=255)                                                                                                           aaaaaTCGTTCTCGGCGTTAGCGGCGGT‑ATTGGTGCC‑TATAAAAACCCCCCAATGGTGCGGCGTTTGTGCGATC‑G‑CGGGGGCGACGTCCCCGTAGCCATGACCGAAGA‑GGGAAAAGACATTATTACCCCACCTAAGTTGCAGGC‑GGgtttt                                              >  2:775484/1‑149 (MQ=255)                                                                                                           aaaaaTCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAAAACCCGAACTGGGGGGGCGTTTGCGCGATC‑G‑CGGGGGCGACGTCCGCGTAGCCATGGCCGAAGC‑GGGAAAAAACTTTTTTACCCCACTTAGCTTGCACGC‑GGgtttt                                              >  1:717598/1‑149 (MQ=255)                                                                                                            aaaaTCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TCTAAAAACCCCGAACTGGGGCGGCGTTTGCGCGGAC‑G‑CGGGGGCGACCGCCGCCCAGGCCTGAACGAAGC‑GGGAAAAGACCTTATCACCCCCCCTAGCTTTCAGGC‑gggt                                                 >  2:653700/1‑145 (MQ=255)                                                                                                              aaTCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGt                                           <  1:1200949/151‑1 (MQ=255)                                                                                                                 cGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTa                                         <  2:559290/150‑1 (MQ=255)                                                                                                                 cGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTAt                                        <  1:657860/151‑1 (MQ=255)                                                                                                                       cGGCGTTAGCGGGGGT‑ATTGCTGCC‑TATAAAGCCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCGAAAGCCTTTATCACCCCTCCTAGCTTGCAGGC‑GGTTTCTGGTTATccc                                     <  1:262965/148‑1 (MQ=255)                                                                                                                       cGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATccc                                     <  1:583219/148‑1 (MQ=255)                                                                                                                       cGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCg                                    <  2:95119/149‑1 (MQ=255)                                                                                                                          cGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTcc                               <  2:1165057/151‑1 (MQ=255)                                                                                                                          cGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTcc                               <  1:433045/151‑1 (MQ=255)                                                                                                                           gTTAGCGGCGGT‑ATTTCTGCC‑TATAAAAACCCCGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCCTGACCGAAGC‑GGCAAAAAACTTTTTCACCCCACTTAGCTTTCAGGCGGGTTTTTGGTTTTCCCCtt                                  >  1:1194296/1‑148 (MQ=255)                                                                                                                                cggcggT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGGGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGGTTTCCCCCTTCCGCCAGt                         >  2:1169939/1‑151 (MQ=255)                                                                                                                                 ggcggT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGGTTCCCAc                            >  2:1220118/1‑147 (MQ=255)                                                                                                                                              gcC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCg             <  2:36567/150‑1 (MQ=255)                                                                                                                                                    taAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTTTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCAGcc       >  1:190423/1‑151 (MQ=255)                                                                                                                                                    taAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCAGcc       >  1:979358/1‑151 (MQ=255)                                                                                                                                                     aaaaCCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCCTTTTCCACAGGCTGCTGGACCCGGGAGcc       >  2:924674/1‑150 (MQ=255)                                                                                                                                                         ccccTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCAGCCGAAGc  <  1:176737/151‑1 (MQ=255)                                                                                                                                                         ccccTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTCCCACAGTCTGCTGGACCCGGCAGCCGAAGc  >  1:1110123/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                           CGCAAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGT‑ATTGCTGCC‑TATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATC‑G‑CGGGGCCGACGTCCGCGTAGCCATGACCGAAGC‑GGCAAAAGCCTTTATCACCCCACTTAGCTTGCAGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCAGCCGAAGC  >  NC_000913/3812642‑3812942

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