BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,102,244	A→C	9.4%	intergenic (+238/+54)	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,244	0	A	C	9.4%	427.3 / 27.1	150	intergenic (+238/+54)	PP_mr66/acsA‑I	P6/acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand): ref base A (93/42); new base C (2/13); total (95/55)
Fisher's exact test for biased strand distribution p-value = 4.25e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

Only 100 of 156 total aligned reads displayed.

TCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGCGTCG  >  NC_002947/5102095‑5102391

tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACCGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTAGCCCGCGCATGCGGGAGGTCAGGGGCCAGCG‑TAGCCTGaa                                                                                                                                                      >  2:6430789/1‑150 (MQ=255)

tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATTCGTGAGGTGAGGGGACAGCT‑TAGCCTGaa                                                                                                                                                      >  2:3163627/1‑150 (MQ=255)

     tGCACCGAAACGGGCCGGACTTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGAc                                                                                                                                                 <  2:5444305/150‑1 (MQ=255)

     tGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGAc                                                                                                                                                 <  1:1470480/150‑1 (MQ=255)

     tGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGAc                                                                                                                                                 <  1:1470676/150‑1 (MQ=255)

     tGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGAc                                                                                                                                                 <  2:244196/150‑1 (MQ=255)

     tGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGAc                                                                                                                                                 <  1:5790312/150‑1 (MQ=255)

     tGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGAc                                                                                                                                                 <  2:4340418/150‑1 (MQ=255)

        aCCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCGCGCATGCGTGAGGTCAGGGGCCAGCT‑TAGCCTGAACTGACGCa                                                                                                                                              >  2:6778825/1‑150 (MQ=255)

             aaCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTcgc                                                                                                                                         <  1:2381229/150‑1 (MQ=21)

              aCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTcgcg                                                                                                                                        <  1:6001794/150‑1 (MQ=21)

              aCGGGCCGGAATTGGGTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTcgcg                                                                                                                                        <  2:4766572/150‑1 (MQ=21)

               cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                                                       >  2:3578442/1‑150 (MQ=255)

               cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCgg                                                                                                                                       <  1:4371051/150‑1 (MQ=21)

               cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCgg                                                                                                                                       <  1:2184440/150‑1 (MQ=21)

               cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCGGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                                                       <  1:1706720/150‑1 (MQ=255)

                gggCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCggg                                                                                                                                      >  2:5997961/1‑150 (MQ=255)

                 ggCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCGGGc                                                                                                                                     <  1:5015328/150‑1 (MQ=21)

                   ccGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCaa                                                                                                                                   <  2:2287833/150‑1 (MQ=255)

                   ccGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCaa                                                                                                                                   <  2:1996520/150‑1 (MQ=255)

                   ccGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCGGGCaa                                                                                                                                   <  1:342006/150‑1 (MQ=21)

                    cGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCG‑TAGCCTGAACTGAAGCATTAGCGGGCAA‑g                                                                                                                                 >  2:1444565/1‑150 (MQ=255)

                      gAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑Gcc                                                                                                                               <  2:1870415/150‑1 (MQ=255)

                            tttCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAAATGAAGCATTCGCGGGCAA‑GCCCGCGcc                                                                                                                         >  2:699231/1‑150 (MQ=255)

                            tttCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCCAATTCGGGAGGTGAGGGGACAGCT‑TATCCGGAACGTACGCAATCGCGGGCAA‑GCCCGCGcc                                                                                                                         >  1:330212/1‑150 (MQ=255)

                             ttCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGACCGCCACTGCGGGAGGTCAGGGGGCAGCG‑TAGCCGGAACCGACGCAGTCGCGGgccc‑gcccccaccg                                                                                                                        >  1:6488073/1‑144 (MQ=255)

                               tcaccggcCCCGCACAGGTTCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCCCGCATGCGTGTGGTCAGGGGGCAGCT‑TTGCCTGAACAGAAGCATTCGCGGGCAA‑GCCCCCTCCGGc                                                                                                                      >  1:3707524/7‑150 (MQ=255)

                                 tGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCCCGAATGCGTGAGGTCAGGGGACAGCG‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAg                                                                                                                    >  2:738482/1‑150 (MQ=255)

                                  ggCGCCCCGAACCTGGGGGGGGCCGGTGACAGAACAGTACCCGATGATTGTGGGGTAACAATGCAAACGGGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTGGCGGGCAA‑GCCCGCTCCTGCAgg                                                                                                                   <  2:6738777/150‑1 (MQ=255)

                                   gcgcCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCCCGAAAGCGTGACGTGAGGGGACAGCT‑TAGCCTGAACTGACGCAATAGCGGGGAA‑GCCCGCGCCAGCAggg                                                                                                                  >  1:4142424/1‑150 (MQ=255)

                                     gcCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCCCAAGCCTGAGGTCAGGGGACAGCG‑TAGCCTGAACTGACGCATTCGCGGGCCA‑GCCCGCGCCTGCAggggg                                                                                                                >  2:4197127/1‑150 (MQ=255)

                                        ccGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTTCCCGCGCAAGGGGGAGGTCCGGGGCCAGCG‑TAGGCTGAACAGAACCATTAGCGGGCAA‑GCCCGCGCCAGCCggggggcg                                                                                                             >  1:621663/1‑147 (MQ=255)

                                        ccGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCCCCAAAGCGGGAGGTCAGGGGACCCCG‑TAGCCGGCACAGACGCAATCCCGGGCCA‑GCCCGCGCCAGCAGGGGGTAt                                                                                                             >  1:6162449/1‑150 (MQ=255)

                                             ccTGTGCGGGGCGGGTGACAGAACATTAACCGATGATTGTCGGTTAAAAATGCAAACGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCg                                                                                                        <  2:6229551/150‑1 (MQ=39)

                                              cgggggggggCCGGTCACACAACATTAACCTTGGATGGTCGGTTAACACTGCAAATGGGGGGTGGCCGGCTATGGCTTGTGGTCCGGGGACAGGT‑TCGCCTGAACTGACGCATCCGGGGGCAC‑CCCCGCTCCTGCAGGGGGTATTGCCgg                                                                                                       <  2:2175859/144‑1 (MQ=255)

                                                    ggggCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCCCGCATGCGGGCGGGCCGGGGGCCGCG‑TAGCCGGGACAGAAGCAGTCGCGGGCCA‑GCCCCCGCCGGCAGGGGGGATTGCCggggcggg                                                                                                 >  1:3001811/1‑145 (MQ=255)

                                                     gggCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTTGCCGCGAATGCGTGAGGGGCGGGGGCAGCG‑TAGGCGGGACTGACGCAGTCGCGGGGAA‑GCCCCCGCCGGCCGGGGGGGTTGCCGGGGAgggg                                                                                                >  1:3626055/1‑149 (MQ=255)

                                                     ggcccGGTTAAAGAACCTTAACCTATTATTGTCGGTAAACAATGAAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTATGGGGACAGAT‑TACCCTGAACTGACGCATTCGCGGGAAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGt                                                                                                <  2:1612100/147‑1 (MQ=255)

                                                        cgggTGACATAACATTAACCGTTGTTTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGTGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCTTTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTa                                                                                             <  2:2740430/148‑1 (MQ=255)

                                                                aGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTTCCCGCCAATTCGTGAGGTCAGGGGGCACCT‑TATCCGGAACGGACGCAATAGCGGGCCA‑CCCCCCGCCTGCAGGGGGTATTGCCGGGGAGGGTTTACGCCGAAg                                                                                     >  1:111352/1‑150 (MQ=255)

                                                                aGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAg                                                                                     >  2:1560024/1‑150 (MQ=255)

                                                                aGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCGGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGGATTGCCGGGCAGGGGTTACGCCGAAg                                                                                     >  2:6501092/1‑150 (MQ=255)

                                                                    ctgTAACCGATGATGGTCGGTTAACAATGCAAATGCGGGCTGGCCCGCGATTGCGTGTGGCCTGGGGACAGCT‑TAGCCTGAACTGAGGCATTCGGGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAg                                                                                 <  2:1057323/148‑1 (MQ=255)

                                                                       tAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGtt                                                                              >  2:3837473/1‑150 (MQ=255)

                                                                       tAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGtt                                                                              >  2:5042795/1‑150 (MQ=255)

                                                                       tAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCAATCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGtt                                                                              >  2:518363/1‑150 (MQ=255)

                                                                       tAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCCCCAATGCCTGAGGTCAGGGGACACCT‑TAGCCGGAACTGACGCAATCGCCGGCAA‑GCCCCCGCCAGCACGGGGTAGTGCCCGGCAGGGTTTACGCCGAAGCCacggt                                                                              >  2:4620677/1‑146 (MQ=255)

                                                                         aCCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCa                                                                            >  2:6090232/1‑150 (MQ=255)

                                                                           cGATGATTGTCGGGTAACAATGCAAAAGAGGGCTTGCCCGCGAATGCGTGTGGTCAGGGGACAGCA‑TAGCCAGAACTGACGCATTCGCGGGCAA‑GCCCGCACCAGCAGGGGGTCTGGCCGGGAAGGGTTTACGACGAAACCAGGGTCATc                                                                          >  1:5679547/1‑150 (MQ=255)

                                                                            gatgatTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCAtag                                                                         >  2:933672/1‑148 (MQ=255)

                                                                            gatgatTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCG‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCCGGTTCATCg                                                                         >  2:1484756/1‑150 (MQ=255)

                                                                            gatgatTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGGGAGGTCAGGGGACAGCT‑TAGCCTGAACTGAAGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTAAGCCGAAGCCAGGTTAAtag                                                                         >  2:229483/1‑148 (MQ=255)

                                                                                atTGTCGGGTACCACCCCAAACGCGGGCTGGCCCGCTATTGGGTGCGGTCAGGGGACAGCT‑TCGCCTGAACTCCCCCATTCGGGGGCAA‑CCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCtt                                                                     <  2:2639282/150‑1 (MQ=255)

                                                                                atTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCCGGGGGTATTGCCGGGGAGGGTTTACGCCGAAGCCCGGTTCATCGCCtt                                                                     >  2:4092368/1‑150 (MQ=255)

                                                                                  tGTCGGGAAACAATGAAAACGCGGGCTGGCCGGCGATTGCGTGTGTTCCGGGGACCGGT‑TCGCCTGACCTGACGAATGCGCGGGCAC‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTtgt                                                                   <  2:3922413/150‑1 (MQ=255)

                                                                                   gtgggttAACAATGCAAACGGGGGCTTGCCCGCGAATGCGTGCGGTCTGGGGACAGCT‑TAGCCTCAACGGACGAATCCGGGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTtgtg                                                                  <  2:4347900/144‑1 (MQ=255)

                                                                                    tCGGTTAACATTGCAAACGCGGGCTTGCCCGCTAATGCGTTAGGTCTGGGTACAGAT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTgg                                                                 <  2:2119238/150‑1 (MQ=255)

                                                                                       ggTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAAATGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTa                                                              >  2:2202833/1‑150 (MQ=255)

                                                                                                 cAAAAGCGGGCTTGCCCCCGAATGCGTGATGTCAGGGGACAGCG‑TAGCCGGCAAATACGCATTAGCGGGCAA‑GCCCCCACCATCCGGGGGTATTTCCGGGAAGGGTTTACGCCGAACCCACGTTCATAGCCTTGTGGGTATCGATCATGt                                                    >  2:5170096/1‑150 (MQ=255)

                                                                                                  attcgCGGGCAAGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTg                                                   >  2:6105896/5‑150 (MQ=255)

                                                                                                  aaaaGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTg                                                   >  2:3707831/1‑150 (MQ=255)

                                                                                                          gCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCGGGCAA‑GCCCGCTTCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACca                                           >  2:1024597/1‑150 (MQ=18)

                                                                                                            ttGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACcaca                                         >  1:966987/1‑150 (MQ=255)

                                                                                                               cccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACAccc                                      >  2:5131755/1‑150 (MQ=255)

                                                                                                               cccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACAccc                                      >  2:5131677/1‑150 (MQ=255)

                                                                                                               cccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACAccc                                      >  1:1980240/1‑150 (MQ=255)

                                                                                                               cccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACAccc                                      >  2:4581113/1‑150 (MQ=255)

                                                                                                               cccGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACcaacaag                                      >  1:4882449/1‑145 (MQ=18)

                                                                                                                cggcttaTGGGTGTGGTCTGGGGACTGAT‑TAGCCTGAACCCACGCATTCGCGGGAAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCg                                     <  2:5318471/144‑1 (MQ=255)

                                                                                                                ccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCg                                     >  1:5949125/1‑150 (MQ=255)

                                                                                                                 cgcgAATGCGTGAGGTCTGGGGACAGCT‑TCGCCTGAACTGCCGCTTTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCgg                                    <  2:6588364/150‑1 (MQ=255)

                                                                                                                  gcgAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCggg                                   >  1:3963743/1‑150 (MQ=255)

                                                                                                                   gttttGCTTGTGGTCGGGGGCCTGCT‑TCGCCCCAAGTGCCGC‑TTCCGGGGCAACCCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGt                                  <  1:1061555/146‑1 (MQ=255)

                                                                                                                    gAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTc                                 >  1:2014613/1‑150 (MQ=255)

                                                                                                                        gCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAgg                                                                                                  <  2:550024/81‑1 (MQ=255)

                                                                                                                        gCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAgg                                                                                                  >  1:550024/1‑81 (MQ=255)

                                                                                                                         cGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCa                            >  2:2893058/1‑150 (MQ=255)

                                                                                                                               gTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCg                      >  2:2404806/1‑150 (MQ=255)

                                                                                                                                tccgggcACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCga                     <  1:6778825/147‑1 (MQ=255)

                                                                                                                                tCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCga                     >  1:4518366/1‑150 (MQ=255)

                                                                                                                                tCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCga                     >  2:2965814/1‑150 (MQ=255)

                                                                                                                                tCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCga                     >  2:4841083/1‑150 (MQ=255)

                                                                                                                                tCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCga                     >  2:6100636/1‑150 (MQ=255)

                                                                                                                                tCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCAAAGGTCga                     >  2:4840930/1‑150 (MQ=255)

                                                                                                                                 ctgggtACAGCT‑TCGCCACAACTCACGCTTTCGCGGGCAA‑CCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCgag                    <  1:4425377/148‑1 (MQ=255)

                                                                                                                                    gggACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGata                 >  1:1727633/1‑150 (MQ=255)

                                                                                                                                         aGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGcc            >  2:1859506/1‑150 (MQ=255)

                                                                                                                                         aGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGcc            >  2:5611299/1‑150 (MQ=255)

                                                                                                                                         aGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGcc            >  2:2506915/1‑150 (MQ=255)

                                                                                                                                          gCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTgg                                                                 >  1:3106697/1‑96 (MQ=255)

                                                                                                                                          gCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTgg                                                                 <  2:3106697/96‑1 (MQ=255)

                                                                                                                                           cT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCa          >  1:2864388/1‑150 (MQ=255)

                                                                                                                                           cT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCa          >  1:2465942/1‑150 (MQ=255)

                                                                                                                                           cT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCa          >  1:5118258/1‑150 (MQ=255)

                                                                                                                                            t‑tAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAg         >  1:2689172/1‑150 (MQ=255)

                                                                                                                                               aGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGATGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAggg       >  2:5010776/1‑150 (MQ=255)

                                                                                                                                                gCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGc      >  2:3669378/1‑150 (MQ=255)

                                                                                                                                                  cttaaCTGACGAATTCGAGGGCAA‑GCCAGCTCCTTCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGACAGGTTCAACGCAATGAGGGTATCGATCACGAGCTGAACCACACACGGGAAAGCCAAGCTCGAGATATCTCCCAGGGcga    >  2:2190189/4‑149 (MQ=255)

                                                                                                                                                   tGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTACACCGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcgtc   >  1:5077806/1‑150 (MQ=255)

                                                                                                                                                    gAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcgtcg  >  1:2660788/1‑150 (MQ=255)

TCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGCGTCG  >  NC_002947/5102095‑5102391

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