BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,102,229	G→T	10.2%	intergenic (+223/+69)	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,229	0	G	T	10.2%	415.8 / 26.7	143	intergenic (+223/+69)	PP_mr66/acsA‑I	P6/acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand): ref base G (78/49); new base T (2/14); total (80/63)
Fisher's exact test for biased strand distribution p-value = 2.54e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 150 total aligned reads displayed.

CACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATC  >  NC_002947/5102081‑5102378

cacaAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGa                                                                                                                                                        <  1:6377851/150‑1 (MQ=255)

cacaAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGGGAGGTCAGGGGa                                                                                                                                                        >  2:5278987/1‑150 (MQ=255)

   aaaaTGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAg                                                                                                                                                     <  1:1734860/150‑1 (MQ=255)

     aaTGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAATGCGGGCTTGCCCGCGAATGCGGGAGGTCAGGGGGCAGCt                                                                                                                                                   >  2:1704489/1‑150 (MQ=255)

     aaTGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCCCGAAGGCGGGAGGTCAGGGGACAGcg                                                                                                                                                   >  2:5051086/1‑149 (MQ=255)

       tGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑Ta                                                                                                                                                <  1:5364165/150‑1 (MQ=255)

          tGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGcc                                                                                                                                             <  1:6114887/150‑1 (MQ=255)

          tGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGcc                                                                                                                                             <  1:6114791/150‑1 (MQ=255)

          tGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGcc                                                                                                                                             <  2:4914340/150‑1 (MQ=255)

            aCTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTg                                                                                                                                           <  1:5936346/150‑1 (MQ=255)

              tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTTCCCGCGGATGCGTGAGGTCAGGGGACAGCT‑TAGCCTgta                                                                                                                                         >  2:3163541/1‑148 (MQ=255)

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

              tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGaa                                                                                                                                         <  2:2275505/150‑1 (MQ=255)

              tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGaa                                                                                                                                         <  2:6264366/150‑1 (MQ=255)

              tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCGCGCATGCGTGAGGTCAGGGGACAGCG‑TAGccggga                                                                                                                                         >  2:760574/1‑146 (MQ=255)

              tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTCCGGGGCC‑GGTGACAGAACAGTAACCGATGACAGACGGGTAACAATGCAATTGCGGCCGTTCCCGCGAAAGCCTGAGGGCCCGGGACAGCT‑AAGccataa                                                                                                                                         >  2:4034613/1‑146 (MQ=255)

                cGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGGCCGCGCAAGCGTGAGGTCAGGGGACAGCT‑TGTCCGGAACt                                                                                                                                       >  1:2328048/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                             cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTCGCgg                                                                                                                          <  2:2122742/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:1996520/150‑1 (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)

                                   ggAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTAGCGGGCAA‑Gc                                                                                                                   >  2:1450492/1‑150 (MQ=255)

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

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

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

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

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

                                                   gcCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCCCAAGCCTGAGGTCAGGGGACAGCG‑TAGCCTGAACTGACGCATTCGCGGGCCA‑GCCCGCGCCTGCAggggg                                                                                                   >  2:4197127/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)

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

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

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

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

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

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

                                                                               aGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCG‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTCTTGCCGGGCAGGGTTTAAGCCGTAg                                                                        >  2:6500934/1‑150 (MQ=39)

                                                                               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‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGtt                                                                 >  2:3837473/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)

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

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

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

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

                                                                                                   tCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTgg                                                    >  2:6791274/1‑150 (MQ=255)

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

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

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

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

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

                                                                                                                              cccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACAccc                         >  2:5131677/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)

                                                                                                                                 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)

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

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

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

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

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

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

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

                                                                                                                                                     gACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATc  >  1:1855729/1‑150 (MQ=255)

                                                                                                                                                     gACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATc  >  2:762627/1‑150 (MQ=255)

CACAAAATGCTGACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCC‑GGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAA‑GCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATC  >  NC_002947/5102081‑5102378

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