BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,986,870	C→T	5.6%	M76I (ATG→ATA) ‡	ycaC‑II ←	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,986,870	0	C	T	5.6%	1064.9 / 30.8	320	M76I (ATG→ATA) ‡	ycaC‑II	putative hydrolase
Reads supporting (aligned to +/- strand): ref base C (136/165); new base T (9/9); total (145/175)
Fisher's exact test for biased strand distribution p-value = 8.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 324 total aligned reads displayed.

GACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGA  >  NC_002947/5986721‑5987019

gACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAAc                                                                                                                                                       >  1:7120367/1‑150 (MQ=255)

 aCAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACa                                                                                                                                                      <  1:8573291/150‑1 (MQ=255)

      gTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATctct                                                                                                                                                 >  1:5505722/1‑150 (MQ=255)

       tCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCtt                                                                                                                                                >  1:6440945/1‑150 (MQ=255)

       tCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCtt                                                                                                                                                >  1:6126246/1‑150 (MQ=255)

          ggAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCa                                                                                                                                             <  1:1260773/150‑1 (MQ=255)

          ggAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCGTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCa                                                                                                                                             >  2:330098/1‑150 (MQ=255)

               gCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCg                                                                                                                                        <  1:5474905/150‑1 (MQ=255)

                  aCGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGc                                                                                                                                     >  2:1841543/1‑150 (MQ=255)

                  aCGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGc                                                                                                                                     >  2:6959611/1‑150 (MQ=255)

                    gCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCAc                                                                                                                                   >  2:7020597/1‑150 (MQ=255)

                      acacaTCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGGTGTCCCAGGCGTGGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTGTCAGCTCGGGCACCa                                                                                                                                 >  2:3952408/1‑150 (MQ=255)

                         caTCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCg                                                                                                                              >  2:5878592/1‑150 (MQ=255)

                          aTCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCgg                                                                                                                             <  2:2867869/150‑1 (MQ=255)

                          aTCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCgg                                                                                                                             <  2:7486812/150‑1 (MQ=255)

                             ggTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGcc                                                                                                                          <  2:2102727/150‑1 (MQ=255)

                              gTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCa                                                                                                                         >  1:4472995/1‑150 (MQ=255)

                                   cacaCCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCgg                                                                                                                    >  1:1597012/1‑150 (MQ=255)

                                      acCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGcc                                                                                                                 <  1:386457/150‑1 (MQ=255)

                                        cGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCtt                                                                                                               <  2:3077491/150‑1 (MQ=255)

                                            gatgatCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTc                                                                                                           >  2:7694894/1‑150 (MQ=255)

                                            gatgatCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTc                                                                                                           >  1:1933730/1‑150 (MQ=255)

                                            gatgatCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTc                                                                                                           >  2:4436954/1‑150 (MQ=255)

                                              tgatCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGa                                                                                                         >  1:4338227/1‑150 (MQ=255)

                                               gatCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGaa                                                                                                        <  1:4167122/150‑1 (MQ=255)

                                               gatCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGaa                                                                                                        <  2:2378896/150‑1 (MQ=255)

                                               gatCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGaa                                                                                                        <  1:2965917/150‑1 (MQ=255)

                                                atCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAg                                                                                                       >  1:2766618/1‑150 (MQ=255)

                                                 tCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGc                                                                                                      >  1:6849007/1‑150 (MQ=255)

                                                  cAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCt                                                                                                     <  1:4566758/150‑1 (MQ=255)

                                                   aGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtg                                                                                                    <  1:9010169/150‑1 (MQ=255)

                                                     ttgcttgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtggt                                                                                                  >  2:783113/1‑150 (MQ=255)

                                                      tgcttgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtggtg                                                                                                 <  1:2343109/150‑1 (MQ=255)

                                                       gcttgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtggtgg                                                                                                >  2:4817021/1‑150 (MQ=255)

                                                         ttcttgCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGCCGTTTTCAAAACTggtggtc                                                                                              >  1:8895561/6‑149 (MQ=2)

                                                          tgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGa                                                                                             <  1:3293593/150‑1 (MQ=255)

                                                          tgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGa                                                                                             <  2:3708978/150‑1 (MQ=255)

                                                           gcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAg                                                                                            >  2:8487287/1‑150 (MQ=255)

                                                             ggCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGa                                                                                          >  2:6054696/1‑150 (MQ=255)

                                                             ggCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGa                                                                                          <  1:1012923/150‑1 (MQ=255)

                                                                 cGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGt                                                                                      >  1:1718980/1‑150 (MQ=255)

                                                                  gTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTc                                                                                     <  2:4940493/150‑1 (MQ=255)

                                                                  gTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGCCGTTTTCTAAACTGGTGGTCAgaacagtc                                                                                     <  2:1030459/150‑8 (MQ=2)

                                                                    ggCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCgg                                                                                   <  2:6473443/150‑1 (MQ=255)

                                                                     gCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGc                                                                                  <  2:1993407/150‑1 (MQ=255)

                                                                     gCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGCCGTTTTCAAAACTGGTGGTCAGAACAGTCGGc                                                                                  >  1:1345669/1‑150 (MQ=2)

                                                                      ccTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCa                                                                                 <  2:3185543/150‑1 (MQ=255)

                                                                       cttcactGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGCCGTTTTCAAAACTGGTGGTCAGAACAGTCGGCAg                                                                                >  1:7912117/7‑150 (MQ=2)

                                                                          gATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGcc                                                                             <  1:300449/150‑1 (MQ=255)

                                                                             tGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCaaa                                                                          <  2:5505722/150‑1 (MQ=255)

                                                                               ccTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGa                                                                        <  2:4662381/150‑1 (MQ=255)

                                                                                cTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGaa                                                                       <  2:1315379/150‑1 (MQ=255)

                                                                                cTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGaa                                                                       >  1:9052057/1‑150 (MQ=255)

                                                                                 ttGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAAc                                                                      >  2:8762389/1‑150 (MQ=255)

                                                                                 ttGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAAc                                                                      <  2:6461741/150‑1 (MQ=255)

                                                                                   gACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACtt                                                                    <  2:1528957/150‑1 (MQ=255)

                                                                                   gACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACtt                                                                    <  1:7500451/150‑1 (MQ=255)

                                                                                    aCGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTg                                                                   <  2:6146069/150‑1 (MQ=255)

                                                                                    aCGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTg                                                                   <  2:5153153/150‑1 (MQ=255)

                                                                                    aCGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTg                                                                   <  2:5790608/150‑1 (MQ=255)

                                                                                    aCGAAGTCTTCGTTGTCCCAGGCGTTGATCTGACCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTg                                                                   <  1:662817/150‑1 (MQ=255)

                                                                                      gAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGc                                                                 >  2:1342329/1‑150 (MQ=255)

                                                                                      gAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGc                                                                 >  2:7590254/1‑150 (MQ=255)

                                                                                        aGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGCGAGGATGGTCGGCAGGCCAAAGAACTTGGCCa                                                               <  2:7235462/150‑1 (MQ=255)

                                                                                         gTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAg                                                              <  2:424936/150‑1 (MQ=255)

                                                                                          tCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAgg                                                             >  1:8747170/1‑150 (MQ=255)

                                                                                            ttCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTc                                                           >  2:5982808/1‑150 (MQ=255)

                                                                                                 tGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCa                                                      >  2:5100602/1‑150 (MQ=255)

                                                                                                     ccAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCCCCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGc                                                  >  1:1978544/1‑150 (MQ=255)

                                                                                                         gCGTTGATCTGGCCGGGACGGGCAATTTACGGCGCATCCGGGAACATCTCTTTCAGCTCTGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATTGTCGGCTGGCCAACGAACTTGGCCAGGTCGACCAAGGCCAGc                                              <  1:363483/150‑1 (MQ=18)

                                                                                                              gATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACgtt                                         <  1:1413195/150‑1 (MQ=255)

                                                                                                              gATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACgtt                                         <  1:8774647/150‑1 (MQ=255)

                                                                                                                   ggCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCt                                    <  1:4712669/150‑1 (MQ=255)

                                                                                                                    gCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCtt                                   <  2:3567784/150‑1 (MQ=255)

                                                                                                                    gCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCtt                                   <  1:5056447/150‑1 (MQ=255)

                                                                                                                            gggCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGt                           >  1:5191274/1‑150 (MQ=255)

                                                                                                                             ggCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGtt                          <  1:4404796/150‑1 (MQ=255)

                                                                                                                             ggCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGtt                          <  2:9094290/150‑1 (MQ=255)

                                                                                                                             ggCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGtt                          <  1:8190792/150‑1 (MQ=255)

                                                                                                                             ggCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGtt                          <  1:6348515/150‑1 (MQ=255)

                                                                                                                             ggCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGtt                          <  1:7020597/150‑1 (MQ=255)

                                                                                                                             ggCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGtt                          >  1:6966990/1‑150 (MQ=255)

                                                                                                                              gCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTg                         <  2:1687334/150‑1 (MQ=255)

                                                                                                                              gCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTg                         <  2:4309977/150‑1 (MQ=255)

                                                                                                                              gCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTg                         <  1:330098/150‑1 (MQ=255)

                                                                                                                              gCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTg                         >  2:2594795/1‑150 (MQ=255)

                                                                                                                                 aTGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAATAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGc                      <  1:7186564/150‑1 (MQ=255)

                                                                                                                                 aTGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGc                      <  1:5313367/150‑1 (MQ=255)

                                                                                                                                       ggCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGaaaaa                >  1:1489691/1‑150 (MQ=255)

                                                                                                                                         cgcATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATc              >  2:5356778/1‑150 (MQ=255)

                                                                                                                                         cgcATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATc              <  1:5878592/150‑1 (MQ=255)

                                                                                                                                            aTCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTg           <  2:2766618/150‑1 (MQ=255)

                                                                                                                                            aTCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTg           <  1:7694894/150‑1 (MQ=255)

                                                                                                                                              ccGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCa         <  2:8136415/150‑1 (MQ=255)

                                                                                                                                                gggAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCAcc       >  2:6708301/1‑150 (MQ=255)

                                                                                                                                                gggAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCAcc       >  1:4847933/1‑150 (MQ=255)

                                                                                                                                                gggAACATCTCTTTCAGCACGGGCACCAGCGGGCCATTAGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCAcc       >  1:8348955/1‑150 (MQ=255)

                                                                                                                                                 ggAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCa      >  2:6383918/1‑150 (MQ=255)

                                                                                                                                                 ggAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCa      >  2:5639370/1‑150 (MQ=255)

                                                                                                                                                     cATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCga  <  2:4112214/150‑1 (MQ=255)

GACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAGGTCGCCCAAGGCCAGCACGTTGTTCTTGAATTCGTTGGGCGAAAAATCCTGCACCAGCGA  >  NC_002947/5986721‑5987019

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