BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,099,217	A→C	8.8%	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,099,217	0	A	C	8.8%	356.1 / 11.5	123	intergenic (+238/+54)	PP_mr66/acsA‑I	P6/acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand): ref base A (81/29); new base C (2/11); total (83/40)
Fisher's exact test for biased strand distribution p-value = 6.79e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 130 total aligned reads displayed.

TCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGCGTCGT  >  NC_002947/5099068‑5099365

tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGaa                                                                                                                                                      <  1:495335/150‑1 (MQ=255)

tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGaa                                                                                                                                                      <  1:1617320/150‑1 (MQ=255)

  cGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTTCCCGCGAATGCGTGAGGTCAGGGGACACCA‑TAGCCGGAACt                                                                                                                                                    >  2:1924054/1‑150 (MQ=255)

     tGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGAc                                                                                                                                                 <  1:2339620/150‑1 (MQ=255)

       cACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGc                                                                                                                                               <  2:762078/150‑1 (MQ=255)

        aCCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGACTGCGTGGGGTCAGGGGCCAGCT‑TAGCCTGGACTGACGCa                                                                                                                                              >  1:6270894/1‑150 (MQ=255)

        aCCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAAAGCGTGCGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCa                                                                                                                                              >  2:3103165/1‑150 (MQ=255)

        aCCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAAAGCGTGAGGTCAGGGGGCAGCT‑TAGGCGGAAAGGACGCa                                                                                                                                              >  1:6362618/1‑150 (MQ=255)

          cGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAATAGCGGGCTTTCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGGCGCAtt                                                                                                                                            >  1:5176460/1‑150 (MQ=255)

           gAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTTCCCGCGAAAGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTc                                                                                                                                           >  1:7418209/1‑150 (MQ=255)

           gAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTc                                                                                                                                           >  1:7418203/1‑150 (MQ=255)

             aaCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGTGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTcgc                                                                                                                                         <  1:2343429/150‑1 (MQ=255)

             aaCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTcgc                                                                                                                                         <  1:8884620/150‑1 (MQ=21)

             aaCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTcgc                                                                                                                                         <  1:7322925/150‑1 (MQ=21)

              aCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAGGAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTcgcg                                                                                                                                        <  2:525500/150‑1 (MQ=18)

              aCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTcgcg                                                                                                                                        <  2:6637919/150‑1 (MQ=255)

              aCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTcgcg                                                                                                                                        <  2:2567820/150‑1 (MQ=21)

               cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                                                       <  2:4582175/150‑1 (MQ=255)

               cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                                                       <  2:7001024/150‑1 (MQ=255)

               cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                                                       <  2:7063782/150‑1 (MQ=255)

               cGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCgg                                                                                                                                       <  1:1545442/150‑1 (MQ=255)

                 ggCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTGAGGGGCCAGCT‑TAGCCGGAACTGACGCATTCGCGGGc                                                                                                                                     >  1:4021283/1‑150 (MQ=255)

                   ccGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGGCAGCT‑TATCCTGGACAGACGCCTTCGCGGGCaa                                                                                                                                   >  1:1739008/1‑150 (MQ=255)

                      gAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAAAAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGcc                                                                                                                                <  1:2367655/150‑1 (MQ=255)

                       atttGTTTCCTGGCGCCCCGAACCTGTGGGGGGCCGGTT‑AAATAACAGTAACCTATGATTGTCGGTTAACAATGCAAAAGCGGGCTTGCCCGCTATTGCGTGAGTTCTGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGccc                                                                                                                               <  1:1161979/148‑1 (MQ=255)

                                cTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCa                                                                                                                      >  2:3777522/1‑150 (MQ=255)

                                 tGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAg                                                                                                                     >  2:3964524/1‑150 (MQ=255)

                                   gcgcCCCGAACCTGTGCGGGG‑CGGTGTACTGAACAGTACCCGTTTCTTGGGGGTTAACATTGCAACCGCGGGCTTGCCCGCCATTGCTTGTGGTCTGGGTACAGCT‑TAGCCTGAACTGCCGATTTCGGGGGCACGCCCGCTCCTGCAggg                                                                                                                   <  2:2961975/150‑1 (MQ=255)

                                      ccccGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACGGACGCATTCGCGGGCCAGCCCCCTCCTGCAGGGGGt                                                                                                                >  2:6212376/1‑150 (MQ=255)

                                        ccGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCGCGAATGCGTGTGGGCAGGGGACAGCT‑TAGCCGGAACTGACGCAGTCGCGGGCCAGCCCGCGCCGGCAggggggct                                                                                                              >  2:1570044/1‑147 (MQ=255)

                                                   cGGGGCCGGTG‑ACAGACCTGAAACCGATGATTGTCGGTTAACCATGCAAACGCGGGGTGGCCCGCGAATGGGTGCGGTCTGGGGACCGCT‑TAGCCTGACCGGCCGCATCCGCGGGCACGCCCGCTCCTGCAGGGGGTATTGCCGGGAAgg                                                                                                   <  2:5140372/150‑1 (MQ=255)

                                                            tG‑ACAGAACCGTAACCGTTGATTGTCGGTTAACAATGCAAAAGCGGGCTTGCCCGCGATTGCGTGAGGTCAGGGGACAGCT‑TAGCCTTAACTGACGCTTTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGcc                                                                                          <  2:3668417/150‑1 (MQ=255)

                                                            tG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGGGTTGCCGCCAATTCGTGTGGTGAGGGGGCAGCT‑TAGCCGGGACTGACGCAGTAGCGGGCAAGCCCGCGCCTGCCGGGGGGGTTGCCGGGGAGGGGTTACGcc                                                                                          >  1:919833/1‑150 (MQ=255)

                                                                 aGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAAAGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAg                                                                                      >  2:8049199/1‑150 (MQ=255)

                                                                 aGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTTCCCCCGAAAGCGTGAGGTCCGGGGACAGCT‑TAGCCAGAACTGACGCATTCGCCGGCAAGCCCGCACCTGCCGGGGGTATTGCCGGGAAGGGTTTACGCCGAAg                                                                                      >  1:3550716/1‑150 (MQ=255)

                                                                            cGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCCCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCGCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCCAAGCCAGGTTCATc                                                                           >  2:7860457/1‑150 (MQ=255)

                                                                             gatgatTGTCGGGTAACAATGCAAAAGCGGGCTTTCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTAAGCCCAAGCCAGGTTCATCg                                                                          >  2:8471386/1‑150 (MQ=255)

                                                                             gatgatTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCg                                                                          >  2:4562511/1‑150 (MQ=255)

                                                                             gatgatTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGGTTACGCCGAAGCCAGGTTCATCg                                                                          >  2:8692789/1‑150 (MQ=255)

                                                                             gatgatTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAAAGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAAATGACGCATTCGCCGGCCAGCCCGCTCCTGCAGGGGGGATTGCCCGGAAGGGTTTACGCCGAAGCCCGGTTCAACg                                                                          >  2:8692752/1‑150 (MQ=255)

                                                                             gatgatTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCCCGAAAGCGTGAGGTGGGGGGGCAGCT‑TAGCCTGAACTGACGCAATAGCGGGCAAGCCCGCTCCTGCAGGGGGTATTTCCGGGAACGGTTTAAGCCCAAGCCACGGTCAACg                                                                          >  2:3988608/1‑150 (MQ=255)

                                                                              atgatTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACGGACGCATTCGCGGGCAAGCCCGCGCCGGCAGGGGGTATTGCCGGGAAGGGGTTAAGCCGAAGCCAGGTTAATCGc                                                                         >  1:3735377/1‑150 (MQ=39)

                                                                                        gttAACATTGAAAACGCGGGCTGGCCGGCTAATGCTTGTGGTCTGGGTACAGCT‑TAGCCTCACCTGCCGCATTCGGGGCAAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTa                                                               <  2:5911618/148‑1 (MQ=255)

                                                                                           aaCAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGaa                                                                                                                                                      <  1:1157370/60‑1 (MQ=255)

                                                                                           aaCAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGaa                                                                                                                                                      >  2:1157370/1‑60 (MQ=255)

                                                                                             cAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGAt                                                          >  2:2768384/1‑150 (MQ=255)

                                                                                             cAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCGGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTAGTGCCGGGCAGGGTTTACGCCGAAGCCAGGGTCATCGCCGTGTGGGTATCGAt                                                          >  2:9004293/1‑150 (MQ=39)

                                                                                              aaTGCAAACGGGGGCTGGCCGGCGATTGCTTGTGTTCTGGGGACTGCT‑TAGCCTGAACTGACGCTTTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATc                                                         <  2:1508312/150‑1 (MQ=255)

                                                                                              aaTGCAAAAGCGGGCTTTCCCGCGAAAGCGTGAGGTCAGGGGCCCGCG‑TAGCCGGAACAGCCGCAATCGCGGGCACGCCCGCTCCAGCAGGGGGGCTTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATc                                                         >  2:7848413/1‑150 (MQ=255)

                                                                                               aTGCAACCGGGGGGTGGCCCGCGATTGCGTGCGGTCTGGGGACCGCT‑TCGCCGCAACGGACGCTTCGGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCa                                                        <  1:7571514/150‑1 (MQ=255)

                                                                                               aTGAAAACGGGGGCTTGCCGGCTATGGCTTGTGGTGTGGGGACAGCT‑TCGCCTGACCTGACGCCTGCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCa                                                        <  2:4398748/150‑1 (MQ=255)

                                                                                                  cAAACGCGGGCTTGCCCGCTAATGCGTGAGGTCCGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGt                                                     <  2:680046/150‑1 (MQ=255)

                                                                                                   aaaaGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCCCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTg                                                    >  2:3847352/1‑150 (MQ=255)

                                                                                                     aaGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCt                                                  >  2:3360165/1‑150 (MQ=255)

                                                                                                     aaGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCt                                                  >  1:3651926/1‑150 (MQ=255)

                                                                                                      aGCGGGCTTGCCCGCGAATTCGTGAGGTAAGGGGACAGCT‑TAGCCTGAAATGACGCATTAGCGGGCAAGCCCGCTCCTGCAGGGGGTATTTCCGGGAAGGGTTTTAGCCGAAGCCAGGTTCATAGCCTTGTGGGTATAGATCAGGTGCTg                                                 >  1:1141841/1‑150 (MQ=255)

                                                                                                         gggCTTGCCCGCGACTGCGGGTGGTCCGGGGACTGCT‑TCGCCTGCACGGCCGCTTGCGGGGGCACGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTAc                                              <  2:8981594/150‑1 (MQ=255)

                                                                                                           gCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACca                                            >  2:5467394/1‑150 (MQ=255)

                                                                                                           gCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACca                                            >  1:2729619/1‑150 (MQ=255)

                                                                                                             ttGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACcaca                                          >  1:1429173/1‑150 (MQ=255)

                                                                                                                cccGCTATTGCGTGTGGTCAGGGCACAGCT‑TCGCCTCAAGTGAGGCTTTCGCGGCCAAGCCCGCTCCTGCAGGGGGTATTGCCAGGAAGGGTTTACGCAGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACAccc                                       <  1:9098839/150‑1 (MQ=255)

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

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

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

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

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

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

                                                                                                                   gcgAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCggg                                    >  2:3381899/1‑150 (MQ=255)

                                                                                                                     gAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTc                                  >  2:3070355/1‑150 (MQ=255)

                                                                                                                     gAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTc                                  >  1:6185664/1‑150 (MQ=255)

                                                                                                                      aaggCGTGAGGTCTGGGGACTGTT‑TAGCATGAACTCCCCCTTGCGCGGCCACGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCa                                 <  1:1452647/147‑1 (MQ=255)

                                                                                                                       ctgcgtgcgGTCTGGGGCCCGCT‑TCGCCTGACCTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAg                                <  2:9015764/149‑1 (MQ=255)

                                                                                                                           gTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGt                                                                <  1:8156293/114‑1 (MQ=255)

                                                                                                                           gTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGt                                                                >  2:8156293/1‑114 (MQ=255)

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

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

                                                                                                                            tGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTgg                                                                  <  1:3258492/111‑1 (MQ=255)

                                                                                                                                 tcgggggACTGCT‑TCGCCTGACCTGAGGCTTGGGGGGGCACGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCga                      <  2:1739008/147‑1 (MQ=255)

                                                                                                                                 tccgggcACCGCT‑TACCCTCAAGTCCCGCTTGCGCGGCAAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCga                      <  1:3651754/147‑1 (MQ=255)

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

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

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

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

                                                                                                                                   aGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGaa                                                                                                     >  1:836161/1‑69 (MQ=255)

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

                                                                                                                                    ggggACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGat                   >  2:6017262/1‑150 (MQ=255)

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

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

                                                                                                                                     gggACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGCGCTGTACCACACCCGGGTCAGCCAAGGTCGAGata                  >  2:1845413/1‑150 (MQ=255)

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

                                                                                                                                         ccgCT‑TAGCATCAACTGCCGCTTCCGGGGGCACGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGc              <  2:3440504/148‑1 (MQ=255)

                                                                                                                                           gCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGccc            >  2:404386/1‑150 (MQ=255)

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

                                                                                                                                                 gCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATc                                                         >  1:5025379/1‑100 (MQ=255)

                                                                                                                                                 gCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGc       >  1:9028370/1‑150 (MQ=255)

                                                                                                                                                  ccTGACCTCACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcg      >  2:1147663/1‑150 (MQ=255)

                                                                                                                                                  ccTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcg      >  2:3612068/1‑150 (MQ=255)

                                                                                                                                                  ccTGAACTGACGCATACGCGGGCAATCCCGCTCCTGCAGGGGGTATTGACGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGCGTATCGATCAGGTGCTGTACCACACCCGGGACAGCCAAGGTCGAGATATAGCCCAGGGcg      >  1:222076/1‑150 (MQ=255)

                                                                                                                                                   cTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcgt     >  1:833167/1‑150 (MQ=255)

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

TCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGCGTCGT  >  NC_002947/5099068‑5099365

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