BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,032,405	G→A	5.1%	V347M (GTG→ATG)	uacT →	uric acid permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,032,405	0	G	A	5.1%	368.0 / 5.6	116	V347M (GTG→ATG)	uacT	uric acid permease
Reads supporting (aligned to +/- strand): ref base G (57/53); new base A (4/2); total (61/55)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.39e-01

Only 100 of 117 total aligned reads displayed.

TCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATA  >  NC_000913/3032256‑3032549

tCGATGGCGTATGGACAATGATAGGAGGCACTTTTATCAGCTTGCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATgg                                                                                                                                                  <  2:275989/150‑1 (MQ=255)

 cGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGt                                                                                                                                                 >  2:1026920/1‑150 (MQ=255)

 cGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGt                                                                                                                                                 >  2:1608537/1‑150 (MQ=255)

  gATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTg                                                                                                                                                >  2:86807/1‑150 (MQ=255)

     ggCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCa                                                                                                                                             <  1:1710012/150‑1 (MQ=255)

     ggCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCa                                                                                                                                             >  1:124708/1‑150 (MQ=255)

         tAGGGACAATGGTAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCACCATGCCaaaca                                                                                                                                         <  2:453704/150‑3 (MQ=255)

          aGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAACCGTTTGCCTGGTTAGCGTGACGCGCGTTCATAGCAGCCGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAAt                                                                                                                                        <  1:159030/150‑1 (MQ=255)

          aGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCACCATGCCaaacag                                                                                                                                        <  1:1315424/150‑4 (MQ=255)

             gACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGc                                                                                                                                     <  1:1131539/150‑1 (MQ=255)

             gACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGGTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTAAAATAATCCTGTTTGGCATGGTGCCAAAAATgga                                                                                                                                     >  2:311098/1‑149 (MQ=255)

                   gATAGGCGGCACGTTTAACATCTTCGCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGTGTGTGTATTTCTTCGGTAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCt                                                                                                                               <  1:1204728/150‑1 (MQ=255)

                   gATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCt                                                                                                                               <  1:1555460/150‑1 (MQ=255)

                   gATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCt                                                                                                                               <  2:385845/150‑1 (MQ=255)

                       ggcggcACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCATAAATGGCGGTGCTGGTa                                                                                                                           >  1:96820/1‑150 (MQ=255)

                          ggcACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGcc                                                                                                                        >  2:388966/1‑150 (MQ=255)

                            cACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTc                                                                                                                      >  1:1331859/1‑150 (MQ=255)

                               gTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCAt                                                                                                                   <  1:1585340/150‑1 (MQ=255)

                                 ttAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGTCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGAGTGTATTTCTTCGGGAATTATATTACTCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTc                                                                                                                 >  2:425616/1‑150 (MQ=255)

                                  tAACAGCTTCCCCCACACTTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTGCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTcc                                                                                                                <  2:303651/150‑1 (MQ=255)

                                  tAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTcc                                                                                                                >  1:584348/1‑150 (MQ=255)

                                  tAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTcc                                                                                                                >  2:646580/1‑150 (MQ=255)

                                  tAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTcc                                                                                                                >  1:200555/1‑150 (MQ=255)

                                   aaCAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCg                                                                                                               >  2:1521799/1‑150 (MQ=255)

                                   aaCAGAGTCCCCCACACGTAATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATGTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCg                                                                                                               <  2:72056/150‑1 (MQ=255)

                                           cccccACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTtgt                                                                                                       >  2:78883/1‑150 (MQ=255)

                                              ccACACGTCATTTTCTCAAAACGTTGGCCTGGATAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGAATGGTGCCAAAAATGGCGGTGCTGGTAGACTCCATTCCGCAATTTGTGCt                                                                                                    >  1:851617/1‑150 (MQ=255)

                                               cacaCGTCATTTGCTCAAAACGTTGGCCTGGGTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTGCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTg                                                                                                   <  2:1030602/150‑1 (MQ=255)

                                                 cacGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTggg                                                                                                 >  2:457611/1‑150 (MQ=255)

                                                  acGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGggc                                                                                                <  1:1154470/150‑1 (MQ=255)

                                                  acGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGggc                                                                                                >  1:1273357/1‑150 (MQ=255)

                                                   cGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGggcg                                                                                               >  2:1438428/1‑150 (MQ=255)

                                                        ttttCTCAAAACGTTGGCCTGGGTAGCGTGACGCGCGTTCATAGCCGCTGCGTGTGTATAGCTTCGGGAATTATATTCATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCt                                                                                          <  2:1685602/150‑1 (MQ=255)

                                                         tttCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCta                                                                                         >  2:624232/1‑149 (MQ=255)

                                                          ttCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTgg                                                                                        <  1:1026920/150‑1 (MQ=255)

                                                            ctcAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTAGGCATGGTGCCAAAAATGGCGGTGCTGGTAGACTCCATTCCGAAATTTGTGCTGTGCGGCGCTGGTc                                                                                      >  1:824357/1‑150 (MQ=255)

                                                             tcAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCt                                                                                     >  1:1491317/1‑150 (MQ=255)

                                                             tcAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGAGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCt                                                                                     >  2:113476/1‑150 (MQ=255)

                                                              cAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTa                                                                                    >  2:577402/1‑150 (MQ=255)

                                                              cAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTa                                                                                    >  2:748780/1‑150 (MQ=255)

                                                              cAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTa                                                                                    >  2:808065/1‑150 (MQ=255)

                                                               aaaaCGTTGGCCTGGTAAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCATAAATGGCGGTGCTGGTAGCCTCCATTCCGCATTTTGTGCTGTGCGGCGCTGGTCTAg                                                                                   <  1:979762/150‑1 (MQ=255)

                                                                aaaCGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGt                                                                                  >  1:649522/1‑150 (MQ=255)

                                                                aaaCGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGCGTGTGTATTACTTCGGGAATTATATTAATCCTGTTTGGCATGGTCCCAAAAATGGCGGTCCTGGTAGCCTCCATTCCGCAATTTGTGATGTGCGGCGCTGGTATAGt                                                                                  >  1:167986/1‑150 (MQ=255)

                                                                 aaCGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTg                                                                                 <  2:573469/150‑1 (MQ=255)

                                                                  aCGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGa                                                                                <  1:1183146/150‑1 (MQ=255)

                                                                    gTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGGTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCTTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATg                                                                              <  2:1455765/150‑1 (MQ=255)

                                                                       ggCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTc                                                                           >  1:1690278/1‑150 (MQ=255)

                                                                             gTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATg                                                                     <  2:1561598/150‑1 (MQ=255)

                                                                              ttAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATgg                                                                    <  1:456188/150‑1 (MQ=255)

                                                                              ttAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATgg                                                                    <  2:1639542/150‑1 (MQ=255)

                                                                                     gACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGc                                                             <  2:200555/150‑1 (MQ=255)

                                                                                        gcgcgTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCAAAAAATGCCGGTGCTGGTAGCATCCTTTCCGAAATTTGTCATGGGCGGAGCTGGTCTATTGTTGTTCGGATTGGTACTGGCGaa                                                          >  1:725353/1‑149 (MQ=255)

                                                                                           cgTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAgg                                                       <  1:1608537/150‑1 (MQ=255)

                                                                                           cgTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAgg                                                       <  2:401049/150‑1 (MQ=255)

                                                                                            gTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTTATCCTGTTTGGCATGGTGACAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGATGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAggg                                                      >  1:451656/1‑150 (MQ=255)

                                                                                             ttCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGa                                                     >  2:714844/1‑150 (MQ=255)

                                                                                                 tAGCCGCTGGGTGTGTAGTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCg                                                 <  2:711631/150‑1 (MQ=255)

                                                                                                      gCTGGGTGTGTATTTCTTCGGTAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTc                                            <  2:1628720/150‑1 (MQ=255)

                                                                                                        tGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGgcgc                                                                                           >  2:1431297/1‑101 (MQ=255)

                                                                                                        tGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGgcgc                                                                                           <  1:1431297/101‑1 (MQ=255)

                                                                                                        tGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTg                                          <  1:985296/150‑1 (MQ=255)

                                                                                                          ggTGTGTATTTCTTCGTGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCTCTGGTCTAGTGATGTTCTGAATGGTACTGGCGACAGGGATTCGAATTCTGTc                                        >  2:1110167/1‑150 (MQ=255)

                                                                                                          ggTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTc                                        >  1:363144/1‑150 (MQ=255)

                                                                                                           gtgtgtATTTCTTCGGGAATTATCTTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTcg                                       <  1:1086628/150‑1 (MQ=255)

                                                                                                           gtgtgtATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTcg                                       <  1:388966/150‑1 (MQ=255)

                                                                                                            ggtgtATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTcgc                                      <  1:159913/149‑1 (MQ=255)

                                                                                                                tATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTg                                  <  2:1331859/150‑1 (MQ=255)

                                                                                                                tATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTg                                  <  1:714844/150‑1 (MQ=255)

                                                                                                                    tcttcGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAAc                              >  2:1257125/1‑150 (MQ=255)

                                                                                                                    tcttcGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGGCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAAc                              <  2:921605/150‑1 (MQ=255)

                                                                                                                     cttcGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACt                             >  1:1600709/1‑150 (MQ=255)

                                                                                                                     cttcGGGAATTATATTAATCCTGTTTGGCAGGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACt                             <  1:1463374/150‑1 (MQ=255)

                                                                                                                     ctgagggAATTATATTAATCCTGTTTGGCATGGGGCCAAATATGGCGGTGCTGGTATCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGGCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACt                             <  1:311098/146‑1 (MQ=255)

                                                                                                                       tcGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTac                           >  1:679677/1‑150 (MQ=255)

                                                                                                                       tcGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTATGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTac                           >  1:1062098/1‑150 (MQ=255)

                                                                                                                         gggAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTacac                         >  1:553073/1‑150 (MQ=255)

                                                                                                                          ggAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGGGCTGGTAGCCTCCATTCCGCTATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTAcacc                        <  2:124708/150‑1 (MQ=255)

                                                                                                                           gAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTAcacca                       <  1:748780/150‑1 (MQ=255)

                                                                                                                           gAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCGGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTAcacca                       <  1:1057486/150‑1 (MQ=255)

                                                                                                                             aTTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACaccacc                     <  2:1388102/150‑1 (MQ=255)

                                                                                                                                atatTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAAc                  >  2:1200386/1‑150 (MQ=255)

                                                                                                                                 aatTAATCCAGTTTGGCATGGTGCCAAAAATGGCGTTGCAGGTAGCCTCAATTCCGCAATATTAGCTGTGCGGCGATGGTCTAGTGAGGTTCGGCATGGTACTGGCGACCGGGATTCGAATTATCTCGCGCTGGAACTACACCTCCAAca                 >  1:1078993/2‑149 (MQ=255)

                                                                                                                                  atTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCg                >  2:308459/1‑150 (MQ=255)

                                                                                                                                   ttACTCCGGTTCGGCAAGGGGCCAAAAATGGCGGGGCTGGTAGCCTCCATTCCGAAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGt               <  2:846573/150‑1 (MQ=255)

                                                                                                                                     aaTCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTa             >  1:1559577/1‑150 (MQ=255)

                                                                                                                                      gccatttttGGCACCATGCCAAACATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACAACCAACCGTTAc            >  1:843139/7‑150 (MQ=255)

                                                                                                                                      gccatttttGGCACCATGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTATCTACACCACCAACCGTTAc            >  1:1626280/7‑150 (MQ=255)

                                                                                                                                      gccatttttGGCACCATGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTAc            >  2:1293807/7‑150 (MQ=255)

                                                                                                                                      aTCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTAc            <  2:1690278/150‑1 (MQ=255)

                                                                                                                                        catttttGGCACCATGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACaa          >  2:1689617/5‑150 (MQ=255)

                                                                                                                                           gTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACct       >  2:946870/1‑150 (MQ=255)

                                                                                                                                             ttGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACctct     <  2:96820/150‑1 (MQ=255)

                                                                                                                                             ttGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACctct     <  2:1138236/150‑1 (MQ=255)

                                                                                                                                             ttGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACctct     >  2:1293895/1‑150 (MQ=255)

                                                                                                                                              tGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCta    <  1:1716103/150‑1 (MQ=255)

                                                                                                                                              tGGCATGGTGCCAAAAATGGCGGGGCTGGTAGCCTCCATTCCGCAATTTGTGCGGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCta    <  2:684767/150‑1 (MQ=255)

                                                                                                                                              tGGCATGGGGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCta    <  2:336007/150‑1 (MQ=255)

                                                                                                                                               ggCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCtat   <  2:613273/150‑1 (MQ=255)

                                                                                                                                                gCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTTTCGCGCTGTAACTACACCACCAACCGTTACAACCTCtata  <  1:728724/150‑1 (MQ=255)

TCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATA  >  NC_000913/3032256‑3032549

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