BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,337,825	G→C	6.5%	G329A (GGG→GCG)	PP_2934 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,337,825	0	G	C	6.5%	398.1 / 15.2	123	G329A (GGG→GCG)	PP_2934	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (58/57); new base C (5/3); total (63/60)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.50e-01

Only 100 of 123 total aligned reads displayed.

GGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACC  >  NC_002947/3337680‑3337968

gggCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCg                                                                                                                                              >  2:160559/1‑150 (MQ=255)

  gCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGAAGCCAGACACGACAGCGTTCTTCAGCCAGACGCTGGGGTCGcc                                                                                                                                            >  2:1123934/1‑150 (MQ=255)

       tCGAATGCATTGCCCATTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAgcgc                                                                                                                                       <  2:262484/150‑1 (MQ=255)

        cGTATTCATGGCCCAAGGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGNAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGCGACCCCAGCGcc                                                                                                                                      <  1:1442420/150‑1 (MQ=255)

          tATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑gg                                                                                                                                   <  2:564013/150‑1 (MQ=255)

              cATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGa                                                                                                                               <  2:3089/150‑1 (MQ=255)

                  gCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGt                                                                                                                           >  2:1318422/1‑150 (MQ=255)

                     cACTGTCATGCCGAGCAGTACGCGCGAATGCCGGCGCTGGGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGTGCGTTCa                                                                                                                        <  2:1344389/150‑1 (MQ=255)

                     cACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCa                                                                                                                        <  1:743448/150‑1 (MQ=255)

                      aCTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑g                                                                                                                                    >  1:1599488/1‑137 (MQ=255)

                      aCTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAg                                                                                                                       <  1:1123934/150‑1 (MQ=255)

                       cTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGc                                                                                                                      <  2:2036110/150‑1 (MQ=255)

                         gCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTg                                                                                                                    <  1:228855/150‑1 (MQ=255)

                             tGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGATGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGAc                                                                                                                >  2:1817228/1‑150 (MQ=255)

                                 gAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGgc                                                                                                            >  1:1089931/1‑150 (MQ=255)

                                       tatgcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGGC‑GGCTGAGCGTTCAGCTGCGACGGGcgcgccg                                                                                                      >  1:953352/4‑150 (MQ=255)

                                        aCGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGcgcgccgc                                                                                                     <  2:1089931/150‑1 (MQ=255)

                                          gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCa                                                                                                                        <  1:43016/129‑1 (MQ=255)

                                          gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCa                                                                                                                        >  2:43016/1‑129 (MQ=255)

                                          gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGcgcgccgcgc                                                                                                   >  2:899390/1‑150 (MQ=255)

                                          gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGcgcgccgcgc                                                                                                   >  1:1799052/1‑150 (MQ=255)

                                              gAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCgccagc                                                                                               >  1:649909/1‑150 (MQ=255)

                                                 tGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCcac                                                                                            >  2:222648/1‑150 (MQ=255)

                                                   ccGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCTATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACaa                                                                                          <  1:2068151/150‑1 (MQ=255)

                                                   ccGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACaa                                                                                          >  2:1503000/1‑150 (MQ=255)

                                                    cGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAAc                                                                                         >  2:1630088/1‑150 (MQ=255)

                                                     ggCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACg                                                                                        >  1:915163/1‑150 (MQ=255)

                                                          tgtgCAAGGACGCCGGCGGAGACCAGCCGCTGTATATCCATCATTAGCTACAGCGCGAACCCACAGCTGGCTGCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAgg                                                                                   <  1:44472/150‑1 (MQ=255)

                                                          tgtgCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAgg                                                                                   <  1:1318422/150‑1 (MQ=255)

                                                            tgCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGcc                                                                                 <  2:72298/150‑1 (MQ=255)

                                                                tcgACGCCGTCGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGAcc                                                                             >  2:2325104/3‑150 (MQ=255)

                                                                tcgACGCCGTCGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGACGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACGACGCTAGGCCGAcc                                                                             >  2:1126018/3‑150 (MQ=255)

                                                                 ggACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGAccc                                                                            >  2:1362265/1‑150 (MQ=255)

                                                                       cGACGGCGTCGACCGCCGGTACATACACGAATCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCa                                                                      >  1:65048/1‑150 (MQ=255)

                                                                        gACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCaa                                                                     >  1:2032291/1‑150 (MQ=255)

                                                                         aCGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCggactac                                                                    <  1:1069523/150‑6 (MQ=255)

                                                                         aCGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAg                                                                    <  2:654210/150‑1 (MQ=255)

                                                                              gTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAg                                                               >  2:1867665/1‑150 (MQ=255)

                                                                               tCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGACAAGGCGAGc                                                              >  1:1844670/1‑150 (MQ=255)

                                                                                  aCCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCaaa                                                           <  2:1844670/150‑3 (MQ=255)

                                                                                   ccgccgGTATATCCATACTTCGCTGCAGCGAGAACCGGCAGCGTTCTTCAGCCAGGCGCTGGGGTGCCCAGCGGG‑GGCTGGGCGTACACCCGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAggg                                                          <  2:2407781/150‑1 (MQ=255)

                                                                                    cgccgGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGc                                                         >  2:989937/1‑150 (MQ=255)

                                                                                      ccgGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAg                                                                                    >  2:1324228/1‑121 (MQ=255)

                                                                                      ccgGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAg                                                       >  2:1698564/1‑150 (MQ=255)

                                                                                       cgGTATCTCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGc                                                      <  1:1565677/150‑1 (MQ=255)

                                                                                       cgGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGc                                                      <  2:915163/150‑1 (MQ=255)

                                                                                        ggTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGcc                                                     <  2:836347/150‑1 (MQ=255)

                                                                                          tataTCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCAc                                                   <  2:410280/150‑1 (MQ=255)

                                                                                             atCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCa                                                <  2:1521198/150‑1 (MQ=255)

                                                                                             atCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCa                                                >  1:404464/1‑150 (MQ=255)

                                                                                              tCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCt                                                                                                                     <  1:1382383/80‑1 (MQ=255)

                                                                                              tCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCag                                               >  2:912353/1‑150 (MQ=255)

                                                                                              tCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCag                                               >  2:409727/1‑150 (MQ=255)

                                                                                                catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAAc                                                                                         <  1:76482/106‑1 (MQ=255)

                                                                                                catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAAc                                                                                         >  2:76482/1‑106 (MQ=255)

                                                                                                catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCagaa                                             >  1:2286183/1‑150 (MQ=255)

                                                                                                catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCa                                                <  1:1631759/147‑1 (MQ=255)

                                                                                                catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCa                                                >  2:1631759/1‑147 (MQ=255)

                                                                                                 atcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAgg                                                                                   <  2:908418/111‑1 (MQ=255)

                                                                                                 atcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAgg                                                                                   >  1:908418/1‑111 (MQ=255)

                                                                                                 atcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCagaag                                            >  2:1825030/1‑150 (MQ=255)

                                                                                                  tcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCgccag                                                                                                >  1:588132/1‑97 (MQ=255)

                                                                                                     ttCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACggg                                                                                                             <  2:830791/81‑1 (MQ=255)

                                                                                                     ttCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACggg                                                                                                             >  1:830791/1‑81 (MQ=255)

                                                                                                     ttCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAg                                                                    >  1:1836175/1‑122 (MQ=255)

                                                                                                     ttCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAg                                                                    <  2:1836175/122‑1 (MQ=255)

                                                                                                     ttCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAg                                                            >  1:1864116/1‑130 (MQ=255)

                                                                                                      tCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGa                                       <  1:205528/150‑1 (MQ=255)

                                                                                                       cGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGGGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAg                                      <  2:2043601/150‑1 (MQ=255)

                                                                                                       cGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAg                                      <  1:156242/150‑1 (MQ=255)

                                                                                                       cGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAg                                      <  1:1503000/150‑1 (MQ=255)

                                                                                                       cGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAg                                      <  2:1568309/150‑1 (MQ=255)

                                                                                                             aGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATAc                                <  1:660959/150‑1 (MQ=255)

                                                                                                              tgtagaAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGTCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGCCAGCCACGCAGAAGATCCAGGCATAcc                               >  2:1150834/4‑150 (MQ=255)

                                                                                                              gCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATAcc                               <  2:1511829/150‑1 (MQ=255)

                                                                                                                  gaAACGGCAGCGTTCTTCAGCCATGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAg                           <  1:2089365/150‑1 (MQ=255)

                                                                                                                   aaaCGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCT‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAg                           <  2:2188945/149‑1 (MQ=255)

                                                                                                                     aCGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCCGCAACAACGTTCACCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGgtg                        >  1:286449/1‑150 (MQ=255)

                                                                                                                      cGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGcc                                                                                                                                      >  2:655008/1‑40 (MQ=255)

                                                                                                                      cGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGcc                                                                                                                                      <  1:655008/40‑1 (MQ=255)

                                                                                                                       ggCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCcac                                                                                            <  1:1218018/80‑1 (MQ=255)

                                                                                                                       ggCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCcac                                                                                            >  2:1218018/1‑80 (MQ=255)

                                                                                                                         cAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTAcc                    <  2:1905953/150‑1 (MQ=255)

                                                                                                                          aGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCg                   >  1:142830/1‑150 (MQ=255)

                                                                                                                           tgaacgctcAGCC‑GGCGCTGGCGACCCCAGCGCCTGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGc                  >  2:1976105/8‑150 (MQ=37)

                                                                                                                           gCGTTCTTCAGCCAGGCGATGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGc                  <  1:1817228/150‑1 (MQ=255)

                                                                                                                             gTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCa                >  1:1783284/1‑150 (MQ=255)

                                                                                                                             gTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGAC‑GGCTGAGCCGAGAGGCTCAACAAGATCGTCTCTCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCa                >  1:1855359/1‑150 (MQ=255)

                                                                                                                                gctcAGCC‑GGCGCTGGCGACCCCAGCGCCTGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACgg             >  1:296251/3‑150 (MQ=38)

                                                                                                                                   cAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGgcgc          <  1:1630088/150‑1 (MQ=255)

                                                                                                                                     gCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGcgcgccgcgc                                                                                                   >  2:2186076/1‑59 (MQ=255)

                                                                                                                                     gCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGcgcgccgcgc                                                                                                   <  1:2186076/59‑1 (MQ=255)

                                                                                                                                      ccAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCa       >  1:809130/1‑150 (MQ=255)

                                                                                                                                       ccggCGCTGGCGACCCCAGCGCCTGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAgg                                     <  2:1875595/118‑1 (MQ=38)

                                                                                                                                       ccggCGCTGGCGACCCCAGCGCCTGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAgg                                     >  1:1875595/3‑120 (MQ=38)

                                                                                                                                       ccggCGCTGGCGACCCCAGCGCCTGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCa       >  2:1703363/3‑150 (MQ=255)

                                                                                                                                       cAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATAcccc                             >  1:94999/1‑127 (MQ=255)

                                                                                                                                       cAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATAcccc                             <  2:94999/127‑1 (MQ=255)

                                                                                                                                       cAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCAt      >  2:867913/1‑150 (MQ=255)

                                                                                                                                           cgcTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCAcc  <  1:475662/150‑1 (MQ=255)

GGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACC  >  NC_002947/3337680‑3337968

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`