Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,008,075	A→G	D189G (GAC→GGC)	ygeX →	2,3‑diaminopropionate ammonia lyase, PLP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,008,075	0	A	G	98.6%	232.7 / NA	71	D189G (GAC→GGC)	ygeX	2,3‑diaminopropionate ammonia lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (46/23);  total (48/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ACATGCCGAAAGGTTCTGCTCAGGAACGCGTTGACGCCATTCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGCGTAACCCCGACGCACGTTCTGCTGCAAGCCGGTGTCGGAGCAATGG  >  NC_000913/3007931‑3008218                                                                                                                                                 |                                                                                                                                                aCATGCCGAAAGGTTCTGCTCAGGAACGCGTTGACGCCATTCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGc                                                                                                                                           <  1:1233968/151‑1 (MQ=255) aCATGCCGAAAGGTTCTGCTCAGGAACGCGTTGACGCCATTCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCGGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGc                                                                                                                                           <  2:561866/151‑1 (MQ=255)     gCCGAAAGGTTCTGCTCAGGAACGCGTTGACGCCCTTCTGAACCTCGGCGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                                                                       >  2:379426/1‑151 (MQ=255)         aaaGGTTCTGCTCAGGAACGCGTTGACGCCATTCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGaa                                                                                                                                    <  2:186872/150‑1 (MQ=255)         aaaGGTTCTGCTCAGGAACGCGTTGACGCCATTCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGaa                                                                                                                                    <  2:594114/150‑1 (MQ=255)         aaaGGTTCTGCTCAGGAACGCGTTGACGCCATTCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGaa                                                                                                                                    <  2:593617/150‑1 (MQ=255)         aaaGGTTCTGCTCAGGAACGCGTTGACGCCATTCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAg                                                                                                                                   <  1:1245472/151‑1 (MQ=255)          aaGGTTCTGCTCAGGAACGCGTTGACGCCATTCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAgg                                                                                                                                  <  2:1051920/151‑1 (MQ=255)           aGGTTCTGCTCAGGAACGCGTTGACGCCATTCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGt                                                                                                                                 >  1:1182498/1‑151 (MQ=255)             gTTCTGCTCAGGAACGCGTTGACGCCATTCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGt                                                                                                                                 <  2:845881/149‑1 (MQ=255)             gTTCTGCTCAGGAACGCGTTGACGCCATTCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGt                                                                                                                                 <  1:442038/149‑1 (MQ=255)             gTTCTGCTCAGGAACGCGTTGACGCCATTCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTa                                                                                                                               >  2:859693/1‑151 (MQ=255)                  gCTCAGGAACGCGTTGACGCCATTCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCa                                                                                                                          >  1:1067134/1‑151 (MQ=255)                  gCTCAGGAACGCGTTGACGCCATTCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGGTACACCa                                                                                                                          >  2:304867/1‑151 (MQ=255)                         aaCGCGTTGACGCCATTCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATcc                                                                                                                    <  1:601448/150‑1 (MQ=255)                             cgTTGACGCCATTCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAAcc                                                                                                               <  2:295102/151‑1 (MQ=255)                              gTTGACGCCATTCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAAcc                                                                                                               <  2:1225149/150‑1 (MQ=255)                                  aCGCCATTCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGAt                                                                                                          <  1:974546/151‑1 (MQ=255)                                  aCGCCATTCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGAt                                                                                                          <  1:491149/151‑1 (MQ=255)                                  aCGCCATTCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGAt                                                                                                          <  1:151182/151‑1 (MQ=255)                                       aTTCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATg                                                                                                      >  1:942720/1‑150 (MQ=255)                                         tCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGcga                                                                                                   >  2:871258/1‑149 (MQ=255)                                          cTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAg                                                                                                  >  2:196395/1‑151 (MQ=255)                                          cTGAACCTCGGTGCCGAGTGCATCGGCACGGATATGAACTATGACGATACCGTTCGCCTGACAATGCAACACGCGCAGGAGCACGGCGGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAg                                                                                                  <  1:576173/151‑1 (MQ=255)                                               ccTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTa                                                                                              >  1:813614/1‑150 (MQ=255)                                                    gTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCaa                                                                                         >  2:975874/1‑150 (MQ=255)                                                      gCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAAcc                                                                                       >  1:685571/1‑150 (MQ=255)                                                           atGCATCGTCCCGGATATGAACTATGTCGAGACCGTACGCCTGACCATGCAACACGCGTAGCAGCACGGCTGGGAAGGGGTACAGCACACGGCATGGGAAGGTTACACCAAAATCCCAACCGGGATCATGCCAGGCTATGCAACCGTGGCa                                                                                 >  2:37372/2‑151 (MQ=255)                                                            tGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGc                                                                                  <  2:685262/149‑1 (MQ=255)                                                                tCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGa                                                                            >  1:165402/1‑151 (MQ=255)                                                                  gTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAg                                                                          >  1:422389/1‑151 (MQ=255)                                                                    cACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGcc                                                                        <  1:1115398/151‑1 (MQ=255)                                                                      cGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGCTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGt                                                                      >  1:103171/1‑151 (MQ=255)                                                                        gATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCg                                                                    >  1:1158390/1‑151 (MQ=255)                                                                        gATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCg                                                                    >  1:122239/1‑151 (MQ=255)                                                                        gATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCg                                                                    >  2:960687/1‑151 (MQ=255)                                                                        gATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCg                                                                    >  2:1136875/1‑151 (MQ=255)                                                                        gATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCg                                                                    >  2:121631/1‑151 (MQ=255)                                                                        gATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCg                                                                    >  1:251206/1‑151 (MQ=255)                                                                             gAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGCCGAGCaa                                                               >  1:37408/1‑151 (MQ=255)                                                                             gAAATATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGACACACCAAAATCCCAAGCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCaa                                                               >  1:1059084/1‑151 (MQ=255)                                                                              aaCTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCaaa                                                              >  1:1635/1‑151 (MQ=255)                                                                                cTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATg                                                            <  2:938606/151‑1 (MQ=255)                                                                                   tGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCgg                                                         >  2:340341/1‑150 (MQ=255)                                                                                   tGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGt                                                         >  1:484678/1‑151 (MQ=255)                                                                                    gACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTg                                                        >  2:1202005/1‑151 (MQ=255)                                                                                        aTACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAAt                                                    >  1:444602/1‑151 (MQ=255)                                                                                          aCCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATgg                                                  >  2:961706/1‑151 (MQ=255)                                                                                           ccGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCCAAATCCAAACCTGGATCATGCCAGGCTACGCAACCCTGGCAGCTGCAGCCGTCGAGCAAATGCGTGAAATggg                                                 >  2:774925/1‑151 (MQ=255)                                                                                            cGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGc                                                <  1:385934/151‑1 (MQ=255)                                                                                                cGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGCGTaa                                            >  2:1775/1‑151 (MQ=255)                                                                                                   cTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGCGTAAccc                                         <  1:897259/151‑1 (MQ=255)                                                                                                      aCCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGCGAAATGGGCGTAACCCCGa                                      >  2:64974/1‑151 (MQ=255)                                                                                                         aTGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGCGTAACCCCGACGc                                   <  1:894268/151‑1 (MQ=255)                                                                                                          tGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGCGTAACCCCGACGCa                                  >  1:708790/1‑151 (MQ=255)                                                                                                                acGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGAGAAATGGGCGTAACCCCGAAGCACGtt                              >  2:677893/1‑149 (MQ=255)                                                                                                                     cagcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGCGTAACCCCGACGCACGTTctgctg                        >  1:389128/1‑150 (MQ=255)                                                                                                                         agcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGCGTAACCCCGACGCACGTTCTGCTGCAAg                    >  2:498512/1‑150 (MQ=255)                                                                                                                             ctgCTGGGAAGTGGTAGAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTCCGCAACCCTGGCCGAAGAAGCCCTCGAGCCAATGCGTGAAATGAGCGTAACCGCGACGCACGTTCTGCTGCGAGCCggg               >  2:53866/3‑150 (MQ=255)                                                                                                                             cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGCGTAACCCCGACGCCAGTTCTGCTGCAAGCCGgt               >  2:448260/1‑151 (MQ=255)                                                                                                                               gCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGCGTAACCCCGACGCACGTTCTGCTGCAAGCCGgtgt             >  1:39279/1‑151 (MQ=255)                                                                                                                                  gggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGCGTAACCCCGACGCACGTTCTGCTGCAAGCCGGTGTCgg          >  2:24899/1‑151 (MQ=255)                                                                                                                                    gAAGCGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGCGTAACCCCGACGCACGTTCTGCTGCAAGCCGGTGTCGGAg        >  2:631811/1‑151 (MQ=255)                                                                                                                                       gTGGTACAGGTCACGGCATGGGAAGGTTACACAAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGCGTAACCCCGACGCACGTTCTGCTGCACGCCGGTGTCGGAGCaa     >  1:126955/1‑151 (MQ=255)                                                                                                                                       gTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGCGTAACCCCGACGCACGTTCTGCTGCAAGCCGGTGTCGGAGaaa     >  2:960919/1‑151 (MQ=255)                                                                                                                                       gTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGCGTAACCCCGACGCACGTTCTGCTGCAAGCCGGTGTCGGAGCaa     >  1:910704/1‑151 (MQ=255)                                                                                                                                       gTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGCGTAACCCCGACGCACGTTCTGCTGCAAGCCGGTGTCGGAGCaa     >  2:413241/1‑151 (MQ=255)                                                                                                                                       gTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGCGTAACCCCGACGCACGTTCTGCTGCAAGCCGGTGTCGGAGCaa     >  1:641742/1‑151 (MQ=255)                                                                                                                                         ggTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGCGTAACCCCGACGCACGTTCTGCTGCAAGCCGGTGTCGGAGCAATg   >  1:1050621/1‑151 (MQ=255)                                                                                                                                          gTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGCGTAACCCCGACGCACGTTCTGCTGCAAGCCGGTGTCGGAGCAATgg  >  1:127824/1‑151 (MQ=255)                                                                                                                                           tACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGCGTAACCCCGACGCACGTTCTGCTGCAAGCCGGTGTCGGAGCAATgg  <  1:471750/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                ACATGCCGAAAGGTTCTGCTCAGGAACGCGTTGACGCCATTCTGAACCTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGCGTAACCCCGACGCACGTTCTGCTGCAAGCCGGTGTCGGAGCAATGG  >  NC_000913/3007931‑3008218

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