BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	160,517	T→G	5.4%	A63A (GCA→GCC) ‡	PP_0150 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	160,517	0	T	G	5.4%	297.5 / 6.5	93	A63A (GCA→GCC) ‡	PP_0150	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (41/46); new base G (3/2); total (45/48)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

CAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTCGGCCAG  >  NC_002947/160369‑160663

agatctcgcacgaGGTTCTCAGTGCGAGTATGGCTTCAGAGAATACCTGGTCGCTGGCGGCGGGCGCTTCCGGTACTCGGTCGCGCTGGGTTTTGAAGGCCGCCTCGTGTTGTACGCGGCGGATGTCACCGTTGCTGGCGACGAAGCtt                                                                                                                                                     <  4:148294/143‑1 (MQ=255)

 aGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTg                                                                                                                                                    >  2:70565/1‑149 (MQ=255)

    cacgcacgCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAg                                                                                                                                                 <  5:241337/149‑1 (MQ=255)

     acgcacgcGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGc                                                                                                                                                >  3:358120/1‑149 (MQ=255)

     acgcacgcGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGc                                                                                                                                                >  5:162011/1‑149 (MQ=255)

       gcacgcGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCtt                                                                                                                                                     >  3:148294/1‑142 (MQ=255)

       gcacgcGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTAGTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGcgt                                                                                                                                              >  4:234970/1‑149 (MQ=255)

             gcTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGACGCTGCAAGcttcgtcgc                                                                                                                                        <  3:342/149‑2 (MQ=255)

             gcTTCACAGTGCGAGGATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGcgtcgtcgt                                                                                                                                        <  7:58490/149‑1 (MQ=255)

              cTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTg                                                                                                                                       <  1:186051/149‑1 (MQ=255)

               ttCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTgg                                                                                                                                      >  4:116552/1‑149 (MQ=255)

               ttCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTgg                                                                                                                                      >  4:84047/1‑149 (MQ=255)

               ttCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTgg                                                                                                                                      >  6:237234/1‑149 (MQ=255)

                  acaGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGACGCTGCAAGcttcgtcgc                                                                                                                                        <  7:268111/144‑2 (MQ=255)

                  acaGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGACGCTGCAAGcttcgtcgc                                                                                                                                        >  8:268111/1‑143 (MQ=255)

                  acaGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGAAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGATGGCGACGAAGCTTGCAGCGTCGTCGTGGGCt                                                                                                                                   >  3:248018/1‑149 (MQ=255)

                      tGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTcgc                                                                                                                               >  2:104796/1‑149 (MQ=255)

                         gAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGca                                                                                                                            <  6:322565/149‑1 (MQ=255)

                         gAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGca                                                                                                                            <  3:116552/149‑1 (MQ=255)

                               ggCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTAGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTAGCAGCGTCGTCGTGGGCTTCGCGCACCACtt                                                                                                                      >  5:197208/1‑149 (MQ=255)

                                gCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACtt                                                                                                                      <  4:319038/148‑1 (MQ=255)

                                 cTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTg                                                                                                                    <  3:169354/149‑1 (MQ=255)

                                 cTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTg                                                                                                                    <  6:52450/149‑1 (MQ=255)

                                 cTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTg                                                                                                                    <  3:234970/149‑1 (MQ=255)

                                 cTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTg                                                                                                                    >  7:39573/1‑149 (MQ=255)

                                 cTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTg                                                                                                                    >  8:221580/1‑149 (MQ=255)

                                 cTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTg                                                                                                                    <  1:104796/149‑1 (MQ=255)

                                  ttCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTg                                                                                                                    >  8:127482/1‑148 (MQ=255)

                                    cAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGa                                                                                                                  <  1:279629/148‑1 (MQ=255)

                                       ccAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcg                                                                                                              >  8:241965/1‑149 (MQ=255)

                                        cAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGcgtcgtcgt                                                                                                                                        <  2:305615/122‑1 (MQ=255)

                                        cAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGcgtcgtcgt                                                                                                                                        >  1:305615/1‑122 (MQ=255)

                                        cAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcgc                                                                                                             <  6:162011/149‑1 (MQ=255)

                                              cTGGTCGCTGGCGCCGCGCGCTTCCGGTTCTCGCTCGCGCAGGGTATTGAAGGACGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATc                                                                                                        <  7:241965/148‑1 (MQ=255)

                                                  tCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggt                                                                                                   >  6:40087/1‑149 (MQ=255)

                                                     cTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggtggt                                                                                                <  3:84047/149‑1 (MQ=255)

                                                      tGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggtggtg                                                                                               <  3:318585/149‑1 (MQ=255)

                                                      tGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggtggtg                                                                                               >  1:181976/1‑149 (MQ=255)

                                                              cgcgcTTCCGGTACTCGCTGGCGCAGGGTTTTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGTGGAGGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTTGTG‑GTGTCTGa                                                                                      <  6:197208/149‑1 (MQ=255)

                                                               gcgcTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAg                                                                                     <  5:237234/149‑1 (MQ=255)

                                                                cgcTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggtggtg                                                                                               <  4:246674/139‑1 (MQ=255)

                                                                cgcTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggtggtg                                                                                               >  3:246674/1‑139 (MQ=255)

                                                                cgcTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAgg                                                                                    >  3:368091/1‑149 (MQ=255)

                                                                   ttCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTg                                                                                  >  4:145924/1‑148 (MQ=255)

                                                                     ccGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGa                                                                                 >  5:244798/1‑147 (MQ=255)

                                                                     ccGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGa                                                                                 <  6:244798/147‑1 (MQ=255)

                                                                     caggTACTCGCTCGCGCTGGGTAGTGAAGGCCGCCTCGATTTGTACGCCGCGGATGTCACAGTTGCTGGAGACGAATCTTGCATCTTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCTCGGATCGAGGTTGTG‑GTGTCTGAgtttaa                                                                                <  4:248018/146‑6 (MQ=255)

                                                                      cGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCAGCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGt                                                                              <  4:358120/149‑1 (MQ=255)

                                                                      cGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGt                                                                              >  7:44/1‑149 (MQ=255)

                                                                      cGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGt                                                                              >  2:194122/1‑149 (MQ=255)

                                                                        gTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGtc                                                                             >  7:242863/1‑148 (MQ=255)

                                                                           ctcgctcgCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGtcgat                                                                          <  1:194122/148‑1 (MQ=255)

                                                                              gctcgcGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcgc                                                                                                             <  1:150701/111‑1 (MQ=255)

                                                                              gctcgcGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTcgc                                                                                                             >  2:150701/1‑111 (MQ=255)

                                                                              gctcgcGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGtcgatcg                                                                        <  4:368091/147‑1 (MQ=255)

                                                                                  gcgcAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGc                                                                  >  2:281246/1‑149 (MQ=255)

                                                                                    gcAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑gtgt                                                                                          <  8:296477/123‑1 (MQ=255)

                                                                                    gcAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑gtgt                                                                                          >  7:296477/1‑123 (MQ=255)

                                                                                       gggTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAg                                                             <  8:114615/149‑1 (MQ=255)

                                                                                       gggTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAg                                                             <  5:40087/149‑1 (MQ=255)

                                                                                             tGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGgcacgc                                                       >  8:117283/1‑149 (MQ=255)

                                                                                               aaGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGgcacgcac                                                     >  2:87473/1‑149 (MQ=255)

                                                                                               aaGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GGGTCTGAGGTGAAGTCGATCGAACGGCCAAAGgcacgcac                                                     >  1:124797/1‑149 (MQ=255)

                                                                                               aaGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTC‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGgcacgcac                                                     >  5:317731/1‑149 (MQ=255)

                                                                                                    cgccTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACgatga                                                <  5:142747/149‑1 (MQ=255)

                                                                                                     gccTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACgatgat                                               <  1:87473/149‑1 (MQ=255)

                                                                                                      ccTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGatgatg                                              <  6:317731/149‑1 (MQ=255)

                                                                                                       cTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGt                                             >  1:44612/1‑149 (MQ=255)

                                                                                                         cGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGt                                                                              <  8:205942/114‑1 (MQ=255)

                                                                                                          gAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGaa                                                                                <  4:218416/111‑1 (MQ=255)

                                                                                                          gAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGaa                                                                                >  3:218416/1‑111 (MQ=255)

                                                                                                          gAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGtt                                            <  3:152506/147‑1 (MQ=255)

                                                                                                          gAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGc                                          <  8:39573/149‑1 (MQ=255)

                                                                                                            gTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCtg                                        <  7:221580/149‑1 (MQ=255)

                                                                                                             ttGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTGCGTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCtg                                        <  1:281246/149‑1 (MQ=255)

                                                                                                             ttGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCtg                                        <  2:181976/148‑1 (MQ=255)

                                                                                                               gTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCtggt                                      <  7:127482/148‑1 (MQ=255)

                                                                                                                tACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCtggtg                                     <  5:228457/148‑1 (MQ=255)

                                                                                                                  cgccgcGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGc                                   >  7:259589/1‑148 (MQ=255)

                                                                                                                    ccgcGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGccac                                >  3:32078/1‑149 (MQ=255)

                                                                                                                    ccgcGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGccac                                >  1:254082/1‑149 (MQ=255)

                                                                                                                               aCCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGTGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCtgtg                     >  1:202521/1‑149 (MQ=255)

                                                                                                                                ccGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGc                    <  4:32078/149‑1 (MQ=255)

                                                                                                                                 cGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggtggtg‑gt                                                                                            <  7:58951/76‑1 (MQ=255)

                                                                                                                                 cGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAggtggtg‑gt                                                                                            >  8:58951/1‑76 (MQ=255)

                                                                                                                                     gCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTCGCTGGTGGCCACCAGGGTCTGTGCCTgg                <  3:145924/148‑1 (MQ=255)

                                                                                                                                       ccacgacgaCGCTGCAAGCTTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGc              >  4:121276/4‑148 (MQ=255)

                                                                                                                                         gCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAg           <  8:242863/149‑1 (MQ=255)

                                                                                                                                         acgacgaCGCTGCAAGCTTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAg           >  8:289444/2‑149 (MQ=255)

                                                                                                                                           gacgaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCtggtg                                     >  3:188777/1‑121 (MQ=255)

                                                                                                                                           gacgaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCtggtg                                     <  4:188777/121‑1 (MQ=255)

                                                                                                                                           gacgaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGt         >  2:195807/1‑149 (MQ=255)

                                                                                                                                           gacgaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGGCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGt         >  2:159414/1‑149 (MQ=255)

                                                                                                                                            acgaAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTc        <  8:44/149‑1 (MQ=255)

                                                                                                                                                aGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTCGGc     >  7:263459/1‑148 (MQ=255)

                                                                                                                                                   ttGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTCGGCcag  >  4:132173/1‑148 (MQ=255)

CAGGCACGCACGCGCTTCACAGTGCGAGTATGGCTTCAGCCAATACCTGGTCGCTGGCGCCGCGCGCTTCCGGTACTCGCTCGCGCAGGGTATTGAAGGCCGCCTCGAGTTGTACGCCGCGGATGTCACCGTTGCTGGCGACGAAGCTTGCAGCGTCGTCGTGGGCTTCGCGCACCACTTTGGAGTCGCGGATCGAGGTGGTG‑GTGTCTGAGGTGAAGTCGATCGAACGGCCAAAGGCACGCACGATGATGTTGCTGGTGGCCACCAGGGTCTGTGCCTGGGCCAGGTCGGCCAG  >  NC_002947/160369‑160663

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: `‑`