BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,099,220	G→C	7.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,099,220	0	G	C	7.6%	376.3 / 11.8	125	intergenic (+241/+51)	PP_mr66/acsA‑I	P6/acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand): ref base G (82/33); new base C (2/8); total (84/41)
Fisher's exact test for biased strand distribution p-value = 2.05e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

TGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGCGTCGTACT  >  NC_002947/5099073‑5099368

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

  cACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGc                                                                                                                                                  <  2:762078/150‑1 (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)

        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‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTA‑AGCTGTCCCCTGACCTCACGCATTcgcg                                                                                                                                           <  2:2567820/150‑1 (MQ=21)

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

          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)

           gggCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGGTTGCCCCCGAAAGCGTGAGGTCAGGGGACAGCT‑TAGCCTGTACTGACGCAATAGCggg                                                                                                                                         >  1:939777/1‑150 (MQ=255)

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

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

                 gcaTTGTGTCCGGGCGCCCCGCACCTGCGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGGCAGGGGACAGCT‑AAGCCTGAACTGACGCATTCGCGGGCAAGcc                                                                                                                                   <  1:2367763/148‑1 (MQ=255)

                 gAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGGGTGCCCGCGCATGCGTGGGGGCCGGGGCCAGCT‑TAGCCGGGACGGACGCAGTGGCGGGCCAGcc                                                                                                                                   >  1:1752230/1‑150 (MQ=255)

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

                    ttttttCCTGGCGCCCCTAACGTGTGGGGGGCCGGTG‑ACAGAACTTTAACCTTTGATTGTGGGTTAACAATCCAAAAGCGGGCTTGCCCGCGAATGCTTGAGGTCTGGGGACTGCT‑TAGCCTTAACTGACCAATTCGCGGGCAAGCCCGc                                                                                                                                <  2:6871002/147‑1 (MQ=255)

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

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

                             ggCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAAAGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAgg                                                                                                                       >  2:6455859/1‑150 (MQ=255)

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

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

                                  cccGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCCCGCAAGCGGGAGGTGCGGGGGCAGCG‑TAGCCGGAAAAGACGCAATAGCGGGCAAGCCCGCGCCAGCAgggggga                                                                                                                  >  1:1579263/1‑148 (MQ=255)

                                  cccGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCATGCCCCCCAATGCGGGAGGTGAGGGGACACCG‑TAGCCTGAAAAGACGCAATCCCGGGCAACCCCCCACCAGCAgggggga                                                                                                                  >  1:1354510/1‑148 (MQ=255)

                                   ccGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCGTGCCCGCGAATGCGGGAGGTCCGGGGACAGCT‑TAGCCGGAACTGACGCAATCGCGGGCAAGCCCGCGCCTGCCGGGGGTAt                                                                                                                 >  1:3567003/1‑150 (MQ=255)

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

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

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

                                                                cAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCCAATGCGTGAGGTCAGGGGACAGCG‑TAGCCTGTACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAg                                                                                     >  1:8330205/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)

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

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

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

                                                                                  ggttAACAATCAAAACGGGGGCTGGCCCGCCAAGGCTTGAGGCCTGGGGACATTT‑TAGCACCAACACCCGCTTTCGGGGCAAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGt                                                                   <  2:1649419/147‑1 (MQ=255)

                                                                                   gttAACATTGAAAACGCGGGCTGGCCGGCTAATGCTTGTGGTCTGGGTACAGCT‑TAGCCTCACCTGCCGCATTCGGGGCAAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTa                                                                  <  2:5911618/148‑1 (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)

                                                                                             cAAACGCGGGGTTGCCGGCGATGGGGTGGGGTCTGGGTACCGCT‑TCGCCTGACGTGCCGCTTGCGCGGCCACCCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGt                                                        <  1:9025124/150‑1 (MQ=255)

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

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

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

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

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

                                                                                                    gggCTTGCCCGCGACTGCGTGCGGTCCGGGGACAGCT‑TAGCCTGACCTCACGCATTCGCGGGCACGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTAc                                                 <  2:8981516/150‑1 (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)

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

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

                                                                                                            ccGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCg                                         >  1:2657959/1‑150 (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)

                                                                                                             cgctactgcgtgcgGTCCGGGGACCGCT‑TCGCCTGAACTGCCGCTTGCGCGGGCACGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCgg                                        <  2:3595115/144‑1 (MQ=255)

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

                                                                                                              gcgAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCggg                                       >  2:3381899/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                                                                   >  2:8156293/1‑114 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                                                                              ctccccTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcgt        <  1:3788888/145‑1 (MQ=255)

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

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

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

                                                                                                                                                 aaCTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGcgtcgt     >  1:6351000/1‑150 (MQ=255)

                                                                                                                                                    tGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGCGTCGTACt  >  2:1836539/1‑150 (MQ=255)

TGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTG‑ACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCT‑TAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGGGCGTCGTACT  >  NC_002947/5099073‑5099368

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