BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,102,247	G→C	5.9%	intergenic (+241/+51)	PP_mr66 → / ← acsA‑I	P6/acetyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,102,247	0	G	C	5.9%	310.2 / 5.4	107	intergenic (+241/+51)	PP_mr66/acsA‑I	P6/acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand): ref base G (65/34); new base C (3/4); total (69/38)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

ATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGCGTCGTACTC  >  NC_002947/5102099‑5102396

aTGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGa                                                                                                                                                      <  2:4579068/150‑1 (MQ=255)

aTGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCG‑TAGCCGGAACTGa                                                                                                                                                      >  1:2330917/1‑150 (MQ=255)

aTGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCGTTCCCGCGAAATCGGGAGGGCAGGGGGCCGCG‑TAGCCTGAACTGa                                                                                                                                                      >  1:4865364/1‑150 (MQ=255)

 tGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGAc                                                                                                                                                     <  1:1004385/150‑1 (MQ=255)

 tGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGAc                                                                                                                                                     <  2:2424552/150‑1 (MQ=255)

         aaCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTcgc                                                                                                                                             <  1:734083/150‑1 (MQ=255)

         aaCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCGGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTcgc                                                                                                                                             <  1:552710/150‑1 (MQ=255)

          aCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAGGAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTcgcg                                                                                                                                            <  2:628858/150‑1 (MQ=18)

          aCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTcgcg                                                                                                                                            <  2:4278259/150‑1 (MQ=255)

           cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                                                           <  2:1676035/150‑1 (MQ=255)

           cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                                                           <  2:1676715/150‑1 (MQ=255)

           cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                                                           <  2:4496184/150‑1 (MQ=255)

           cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                                                           <  1:1746098/150‑1 (MQ=255)

           cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                                                           <  1:2280064/150‑1 (MQ=255)

            gggCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGGATGCGGGAGGTCAGGGGACAGCT‑TAGCCGGAACGGACGCAGTCGCcgg                                                                                                                                          >  1:4896608/1‑150 (MQ=255)

            gggCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCggg                                                                                                                                          <  1:644707/150‑1 (MQ=255)

              gCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCa                                                                                                                                        <  1:3139812/150‑1 (MQ=255)

               ccGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCaa                                                                                                                                       <  1:2664845/150‑1 (MQ=255)

                 ggAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGCATGCGGGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGc                                                                                                                                     >  2:3001066/1‑150 (MQ=255)

                 ggAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGc                                                                                                                                     >  2:3853404/1‑150 (MQ=255)

                       gTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGCATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTc                                                                                                                               >  1:5094062/1‑150 (MQ=255)

                             tGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGCATGCGTGAGGTCCGGGGACAGCG‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAg                                                                                                                         >  2:3682296/1‑150 (MQ=255)

                             tGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCA‑TAGCCAGAAAAGACGCAATCGCGGGCAAGCCCGCTCCTGCAg                                                                                                                         >  2:1690456/1‑150 (MQ=255)

                              ccggcCCCGCACAGGTTCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGGTTGCCCGCGAAATCGTTATGTGGTGGGTCAGCT‑TATGCTTAAAATACGCAAAAGCGGGCAAACCCGCTCCTGCaag                                                                                                                        >  1:1964990/4‑148 (MQ=255)

                               gcgcCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGGCCGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCGGCAggg                                                                                                                       >  2:3167210/1‑150 (MQ=255)

                               gcgcCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGGCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACCGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCACCTGCAggg                                                                                                                       >  2:3167204/1‑150 (MQ=255)

                                    ccGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCGTGCCCCCGAAAGCGTGTGGTCCGGGGACAGCT‑TAGCCTGGACAGACGCAATAGCCGGCAAACCCCCTCCAGCAGGGGGGAt                                                                                                                  >  1:3909874/1‑150 (MQ=255)

                                      gAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAAAGCGTGAGGTCCGGGGGCAGCG‑TAGCCGGAACTGACGCATTCGCGGGCAAGCCCGCACCAGCAGGGGGTATTg                                                                                                                >  2:490835/1‑150 (MQ=255)

                                           tgtgCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGGCAGGGGACAGCT‑TAGCCGGAACGGACGCATTCGCGGGCCAGCCCGCTCCTGCAGGGGGGATGGCCggg                                                                                                           >  2:5367857/1‑150 (MQ=255)

                                             tgCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGaa                                                                                                         >  2:3747845/1‑150 (MQ=255)

                                             tgCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCAGCAGGGGGTATTGCCGGGaa                                                                                                         >  2:2485257/1‑150 (MQ=255)

                                             ggCGGGGCCGGTGACAGAACCGCACCCGCGGTTTGT‑GGGTAACCCCGCCAACCGGGGGGTGGCCCGCGACGGGGGGCGGTCCGGGGACCGCT‑TCGCCTCACCTGCCGCTTTCGGGGGCACGCCCGCTCCTGCAGGGGGTATTGCCGGGaa                                                                                                         <  2:4982169/149‑1 (MQ=255)

                                                ggggCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAATGCGGGCTTGCCCGCGAATGCGTGTGGGCAGGGGACCGCT‑TAGCCGGAACTGACGCATTCGCGGGCCAGCCCGCGCCGGCCGGGGGTGTTGCGGGGGAggg                                                                                                      >  2:4729954/1‑150 (MQ=255)

                                                ggggCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGGTTGCCCGCGAAAGCGGGAGGTGAGGGGACAGCG‑TATCCTGGAATGACGCATTAGCGGGCAAGCCCGCTCCGGCAGGGGGTATTTGCGGGAAggg                                                                                                      >  1:730128/1‑150 (MQ=255)

                                                      ggTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGGGGTCAGGGGGCAGCT‑TAGCCGGTACTGACGCAATAGCGGGCCAGCCCGCTCCGGCAGGGGGTGTTGCCGGGAAGGGGTTaag                                                                                                >  1:4024954/1‑148 (MQ=255)

                                                          aCAGAACATTCCCCGATTTTTGTCGGTTAACAATG‑AAAACGGGGGGTTGCCCGCGATTGCGTGTGGACTGGGGACAGCT‑TAGCCTGACGTCCCGCTTTCGCGGCCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGa                                                                                            <  2:4919072/150‑1 (MQ=255)

                                                           ccgAACATTAACCTATGTTTGTCGGGCAACAAGC‑CAAACGCGGGCTGGCCCGCGAAGGGGTGAGTTCAGGGGACCGCT‑TAGCCTGAACTGACGCATTCGCGGGCACGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGaa                                                                                           <  2:2729917/148‑1 (MQ=255)

                                                            aGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAAAGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGGATTGCCGGGAAGGGTTTACGCCGAAg                                                                                          >  2:3970917/1‑150 (MQ=255)

                                                            aGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCGTGCCCGCGAATGCGTGAGGTCAGGGGCCACGG‑TAGCCTGACCGGACGCATTCGCGGGCCAACCCGCGCCTGCAGGGGGTTTTTCCGGGGCGGGTTTACGCCGAAg                                                                                          >  1:2937210/1‑150 (MQ=255)

                                                               aCAGTAACCGATGATTGTCGGGTAACAAGG‑CAAACGCGGGCTTGCCCGCGATGGCGTGTGGTCAGGGGACAGCT‑TAGCCTGACCTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCa                                                                                       <  2:4637955/150‑1 (MQ=255)

                                                                 aGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCGTGCCCCCGAATGCGGGAGGTCCGGGGACAGCT‑TAGCCGGAACTGACGCAATAGCGGGCAAGCCCCCTCCGGCAGGGGGTATTGCCGGGAAGGGTTTAAGCCCAAGCCAgg                                                                                     >  1:864290/1‑150 (MQ=255)

                                                                   tAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCCCGAATGCGTGAGGTCAGGGGACAGCG‑TAGCCGGAACAGACGCAAACGCCGGCAAGCCCCCGCCAGCAGGGGGTATTGCCGGGAAGGGGTTACGCCGAAGCCAGGtt                                                                                   >  1:4671513/1‑150 (MQ=255)

                                                                    aaCCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACAGACGCATTCGCGGGCAAGCCCGCTCCAGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTc                                                                                  >  2:2842918/1‑150 (MQ=255)

                                                                     cccGATGATGGTCGGTTAACAATG‑CAAACGGGGGGTTGCCCGCGAAGGCGTGAGGTCCGGGGACATCT‑TAGCCTGAACTGCCGCTTTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCa                                                                                 <  2:2220205/149‑1 (MQ=255)

                                                                      ccGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTcgcg                                                                                                                                            >  1:326303/1‑90 (MQ=255)

                                                                      ccGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTcgcg                                                                                                                                            <  2:326303/90‑1 (MQ=255)

                                                                       cGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTTCCCCCGAATTCGTGAGGTCAGGGGACAGCT‑TAGCCTGAAATGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCCGGAAGGGTTTAAGCCGAAGCCCGGTTAATc                                                                               >  2:1430519/1‑150 (MQ=255)

                                                                       cGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGCAAGCGTGAGGGCAGGGGACAGCT‑TAGCCAGAACTGACCCAATCGCGGGCAAGCCCGCACCAGCAGGGGGTATTGCCGGGCAGGGTTTACGCCGAAGCCAGGTTCaaa                                                                               >  2:1430473/1‑148 (MQ=255)

                                                                       cGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCCCGAAAGCGTGATGTCAGGGGCCAGCG‑TAGCCAGAAAGGACCCATACGCGGGCAAGCCCCCACCAGCACGGGGGATTGCCGGGAAGGGGTTACGCCCAACCCCGGTTCATc                                                                               >  2:3017685/1‑150 (MQ=255)

                                                                                tcggtaaaCATTG‑CAAACGCGGGGTTGCCCGCGAATGGGTGAGGTCAGGGGACAGCT‑TAGCCCCAACTGCCGCATTCGCGGGCACGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTgg                                                                      <  2:3131418/144‑1 (MQ=255)

                                                                                  gggTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGCCAGCT‑TAGCCTGAACTGGCGCATTCGCGGGCAAGCCCGCTCCTGCCGGGGGTCTTGGCGGGAAGGGTTTACGCCGAAGCCCGGTTCATCGCCTTGTGGGt                                                                    >  2:1750876/1‑150 (MQ=255)

                                                                                  gggTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCATGTACTGACGCATTCGCGGGCAAGCCCGCTCCTGCCGGGGGTATTTCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGt                                                                    >  2:1267429/1‑150 (MQ=255)

                                                                                   ggTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACCGCT‑TAGCCTGAACGGACGCATTCGCGGGCAAGCCCGCGCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTa                                                                   >  2:4487554/1‑150 (MQ=255)

                                                                                   ggTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TATCCGGAACAGACGCATTCGCGGGCAAGCCCGCTCCGGGAGGGGGTATTGCCGGGAAGGGGTTACGCCGAAGCCAGGTTCATCGCCGTGTGGGTa                                                                   >  2:1640248/1‑150 (MQ=255)

                                                                                   ggTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACACCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCAGCAGGGGGTATTGCCGGGCAGGGGTTACGCCGAAGCCAGGTTCATAGCCTTGTGGGTa                                                                   >  2:4487327/1‑150 (MQ=39)

                                                                                          attccAAACGCGGGCTGGCCGGCGAATGGGTGCGGCCGGGGGACAGAT‑TCGCCTGAAGCCACGCATGCGGGGGCACGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATc                                                             <  2:3859662/146‑1 (MQ=255)

                                                                                              caaccgCGGGCTTGCCCGCTACTGCGTGCGGTCTGGGGACCGCT‑TAGCCTCACCTGCCGCTTCCGCGGGCACGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGt                                                         <  2:579887/145‑1 (MQ=255)

                                                                                              cAAACGGGGGCTGGCCCGCTAATGGTTGTGGTCAGGGGACCGCT‑TAGCCTTACCTGACGCATTCCCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGt                                                         <  2:3457350/150‑1 (MQ=255)

                                                                                                 acgggggCTTGCCCGCGAATGCGTGTGGTCTGGGGACAGCT‑TAGCCCGACCTCCCGCTTTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCt                                                      <  2:1759125/146‑1 (MQ=255)

                                                                                                       gCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACca                                                >  1:3094746/1‑150 (MQ=18)

                                                                                                            cccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACAccc                                           >  2:3404244/1‑150 (MQ=255)

                                                                                                             ccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCg                                          >  1:652510/1‑150 (MQ=255)

                                                                                                             ccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCg                                          >  1:2401842/1‑150 (MQ=255)

                                                                                                             cagagAATGCGTGAGGTCAGGGGACAGCT‑TATCCTGAACTGACTCATTCGCGGGAAATCCCGCTCATGCAGGGGTTATTGCGGGGAAGGGTTTACGCCGAAGCGAGGTTCATCGACTTGTGGGTATCGATCAGGTGCTGTACCACACACg                                          >  2:530035/5‑150 (MQ=255)

                                                                                                              cgcgAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCgg                                         >  2:2464582/1‑150 (MQ=255)

                                                                                                              cgcgAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGACGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCgg                                         >  1:4080327/1‑150 (MQ=21)

                                                                                                              cgcgAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCgg                                         >  1:1414898/1‑150 (MQ=21)

                                                                                                               gcgAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCTTTCGGGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCggg                                        <  2:2240806/150‑1 (MQ=255)

                                                                                                                cgatggggtgtGGTCGGGGCACCGCT‑TAGCCTGAAGTGCCCCTTGCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGt                                       <  2:603929/143‑1 (MQ=255)

                                                                                                                   aTGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAg                                    >  1:2571579/1‑150 (MQ=255)

                                                                                                                       gTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTCTTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCaa                                >  2:1191128/1‑150 (MQ=255)

                                                                                                                       gTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCaa                                >  2:529906/1‑150 (MQ=255)

                                                                                                                       gTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCaa                                >  2:590087/1‑150 (MQ=255)

                                                                                                                       gTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCaa                                >  2:5355450/1‑150 (MQ=255)

                                                                                                                            gTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCg                           >  1:3653682/1‑150 (MQ=255)

                                                                                                                             tctggggACTGCT‑TCGCCTCAATGGCCGCTTTGGCGGGCACGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCga                          <  1:5082286/147‑1 (MQ=255)

                                                                                                                             tctggggACCGCT‑TCGCCTCACCTGCCGCATGGGCGGGCACGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCga                          <  2:2848845/147‑1 (MQ=255)

                                                                                                                             tCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCga                          >  1:1226334/1‑150 (MQ=255)

                                                                                                                               ggggtACAGTT‑TCGCCTTAACTGCGGCATTCGGGGCAAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCgaga                        <  2:3770987/149‑1 (MQ=255)

                                                                                                                               aGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCgaga                        >  1:307061/1‑150 (MQ=255)

                                                                                                                                 gggACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGata                      >  2:5020782/1‑150 (MQ=255)

                                                                                                                                 gggACAGCT‑TAGCCTGAACTGACGCATTCGCGGGAAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGata                      >  1:802641/1‑150 (MQ=255)

                                                                                                                                  ggACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGatat                     >  2:2295901/1‑150 (MQ=255)

                                                                                                                                      aGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGcc                 >  1:206346/1‑150 (MQ=255)

                                                                                                                                        cT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCa               >  2:523784/1‑150 (MQ=255)

                                                                                                                                        cT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCa               >  2:1433527/1‑150 (MQ=255)

                                                                                                                                        cT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCa               >  2:523985/1‑150 (MQ=255)

                                                                                                                                        cT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCa               >  2:523793/1‑150 (MQ=255)

                                                                                                                                         t‑tAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAg              >  1:4781107/1‑150 (MQ=255)

                                                                                                                                           tAGCGTGAAGTGACTATTTCGGGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGATGGGTTTACGCCGAAGCCTGGTTCATCGCCTTGTGGGTATCTATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAgg             <  2:2923296/150‑1 (MQ=255)

                                                                                                                                           tAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAgg             >  2:4697367/1‑150 (MQ=255)

                                                                                                                                             gCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCAAGGGc           >  2:2051399/1‑150 (MQ=255)

                                                                                                                                                tGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTg                                                     <  2:3449898/105‑1 (MQ=255)

                                                                                                                                                tGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTg                                                     >  1:3449898/1‑105 (MQ=255)

                                                                                                                                                tGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcgtc        >  1:2088277/1‑150 (MQ=255)

                                                                                                                                                tGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcgtc        >  1:2088803/1‑150 (MQ=255)

                                                                                                                                                tGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcgtc        >  2:1082247/1‑150 (MQ=255)

                                                                                                                                                tGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcgtc        >  2:1875435/1‑150 (MQ=255)

                                                                                                                                                 gAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcgtcg       >  2:4127225/1‑150 (MQ=255)

                                                                                                                                                      gACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGCGTCGTACTc  >  2:509973/1‑150 (MQ=255)

ATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATG‑CAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGCGTCGTACTC  >  NC_002947/5102099‑5102396

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