Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,933,239	C→T	9.4%	D456N (GAC→AAC)	edd ←	phosphogluconate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,933,239	0	C	T	9.4%	183.1 / 7.7	64	D456N (GAC→AAC)	edd	phosphogluconate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (28/29);  new base T (2/4);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

CAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCA  >  NC_000913/1933092‑1933388                                                                                                                                                    |                                                                                                                                                      cAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGAAGGCATAACGTCAt                                                                                                                                                     >  1:894178/1‑150 (MQ=255)  aaTTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATg                                                                                                                                                    <  1:1221293/150‑1 (MQ=255)   aTTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGc                                                                                                                                                   <  2:853045/150‑1 (MQ=255)   aTTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGc                                                                                                                                                   >  1:1618618/1‑150 (MQ=255)     tttGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTg                                                                                                                                                 >  2:288808/1‑150 (MQ=255)        gAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCt                                                                                                                                              <  1:464194/150‑1 (MQ=255)           acacCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCAGCATGACGTTATGCCGGCcttt                                                                                                                                           <  2:469202/150‑2 (MQ=255)           acacCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCAGCATGACGTTATGCCGGCcttt                                                                                                                                           <  2:1538755/150‑2 (MQ=255)             acCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCaa                                                                                                                                         >  2:980911/1‑150 (MQ=255)                 gTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAaac                                                                                                                                     <  2:806328/150‑1 (MQ=255)                    cAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTTATGCTGGCTTTCAAAaacaac                                                                                                                                  >  2:822168/1‑150 (MQ=255)                                aaGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTc                                                                                                                      <  2:1313007/150‑1 (MQ=255)                                  gTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCaa                                                                                                                    <  2:927549/150‑1 (MQ=255)                                  gTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCAGCATGACGTTATGCCGGCctttg                                                                                                                                          <  2:1540265/128‑3 (MQ=255)                                  gTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCAGCATGACGTTATGCCGGCctttg                                                                                                                                          >  1:1540265/1‑126 (MQ=255)                                     gcggcATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCa                                                                                                                 <  1:887684/150‑1 (MQ=255)                                       ggcATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCAcc                                                                                                               <  2:542575/150‑1 (MQ=255)                                          aTGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTgg                                                                                                            <  1:705173/150‑1 (MQ=255)                                              gTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTtct                                                                                                        <  2:988840/150‑1 (MQ=255)                                                 tATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCaa                                                                                                     >  2:828051/1‑150 (MQ=255)                                                  aTGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAAc                                                                                                    >  2:681481/1‑150 (MQ=255)                                                     tAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTTATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACcgg                                                                                                 <  1:822168/150‑1 (MQ=255)                                                          cTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACcggcacgg                                                                                            >  2:455008/1‑150 (MQ=255)                                                            ggCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCACCcggcacggca                                                                                          >  1:543978/1‑150 (MQ=255)                                                            ggCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACcggcacggca                                                                                          <  1:94745/150‑1 (MQ=255)                                                               aTTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCagag                                                                                       >  2:478118/1‑150 (MQ=255)                                                                  ccGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGtt                                                                                    >  2:1511120/1‑150 (MQ=255)                                                                  ccGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGtt                                                                                    >  2:1396886/1‑150 (MQ=255)                                                                     ttCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCTGTCCAGCAAACCCGCTTCAAAGGCCGGCATATCGTCATGCTTGCTTTCAAAAAGAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTc                                                                                 <  2:1669926/150‑1 (MQ=255)                                                                      tCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCa                                                                                <  1:34817/150‑1 (MQ=255)                                                                      tCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCa                                                                                <  1:726321/150‑1 (MQ=255)                                                                              ggCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCa                                                                        >  1:281107/1‑150 (MQ=255)                                                                                    tgatgaCGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGccc                                                                  <  1:980911/150‑1 (MQ=255)                                                                                     gatgaCGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCa                                                                 >  1:78373/1‑150 (MQ=255)                                                                                      atgaCGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAg                                                                <  1:681481/150‑1 (MQ=255)                                                                                       tgaCGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAgg                                                               >  1:926478/1‑150 (MQ=255)                                                                                        gaCGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCACCCGGCCGGGCAGAGGTTTTCATAACCGCACGGCCCAGGt                                                              >  2:898433/1‑150 (MQ=255)                                                                                         aCGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGtt                                                             >  1:1217413/1‑150 (MQ=255)                                                                                            gACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCGTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTAcc                                                          >  2:1287028/1‑150 (MQ=255)                                                                                                acaacGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCtt                                                      >  2:1348836/1‑150 (MQ=255)                                                                                                 caacGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTa                                                     >  1:1668172/1‑150 (MQ=255)                                                                                                       cacaATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACtt                                               >  2:581493/1‑150 (MQ=255)                                                                                                       cacaATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACtt                                               >  2:1283652/1‑150 (MQ=255)                                                                                                        acaATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACttt                                              <  1:288808/150‑1 (MQ=255)                                                                                                        acaATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACttt                                              <  1:478118/150‑1 (MQ=255)                                                                                                        acaATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACttt                                              >  1:771960/1‑150 (MQ=255)                                                                                                          aaTCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTg                                            <  1:1511120/150‑1 (MQ=255)                                                                                                             cgcgGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTcc                                         <  1:970750/150‑1 (MQ=255)                                                                                                                ggTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGGCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCccac                                      <  2:1195610/150‑1 (MQ=255)                                                                                                                    cAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCCCCatg                                  >  1:834298/1‑150 (MQ=255)                                                                                                                     aGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCatga                                 >  1:1419854/1‑150 (MQ=255)                                                                                                                       cAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCatgatg                               <  2:254096/150‑1 (MQ=255)                                                                                                                       cAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCatgatg                               <  1:1305099/150‑1 (MQ=255)                                                                                                                         aaCCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATgag                             <  2:909973/150‑1 (MQ=255)                                                                                                                           cccGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATgagag                           >  2:1305222/1‑150 (MQ=255)                                                                                                                           cccGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTACTCAACCGGCACGGCTGAGGTTTTCATATGCGCACGGCCCAGGTTACCGCTTAACACTTTTGTGCCACCATGATgagag                           <  2:99639/150‑1 (MQ=255)                                                                                                                           cccGCTTCAAAGGCCGGCATAACGTAATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCTCCTGGTTCTCAACCGGCACGGGAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATgagag                           <  2:253293/150‑1 (MQ=255)                                                                                                                                   aaaGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGgtgg                   <  2:889458/150‑3 (MQ=255)                                                                                                                                    aaGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGgttgt                  >  2:1368753/1‑150 (MQ=255)                                                                                                                                    aaGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGgttgt                  <  1:1580521/150‑1 (MQ=255)                                                                                                                                    aaGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGgttgt                  >  1:1166473/1‑150 (MQ=255)                                                                                                                                           gCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCgaagg           <  2:926478/150‑1 (MQ=255)                                                                                                                                              tAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCgaaggaag        >  1:640984/1‑150 (MQ=255)                                                                                                                                                   tCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCTATc   <  2:1419854/150‑1 (MQ=255)                                                                                                                                                    cATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTg                                            >  2:1543664/1‑108 (MQ=255)                                                                                                                                                    cATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTg                                            <  1:1543664/108‑1 (MQ=255)                                                                                                                                                    cATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATca  <  1:531922/150‑1 (MQ=255)                                                                                                                                                    cATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATca  <  1:1209802/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CAATTTTGAAACACCGGTCCAATAATACACCAAGTGGCGGCATGAGTTTATGTAATTCTGGCATTCCGTTCGCTTTTGGCCCCTGATGACGGACAACAACGACACAATCGCGGTCCAGCAAACCCGCTTCAAAGGCCGGCATAACGTCATGCTGGCTTTCAAAAACAACCGCTGGCGCTTCAATCACCTGGTTCTCAACCGGCACGGCAGAGGTTTTCATAACCGCACGGCCCAGGTTACCGCTTAACACTTTTGTCCCACCATGATGAGAGAAAGGTTGTTCGAAGGAAGCGATCA  >  NC_000913/1933092‑1933388

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