BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,102,247	G→C	9.6%	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	9.6%	438.6 / 28.8	155	intergenic (+241/+51)	PP_mr66/acsA‑I	P6/acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand): ref base G (94/43); new base C (2/13); total (99/56)
Fisher's exact test for biased strand distribution p-value = 4.59e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

Only 100 of 156 total aligned reads displayed.

TGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGCGTCGTACTC  >  NC_002947/5102100‑5102396

tGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGAc                                                                                                                                                      <  2:4340414/150‑1 (MQ=255)

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

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

   aCCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCa                                                                                                                                                   <  2:3896620/150‑1 (MQ=255)

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

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

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

          cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                                                            <  1:6425583/150‑1 (MQ=255)

          cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                                                            <  2:3673350/150‑1 (MQ=255)

          cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                                                            <  2:4732994/150‑1 (MQ=255)

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

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

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

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

             gCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGCATGCGGGTGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCa                                                                                                                                         >  1:2175859/1‑150 (MQ=255)

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

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

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

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

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

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

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

                          ccTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGc                                                                                                                           >  2:5610052/1‑150 (MQ=255)

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

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

                                gcCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCG‑GAGCCGGAACTTAAGCAGTCGCGGGCAA‑GCCCCCACCAGCAggggg                                                                                                                     >  1:3922413/1‑150 (MQ=255)

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

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

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

                                            tgCGGGGCC‑GGTGACAGAACAGTACCCGATGTTTGGCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGGGTGAGGTCAGGGGACAGCT‑TAGCCTGAACGGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGaa                                                                                                         <  2:3952833/150‑1 (MQ=255)

                                            tgCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGGTTTCCCGCGAAAGCGTGGGGTGAGGGGACAGCT‑TTGGCGGAACTGACGCATTCGCGGGCCA‑GCCCGCGCCGGGCGGGGGGGTTGCCggggc                                                                                                         >  1:4624211/1‑148 (MQ=255)

                                               ggggCCGGGTCACACAACTGTACCGGATGTTTGTCGTTTAAAACTGCAAACGGGGGCTTGCCGGCCATTGGGTGTGGTACGGGTACAGCT‑TCGCCTGAAGGGACGATTTGGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAgg                                                                                                       <  2:1503259/150‑1 (MQ=255)

                                               ggggCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATTCGTGAGGTCAGGGGACAGCT‑TATCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGGGTTGCCGGGAAggg                                                                                                      >  2:1893176/1‑150 (MQ=255)

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

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

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

                                                 ggCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGCCAGCT‑TAGCCTGAACAGACGCATTCGCGGGCAA‑GCCCGCTCCTTCAGGGGGGATTGCCGGGAAGGGtt                                                                                                    >  2:3064226/1‑150 (MQ=255)

                                                   cc‑GGTTACACAACAGTACCCGATGATTGTCGGTTAAAAATGCAAACGGGGGCTGGCCCGCTAATGCTTGAGGTCAGGGGACAGCT‑TCGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTa                                                                                                  <  2:2740326/150‑1 (MQ=255)

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

                                                            aGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGCAAGCGTGAGGTCAGGGGACAGCT‑TAGCCGGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGGATTGCCGGGACGGGTTTACGCCGAAg                                                                                          >  1:5057228/1‑150 (MQ=255)

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

                                                            aGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTTCCCCCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCAGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCCGGGGGTATTTCCGGGAAGGGTTTAAGCCGAAg                                                                                          >  2:24955/1‑150 (MQ=255)

                                                                ctgTACCCTATGTTTGTCGGGTAACAATGCAAACGCGGGCGTGCCCGCGATTGCGTGCGTTCAGGTTACAGCT‑TCGCCTGACCTGCCGCTTTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAg                                                                                      <  2:1057310/148‑1 (MQ=255)

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

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

                                                                     aCCGTTGTTTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGGGTGTGGTCTGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAC‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCa                                                                                 <  2:2805912/150‑1 (MQ=255)

                                                                      ccGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAAATCGTGAGGTAAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCCA‑GCCCCCTCCTGCAGGGGGTATTGCCGGGAAGGGGTTAAGCCGCAGCCAGGTTCaa                                                                                >  2:1834126/1‑149 (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)

                                                                        gatgatTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAAAGCGTGGGGGCAGGGGACAGCG‑TAGCCTGAACTGACGCATTCGCGGGCAA‑CCCCGCTCCTGCAGGGGGTATTGCCGGGGAGGGTTTACGCCGAAGCCAGGTTCaaag                                                                              >  2:229852/1‑147 (MQ=255)

                                                                            ttggTCGGTAAACACTGCAACCGGGGGCTGGCCCGCTATTGGGTGCGGTCCGGGGACCGCT‑TCGCCTTACCTGCCGCATTCGGGGGCAC‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCtt                                                                          <  2:6062555/147‑1 (MQ=255)

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

                                                                               gTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGTCAGCT‑TAGCCTGAACGGTCGCGTTCGCGGGCAA‑GCCCGCTGCGGCAGGGGGTGTTGCCGGGCAGGGTTTACGCCGAAGCCAGGGTAATCGCCTtgtg                                                                       >  2:5357049/1‑150 (MQ=255)

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

                                                                                   gttaaataTGCAAAAGGGGGTTTGCCCGGTTTTGCTTGTGTTCGGGTTAATGAT‑TCGCCTCAAATGCCGAATTCGGGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTa                                                                   <  1:5711540/143‑1 (MQ=255)

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

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

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

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

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

                                                                                                           cccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACAccc                                           >  1:4213730/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:1639150/1‑150 (MQ=255)

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

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

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

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

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

                                                                                                                            tctggggACATCT‑TAGCCTTAACTGACGAATTCGCTGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCATGTGCTGTACCACACCCGTGTCAGCCAAGGtata                          <  2:5054101/147‑4 (MQ=255)

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

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

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

                                                                                                                                  gACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATc                    >  1:1855729/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:2465942/1‑150 (MQ=255)

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

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

                                                                                                                                        t‑tAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTgg                                                                      <  2:446711/94‑1 (MQ=255)

                                                                                                                                        t‑tAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTgg                                                                      >  1:446711/1‑94 (MQ=255)

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

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

                                                                                                                                             ccTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcg          >  2:4767631/1‑150 (MQ=255)

                                                                                                                                              cTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcgt         >  2:2190166/1‑150 (MQ=255)

                                                                                                                                              cTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcgt         >  1:6018682/1‑150 (MQ=255)

                                                                                                                                               tGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcgtc        >  1:6809136/1‑150 (MQ=255)

                                                                                                                                                gcccTGCCGCTTGCGCGGGCAC‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcgtcg       <  1:5051086/147‑1 (MQ=255)

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

                                                                                                                                                    tGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGCGTCGTACt   >  2:4776952/1‑150 (MQ=255)

                                                                                                                                                    tGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGCGTCGTACt   >  1:1201795/1‑150 (MQ=255)

                                                                                                                                                    tGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGCGTCGTACt   >  2:3614033/1‑150 (MQ=255)

TGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGCGTCGTACTC  >  NC_002947/5102100‑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: `‑`