Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,812,792	0	C	A	45.7%	96.8 / 51.3	81	T21N (ACC→AAC)	dfp	fused 4'‑phosphopantothenoylcysteine decarboxylase/phosphopantothenoylcysteine synthetase, FMN‑binding
Reads supporting (aligned to +/- strand):  ref base C (15/29);  new base A (37/0);  total (52/29)
Fisher's exact test for biased strand distribution p-value = 2.90e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

AAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑AGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCAGCCGAAG  >  NC_000913/3812645‑3812941                                                                                                                                                    |                                                                                                                                                        aaaGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAcccc                                                                                                                                                      <  1:633504/151‑1 (MQ=255) aaaGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAcccc                                                                                                                                                      <  1:239165/151‑1 (MQ=255)     tGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCCGa                                                                                                                                                   >  2:43320/1‑150 (MQ=255)       aaTCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATTAAAACCCCGaa                                                                                                                                                  >  2:620380/1‑149 (MQ=255)              ggTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCCGAACTGGGGcg                                                                                                                                          >  1:124858/1‑150 (MQ=255)                tGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGGATTGCTGCCTATAAAAACCCCGAACTGGTGc                                                                                                                                           >  2:435130/1‑147 (MQ=255)                 gAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGcgtcgt                                                                                                                                      <  1:116453/151‑1 (MQ=255)                  aaGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGTTGCCCCTAAAAACCCCGCAATGGGGCGTCGtt                                                                                                                                     >  1:807237/1‑151 (MQ=255)                   aGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCCATAAAAACCCTGAAATGGTGCGTCGttt                                                                                                                                    >  2:536984/1‑151 (MQ=255)                       aCGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTgcgc                                                                                                                                <  1:784267/151‑1 (MQ=255)                        cGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTgcgc                                                                                                                                >  2:957530/1‑150 (MQ=255)                         gACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCCGAACTGGGGCGGCGTTTGcgcg                                                                                                                               >  2:195873/1‑150 (MQ=255)                           caacaaAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATc                                                                                                                            <  1:118809/151‑1 (MQ=255)                           caacaaAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCCGAACTGGTGCGGCGTTTGCGCGAt                                                                                                                             >  2:389551/1‑150 (MQ=255)                            aacaaAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCCGGAATGGTGGGTCGTTTGCGCGa                                                                                                                              >  2:614116/1‑148 (MQ=255)                              caaAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCCGAACTGGTGCGTCGTTTGCGCGATcgc                                                                                                                          >  1:536471/1‑150 (MQ=255)                              caaAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCTGAACTGGTGCGGCGTTTGCGCGATcgc                                                                                                                          >  2:370512/1‑150 (MQ=255)                                     aTTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTTCTGGCTTTAAAACCCCTTAACCTGTGGGGCGGTTTTGGGATCGCgggg                                                                                                                      >  1:692498/1‑147 (MQ=255)                                     aTTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCg                                                                                                                   <  1:705935/150‑1 (MQ=255)                                       tgtgATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGGTGCCTATAAAACCCCCGAACTGGTGCGTCGTTTGCGCGATCGCGGGGGCGa                                                                                                                  >  1:216122/1‑149 (MQ=255)                                       tgtgATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGAc                                                                                                                 <  2:163183/150‑1 (MQ=255)                                         tgATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTTTAAAAACCCCGAACTGGGGGGTCGTTTTGGCGATCGCGGGGGCGGCGTc                                                                                                              >  1:816931/1‑151 (MQ=255)                                         tgATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCTGAACTGGGGCGGCGTTTGCGCGATCGCGGGGGCg                                                                                                                   >  2:447685/1‑146 (MQ=255)                                          gATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGTTGCCTCTAAAAACCCTGAACTGGGGGGTCGTTTGCGCGATCGCGGGGGCGACGGcc                                                                                                             >  1:478165/1‑151 (MQ=255)                                                  cGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTTCTGCCTATAAAAAACCTTAAATTGTGCGTCGCTTTCGCGCTCGCGGGGCCGACGGCCGCGTGGc                                                                                                      >  1:804297/1‑150 (MQ=255)                                                  cGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGGCTATAAAAACCCCGAACTGGGGCGTCGGTTGCGCGATCGCGGGGGCGACGTCcgcg                                                                                                          >  1:913593/1‑146 (MQ=255)                                                   gCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTTCTGCCTATAAAAACCCTGAAACGGTGGGTCGGTTGCGCGGTCGCGGGGGCGACGTCCGCGGAggcc                                                                                                    >  2:46502/1‑148 (MQ=255)                                                    ccAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCAt                                                                                                   <  1:458350/151‑1 (MQ=255)                                                    ccAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGATCTCGGCGTTAGCGGCGGTATTGCTGCCCATAAAAACCCTGAACTGGGGCGGCGGTTGCGCGATCGCGGGGGCCACGTCCGCGCAGcc                                                                                                     >  2:779440/1‑149 (MQ=255)                                                     cAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATg                                                                                                  >  2:210694/1‑151 (MQ=255)                                                      aaCTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGa                                                                                                 <  2:1025731/151‑1 (MQ=255)                                                      aaCTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGGCTTGCTACCTTTAAAAACCCCGAACACGGGCGGCGTTTGCGGGATCGCCGGGCCGACATCcgcccacacctg                                                                                                  >  2:748002/1‑141 (MQ=255)                                                       aCTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGAc                                                                                                <  2:864969/151‑1 (MQ=255)                                                           cTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCCGAACTTGTGCGTCGGTTTCGCGATCGCGGGGCCGACGTCCGCGGAGCCATGACCg                                                                                              >  1:768135/1‑149 (MQ=255)                                                           cTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGGATTGGTGCCTTTAAAACCCCCTCAATTTTGCGTCGGTTGCGCCAGTGGCGGGCCGACGACcgcggtaccagggc                                                                                                >  2:122996/1‑137 (MQ=255)                                                            tGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCTGAACTGGTGCGTCGGTTGCGCGGTCGCGGGGCCGCCGTCCCCGGAGCCATGAcc                                                                                               >  1:363652/1‑147 (MQ=255)                                                             ggTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCCGAACTGGGGCGTCGTTTGCGCGGTCGCGGGGGCGACGTCCGCGTAGCCATGACCGaa                                                                                            >  1:23916/1‑149 (MQ=255)                                                              gtgtCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCg                                                                                         <  1:97184/151‑1 (MQ=255)                                                               tgtCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCTGAACTGGGGCGTCGTTTTCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCgg                                                                                        >  2:29081/1‑151 (MQ=255)                                                                  cACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCa                                                                                      <  1:272879/150‑1 (MQ=255)                                                                   aCCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCTGAACTGGGGCGGCGTTTGTGCGGTCGCGGGGGCGACGTCCCCGGAGCCAAGAcc                                                                                               >  2:826733/1‑140 (MQ=255)                                                                     ccTGCAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAg                                                                                  <  2:112433/151‑1 (MQ=255)                                                                     ccTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTTGTTCCCATAAAAACCCTCAACTTGGGCGCCGGTTTCGCGGTCCCGGGGGCCACGTCCCCGCAGCCCTGACCGAAGGGGCcaaa                                                                                   >  2:258517/1‑150 (MQ=255)                                                                     ccTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCCGAACTGGGGCGGCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCaaaa                                                                                   >  2:689159/1‑150 (MQ=255)                                                                          aGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCttt                                                                             <  2:77626/151‑1 (MQ=255)                                                                          aGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAAACCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGGCGACGTCCGCGTAGCCATGAACGCAGCGGCGAAAGCCttt                                                                             >  1:1042007/1‑151 (MQ=255)                                                                               aaaTCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCTGAACTGGTGCGTCGTTTGGGCGATCGCGGGGCGGACGTCCGCGGAGCCATGACCGAAGCGGCAAAAGCCTTTTTCAc                                                                        >  1:528462/1‑151 (MQ=255)                                                                                   catcatGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGGGCGTCGTTTTCGCGAGCGCGGGGGCGACGTCCGCGTAGCCATGACCGAAGCGGCCAAAACCTTTTTCAccccc                                                                    >  2:321713/1‑150 (MQ=255)                                                                                      catGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACtt                                                                 <  2:124912/151‑1 (MQ=255)                                                                                      catGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCCGAACTGGTGGGTCGTTTGCGCGATCGCGGGGCCGGCGTCCGCCGAGCCATGACCGAAGCGGCCAAAGCCTTTATCAcccccc                                                                   >  1:815905/1‑147 (MQ=255)                                                                                       atGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTTCTGGCCCTAAAAACCCTGAACCGGGGCGTCGTTTGGGCGATCGCGGGGGCGGCGGCCGCGTAGCCATGACCGAAGCGGCCAAAACCTTTATTACCCCAcct                                                                 >  1:927348/1‑148 (MQ=255)                                                                                         gAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAg                                                               <  1:298455/150‑1 (MQ=255)                                                                                           gCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCCATAAAAACCCCGAAATTGTGCGGCGTTTGCGCGGTCGCGGGGGCGGCGTCCCCCGAGCCCTGACCGAAGCGGCCAAAACCCTTTTCAcccccc                                                                   >  1:558788/1‑142 (MQ=255)                                                                                              tGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTTTTGCCTGTAAAAACCCCGAACTAGTGCGTCGGTTGCGCGTTTGCGGGGGCGGCGTCCGCGTAGGCATGGCCGTAAGGGCCAAAGCCTTTATCCCCCCCCCTAGcct                                                            >  2:1019451/1‑146 (MQ=255)                                                                                              tGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCCTAACTGGGGCGGCGTTTGCGCGGTCGCGGGGCCGACGGCCGCGCAGCCATGACCGGCGCGGCGAAAGCCTTTTTCACCCCCCCTAGCTTGc                                                          >  2:975725/1‑150 (MQ=255)                                                                                               ggCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTGTTTCTGCCTATAAAAAACCCGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCCTGGCCGAAGCGGCCAAAACCTTTTTCACCCCACTTAGCt                                                             >  1:734667/1‑146 (MQ=255)                                                                                               ggCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑Ag                                                       <  2:641635/151‑1 (MQ=255)                                                                                               ggCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑Ag                                                       <  2:187680/151‑1 (MQ=255)                                                                                                gCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑Agg                                                      >  2:882972/1‑151 (MQ=255)                                                                                                   ggTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑AGGC‑gg                                                  <  2:741172/151‑1 (MQ=255)                                                                                                      aaaaaaaTCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑AGGC‑GGtt                                                <  1:478047/150‑1 (MQ=255)                                                                                                       aaaaaaTCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAAACCCCCAACTGGTGGGGCGCTTTCGCGATCGCGGGGCCGACGCCCGCGGAGCCATGACCGAAGCGGCAAAAACCCTTTTCACCCCCCATAGCTTTC‑AGGG‑GGgtt                                               >  2:519061/1‑150 (MQ=255)                                                                                                        aaaaaTCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑AGGC‑GGTTTCt                                             <  2:128166/151‑1 (MQ=255)                                                                                                             tCGTTCTCGGCGTTAGCGGCGGTATTTCTGCCCCTAAAACCCCTGAACTGGTGGGGCATCTGGGCGATCGCGGGGCCGCCGTCCGCCTGGCCATGACCGCAGCGGGAAAAGCCTCTTTTACCCCACCTAGCTTCC‑AGGC‑gg                                                  >  1:618876/1‑141 (MQ=255)                                                                                                              cGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑AGGC‑GGTTTCTGGTTa                                        <  1:254025/150‑1 (MQ=255)                                                                                                               gTTCTCGGCGTTAGCGGCGGTATTGCTGGCTATAAAACCCCTGAACTGGTGCGGCGTTTTCGCGGTCGCGGGGGCCACGCCCGCGGAGCCCTGTCCGAAGCGGCCAAAACCCTTTTAACCCCCCTTAGCTTTCAAGGG‑GGTTTTTGtttt                                        >  2:120202/1‑149 (MQ=255)                                                                                                                ttCTCGGCGTTAGCGGCGGTATTGCTGCCCATAAAAACCCCGAACTGGTGCGGCGTTTTCGCGATCGCGGGGGCGACGTCCGCGTAGGCATGACCGAAGCGGCAAAAGCCTTTTTCACCCCACCTAGCTTGC‑AGGC‑GGgtttttg                                           >  2:646089/1‑141 (MQ=255)                                                                                                                  ctcGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCCGAACTGGGGCGTCGTTTTCGCGATCGCGGGGCCGACGGCCGCGGAGCCATGACCGAAGCGGCAAAAGCCTTTTTCACCCCCCTTAGCTTGC‑AGGCGGGTTTTTGGTTTTcc                                     >  2:784961/1‑150 (MQ=255)                                                                                                                     ggCGTTAGCGGCGGTATTTGTGCCCCTAAAACCCCCTAAATGGTGCGGCGTTTGGGCGGTCGCGGGGGCGACGTCCGCGTAGCCATGACCGAAGAGGCCAAAGCCTTTTTCACCCCCCCTAGCTTGC‑CGGCGGGTTTTTGGTTTTCCGGtt                                 >  1:725266/1‑151 (MQ=255)                                                                                                                     ggCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑AGGC‑GGTTTCTGGTTATCCCGttt                                <  2:195786/151‑1 (MQ=255)                                                                                                                      gCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑AGGC‑GGTTTCTGGTTATCCCGTTTc                               <  1:657994/151‑1 (MQ=255)                                                                                                                           aGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑AGGC‑GGTTTCTGGTTATCCCGTTTCCGACa                          <  1:46497/151‑1 (MQ=255)                                                                                                                               gcggTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑AGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTct                      <  2:662871/151‑1 (MQ=255)                                                                                                                               gcggTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑AGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTct                      <  2:258499/151‑1 (MQ=255)                                                                                                                                 ggTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑AGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTctgc                    >  2:343315/1‑151 (MQ=255)                                                                                                                                      tgctgcCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑AGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGa                >  1:9474/1‑150 (MQ=255)                                                                                                                                          gcCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑AGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCg            <  1:43316/150‑1 (MQ=255)                                                                                                                                           ccTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑AGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCgg           >  1:823259/1‑150 (MQ=255)                                                                                                                                              ataAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑AGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCa         <  2:817065/149‑1 (MQ=255)                                                                                                                                              ataAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑AGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCAg        >  2:41471/1‑150 (MQ=255)                                                                                                                                                 aaaCCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑AGGC‑GGTTTCTGGGTATCCCGTTTCCGACAGTCTGCTGGACCCGGc          >  1:1032280/1‑145 (MQ=255)                                                                                                                                                  aaCCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑AGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCAGCCGaa   <  2:768263/151‑1 (MQ=255)                                                                                                                                                   aCCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑AGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCAGCCGAAg  <  2:674034/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                        AAAGTGAATCTCAGGTGAAGGCACGACAACAAAAGTATTGTGATAAAATCGCCAACTTCTGGTGTCACCCTACAGGAAAAATCATCATGAGCCTGGCCGGTAAAAAAATCGTTCTCGGCGTTAGCGGCGGTATTGCTGCCTATAAAACCCCTGAACTGGTGCGTCGTTTGCGCGATCGCGGGGCCGACGTCCGCGTAGCCATGACCGAAGCGGCAAAAGCCTTTATCACCCCACTTAGCTTGC‑AGGC‑GGTTTCTGGTTATCCCGTTTCCGACAGTCTGCTGGACCCGGCAGCCGAAG  >  NC_000913/3812645‑3812941

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