BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,285,528	A→C	5.3%	L615R (CTG→CGG)	acs ←	acetyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,285,528	0	A	C	5.3%	407.6 / 8.4	130	L615R (CTG→CGG)	acs	acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand): ref base A (58/65); new base C (5/2); total (63/67)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 131 total aligned reads displayed.

AACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCGTGATTA  >  NC_000913/4285380‑4285677

aaCCGTTACAGACTCGCATCGGGCAATAGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTACAAGCCGCTTATCGACTACGCCAGGATCGGCATGCCTCGACGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAg                                                                                                                                                      >  2:507667/1‑150 (MQ=255)

 aCCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGa                                                                                                                                                     <  2:566138/150‑1 (MQ=255)

   cGTTACCGACTCGCATCGGGCAATTGGGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGGTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAAt                                                                                                                                                   <  2:1043998/150‑1 (MQ=255)

     ttACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATAc                                                                                                                                                 <  1:234569/150‑1 (MQ=255)

     ttACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCAGCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATAc                                                                                                                                                 <  1:421385/150‑1 (MQ=255)

       aCCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACgg                                                                                                                                               <  2:1214484/150‑1 (MQ=255)

        ccGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGgc                                                                                                                                              <  2:817304/150‑1 (MQ=255)

         cGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGgcg                                                                                                                                             <  2:803388/150‑1 (MQ=255)

         cGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGgcg                                                                                                                                             <  1:1556039/150‑1 (MQ=255)

               gCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAAt                                                                                                                                       <  1:1334172/150‑1 (MQ=255)

                cATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAAtt                                                                                                                                      >  1:583634/1‑150 (MQ=255)

                cATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAAtt                                                                                                                                      >  1:854702/1‑150 (MQ=255)

                  tCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAAtttt                                                                                                                                    <  1:1717943/150‑1 (MQ=255)

                     ggCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGcc                                                                                                                                 >  1:1300787/1‑150 (MQ=255)

                     ggCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGcc                                                                                                                                 <  2:400890/150‑1 (MQ=255)

                       cAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCgg                                                                                                                               >  2:878266/1‑150 (MQ=255)

                              gggCTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGGCGAGGTATCGCCCAGGGTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGt                                                                                                                        <  1:356160/150‑1 (MQ=255)

                                 ttACGATGGCATCGCGATAACCTGCTTCTCTTCAAGCAGCTTCTCGACTACTCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGtttt                                                                                                                     <  2:610210/150‑1 (MQ=255)

                                  tACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCCGTATTCTgcgc                                                                                                                                            >  1:101818/1‑126 (MQ=255)

                                  tACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCCGTATTCTgcgc                                                                                                                                            <  2:101818/126‑1 (MQ=255)

                                  tACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCCGTATTCTgcgc                                                                                                                                            >  1:101819/1‑126 (MQ=255)

                                   aCGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAg                                                                                                                   <  1:878266/150‑1 (MQ=255)

                                    cGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAgg                                                                                                                  <  1:776944/150‑1 (MQ=255)

                                     gATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCTAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGc                                                                                                                 <  2:1417798/150‑1 (MQ=255)

                                        ggCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAgg                                                                                                              >  2:1647714/1‑150 (MQ=255)

                                        ggCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAgg                                                                                                              <  1:1289960/150‑1 (MQ=255)

                                         gCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAggg                                                                                                             <  2:388698/150‑1 (MQ=255)

                                          cATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGa                                                                                                            >  1:1430249/1‑150 (MQ=255)

                                          cATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGa                                                                                                            >  1:1182748/1‑150 (MQ=255)

                                           aTCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAg                                                                                                           >  1:967599/1‑150 (MQ=255)

                                                  tAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGAGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTcc                                                                                                    >  1:1442032/1‑150 (MQ=255)

                                                    gCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAg                                                                                                  >  1:700116/1‑150 (MQ=255)

                                                    gCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAg                                                                                                  >  1:1223886/1‑150 (MQ=255)

                                                      cTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTg                                                                                                <  2:942765/150‑1 (MQ=255)

                                                       tGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCggg                                                                                                                         >  1:603780/1‑124 (MQ=255)

                                                       tGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTgc                                                                                               >  1:480398/1‑150 (MQ=255)

                                                       ggCTTCTCTTCAAGATGCTTCTCGTCTAGGCCAGGATCGGCAAGCTTCGAGGTATCGCCCAGGTTGCTGGTATCGCACGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTgc                                                                                               <  2:283587/149‑1 (MQ=255)

                                                         cTTCTCTTCCAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTgcag                                                                                             <  2:1523541/150‑1 (MQ=255)

                                                         cTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTgcag                                                                                             <  2:216195/150‑1 (MQ=255)

                                                          ttCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTgcagc                                                                                            <  2:949859/150‑1 (MQ=255)

                                                          ttCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTgcagc                                                                                            >  1:1681921/1‑150 (MQ=255)

                                                             tctTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCTGGTTTTAGGCAGGGAGTCGGTCAAGTGCAGCACg                                                                                         >  2:1530577/1‑150 (MQ=255)

                                                             gctTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGTGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACg                                                                                         <  1:1687617/149‑1 (MQ=255)

                                                              ctTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGt                                                                                        >  1:1610181/1‑150 (MQ=255)

                                                               ttCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTc                                                                                       <  1:1115628/150‑1 (MQ=255)

                                                                   agcagcTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGc                                                                                   >  2:695256/1‑150 (MQ=255)

                                                                    gcagcTTCTCGCCTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCg                                                                                  <  1:1271809/150‑1 (MQ=255)

                                                                     cagcTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGt                                                                                 >  2:1626535/1‑150 (MQ=255)

                                                                      agcTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTc                                                                                >  2:1701745/1‑150 (MQ=255)

                                                                           ctcGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAg                                                                           <  2:1539072/150‑1 (MQ=255)

                                                                              gACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCgg                                                                        <  2:1249543/150‑1 (MQ=255)

                                                                               aCTACGACAGGATCGGAAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTAGCCAGCggg                                                                       >  1:187210/1‑150 (MQ=255)

                                                                                   cGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCaa                                                                   >  1:1432361/1‑150 (MQ=255)

                                                                                   cGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCaa                                                                   >  2:588309/1‑150 (MQ=255)

                                                                                     ccAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCATGTCTGGCGTCGCCAGCGGGCCAAtc                                                                 <  2:700116/150‑1 (MQ=255)

                                                                                     ccAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAAtc                                                                 <  2:19171/150‑1 (MQ=255)

                                                                                       aGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAAtctc                                                               <  1:1634148/150‑1 (MQ=255)

                                                                                         gATCGGCAAGCGTCGCGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCtt                                                             <  1:692634/150‑1 (MQ=255)

                                                                                           tCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTa                                                           <  2:636489/150‑1 (MQ=255)

                                                                                           tCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTa                                                           <  1:872223/150‑1 (MQ=255)

                                                                                            cGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTAc                                                          <  2:1223886/150‑1 (MQ=255)

                                                                                              gCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGc                                                        >  2:1418466/1‑150 (MQ=255)

                                                                                              gCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGc                                                        >  1:1339434/1‑150 (MQ=255)

                                                                                               cAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGcc                                                       >  1:87063/1‑149 (MQ=255)

                                                                                                 cgCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCAcc                                                     <  2:1309314/149‑1 (MQ=255)

                                                                                                 aGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCAcc                                                     >  1:1476075/1‑150 (MQ=255)

                                                                                                 aGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGCTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGACCAATCTCTTTACGCAcc                                                     >  2:1055181/1‑150 (MQ=255)

                                                                                                 aGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCAGCAATTATGCGCAGTATACGGCGCATTATTTTGCCGGAGCGCGTTGTAGGCAGGGAGTCGGTGCAGTGCAGCACGACTCGCGTCGCCAGCGGGCCAAGCACTTTACGccac                                                     >  1:699570/1‑147 (MQ=255)

                                                                                                      cGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGtt                                                <  2:583634/150‑1 (MQ=255)

                                                                                                      agAGGTATCGCCCAGGTTGCTTGTATCGCCCGCCGCAATTTTGAGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGtt                                                <  2:854702/149‑1 (MQ=255)

                                                                                                       gAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTg                                               >  1:1665269/1‑150 (MQ=255)

                                                                                                       gAGGTATCGCCCAGGTTGCTGGTATCGACCGCCGCAATTATGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCCGCACGTCAGGCGTCGCCCGCGCGCCAATCACTTTACGCACCCAGTTg                                               >  1:116661/1‑150 (MQ=255)

                                                                                                        aGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGc                                              >  2:542375/1‑150 (MQ=255)

                                                                                                           tATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGACTACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGa                                           <  1:507667/150‑1 (MQ=255)

                                                                                                           tATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCTGTGCAGCACGTCTGGCGTCGCCCGCGGTCCAATCTCTTTACGCGCACATTTGCggc                                           >  2:301996/1‑149 (MQ=255)

                                                                                                           tATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGa                                           >  1:1708398/1‑150 (MQ=255)

                                                                                                           tATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGa                                           >  2:1118837/1‑150 (MQ=255)

                                                                                                           tATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATATTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGa                                           >  1:859445/1‑150 (MQ=255)

                                                                                                              cGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTAGCCTGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACtt                                        >  2:328662/1‑150 (MQ=255)

                                                                                                               gCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTc                                       >  2:349999/1‑150 (MQ=255)

                                                                                                               gCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTc                                       <  1:1626535/150‑1 (MQ=255)

                                                                                                                 ccAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGACCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTg                                     >  1:173634/1‑150 (MQ=255)

                                                                                                                  cAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACACAGTTGCGGACTTCTGc                                    >  2:499976/1‑150 (MQ=255)

                                                                                                                    ggTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTgcagc                                                                                            >  2:663141/1‑92 (MQ=255)

                                                                                                                         cTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGt                             <  2:1585751/150‑1 (MQ=255)

                                                                                                                         cTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGt                             <  2:1681921/150‑1 (MQ=255)

                                                                                                                         cTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGGGTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGt                             >  1:995566/1‑150 (MQ=255)

                                                                                                                             tATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGTTCTg                         >  1:1381279/1‑150 (MQ=255)

                                                                                                                              aTCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGTTCTgg                        <  2:187210/150‑1 (MQ=255)

                                                                                                                              aTCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTAGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGGCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCATTTGCGGACTTCTGCGTACAGTTCTgg                        <  1:1564210/150‑1 (MQ=255)

                                                                                                                                 gcccgccGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGTTCTGGTGa                     >  1:1576957/1‑150 (MQ=255)

                                                                                                                                    cggcaaaaTTATGCGCCGTATTCTGCGCAAAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGTTCTGGTGACgg                  >  2:476819/7‑150 (MQ=255)

                                                                                                                                        gcAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGTTCTGGTGACGGTtcc              <  2:20816/150‑1 (MQ=255)

                                                                                                                                        gcAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCAGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGTTCTGGTGACGGTtcc              <  2:1257707/150‑1 (MQ=255)

                                                                                                                                         cAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGTTCTGGTGACGGTtcct             >  2:622774/1‑150 (MQ=255)

                                                                                                                                          aaTTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGTTCTGGTGACGGTtcctc            <  1:816985/150‑1 (MQ=255)

                                                                                                                                            ttttGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTccc          <  2:1576957/150‑1 (MQ=255)

                                                                                                                                            ttttGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGGTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTccc          <  1:1032876/150‑1 (MQ=255)

                                                                                                                                              atGCGCCGTATTCTGCGCAAAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCg        >  1:524415/2‑150 (MQ=255)

                                                                                                                                                    aGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCGTGATTa  >  1:591026/1‑150 (MQ=255)

AACCGTTACCGACTCGCATCGGGCAATTGTGGGTTACGATGGCATCGCGATAGCCTGCTTCTCTTCAAGCAGCTTCTCGACTACGCCAGGATCGGCAAGCGTCGAGGTATCGCCCAGGTTGCTGGTATCGCCCGCCGCAATTTTGCGCAGAATACGGCGCATAATTTTGCCGGAGCGGGTTTTAGGCAGGGAGTCGGTCCAGTGCAGCACGTCTGGCGTCGCCAGCGGGCCAATCTCTTTACGCACCCAGTTGCGGACTTCTGCGTACAGTTCTGGTGACGGTTCCTCCCCGTGATTA  >  NC_000913/4285380‑4285677

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