BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,032,404	G→C	5.2%	M346I (ATG→ATC) ‡	uacT →	uric acid permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,032,404	0	G	C	5.2%	355.3 / 9.1	116	M346I (ATG→ATC) ‡	uacT	uric acid permease
Reads supporting (aligned to +/- strand): ref base G (57/53); new base C (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 = 4.87e-01

Only 100 of 118 total aligned reads displayed.

GTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATA  >  NC_000913/3032255‑3032549

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

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

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

    aTGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGc                                                                                                                                               <  2:1344984/150‑1 (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)

           aGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAAt                                                                                                                                        >  2:1154470/1‑150 (MQ=255)

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

            gggACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATg                                                                                                                                       >  2:1716103/1‑150 (MQ=255)

              gACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGc                                                                                                                                     >  1:336007/1‑150 (MQ=255)

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

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

                cAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCgg                                                                                                                                   >  1:1561598/1‑150 (MQ=255)

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

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

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

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

                               cGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCa                                                                                                                    >  1:613273/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                                                                                                                >  2:646580/1‑150 (MQ=255)

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

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

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

                                     aCAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGc                                                                                                              <  2:459593/150‑1 (MQ=255)

                                      cAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCa                                                                                                             >  2:159913/1‑150 (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:1273357/1‑150 (MQ=255)

                                                   acGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGggc                                                                                                <  1:1154470/150‑1 (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:808065/1‑150 (MQ=255)

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

                                                               cAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTa                                                                                    >  2:748780/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)

                                                                   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)

                                                                                    gTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTg                                                               >  1:1665757/1‑150 (MQ=255)

                                                                                      gACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGc                                                             <  2:200555/150‑1 (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)

                                                                                              ttCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGa                                                     <  1:78883/150‑1 (MQ=255)

                                                                                                  tAGCCGCTGGGTGTGTAGTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCg                                                 <  2:711631/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                                  <  1:714844/150‑1 (MQ=255)

                                                                                                                 tATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTg                                  <  2:1331859/150‑1 (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:989337/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)

                                                                                                                                      aaTCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGATTTCTGTCGCGCTGTAACGACACCACCAACCGTTa             >  2:345187/1‑150 (MQ=255)

                                                                                                                                       gccatttttGGCACCATGCCAAACATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACAACCAACCGTTAc            >  1:843139/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)

GTCGATGGCGTAGGGACAATGATAGGCGGCACGTTTAACAGCTTCCCCCACACGTCATTTTCTCAAAACGTTGGCCTGGTTAGCGTGACGCGCGTTCATAGCCGCTGGGTGTGTATTTCTTCGGGAATTATATTAATCCTGTTTGGCATGGTGCCAAAAATGGCGGTGCTGGTAGCCTCCATTCCGCAATTTGTGCTGGGCGGCGCTGGTCTAGTGATGTTCGGCATGGTACTGGCGACAGGGATTCGAATTCTGTCGCGCTGTAACTACACCACCAACCGTTACAACCTCTATA  >  NC_000913/3032255‑3032549

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