BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,389,946	C→G	5.6%	A23G (GCC→GGC)	PP_1212 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,389,946	0	C	G	5.6%	292.8 / 5.1	89	A23G (GCC→GGC)	PP_1212	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (49/35); new base G (3/2); total (52/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.77e-01

CGCAGTGCTTGCTAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAGGCC  >  NC_002947/1389800‑1390094

cgcAGTGCTTGCTAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAgcgccg                                                                                                                                                    >  3:31886/1‑149 (MQ=255)

cgcAGTGCTTGCTAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAgcgc                                                                                                                                                      >  5:23395/1‑147 (MQ=255)

cgcAGTGCTTGCTAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAgcgc                                                                                                                                                      <  6:23395/147‑1 (MQ=255)

 gcAGTGCTTGCTAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAgcggcgc                                                                                                                                                   >  7:208486/1‑149 (MQ=255)

       cttgctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCt                                                                                                                                             <  4:31886/149‑1 (MQ=255)

          gctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGa                                                                                                                                           <  6:356035/148‑1 (MQ=255)

           ctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGGCGCGGCGCTGATc                                                                                                                                         <  8:208486/149‑1 (MQ=255)

                cATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGGCGCGGCGCTGATcttctg                                                                                                                                    <  1:27481/149‑6 (MQ=255)

                 atatatAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTAGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTg                                                                                                                                    <  1:300359/149‑1 (MQ=255)

                 atatatAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGc                                                                                                                                   >  3:151577/1‑149 (MQ=255)

                        aTCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCt                                                                                                                            >  2:261610/1‑149 (MQ=255)

                         tCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTg                                                                                                                           >  1:54548/1‑149 (MQ=255)

                          cccGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGt                                                                                                                          >  7:274516/1‑149 (MQ=255)

                           ccGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGCCACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGAGGAGCGATTGCCCGGCATGTc                                                                                                                         >  4:273790/1‑149 (MQ=255)

                           ccGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGt                                                                                                                          >  1:251892/1‑148 (MQ=255)

                           ccGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGt                                                                                                                          >  4:181621/1‑148 (MQ=255)

                           ccGCTTTTTGCTGCCCCCAGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGAGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATAAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGt                                                                                                                          >  4:273771/1‑148 (MQ=255)

                            cGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGt                                                                                                                          >  6:187011/1‑147 (MQ=255)

                            cGCTTTTTGCTGACCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCAAGCCCCTTTATGACTATCAATGTGCGTCTT‑GCGAGCAACGAATGGAAATGCTGCAGAAGATCAGCGCCGCGCCGCTCATCGATTGCACGGACTgacc                                                                                                                        >  8:156395/1‑146 (MQ=255)

                                ttttGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCGGATCGATTGCCCGGCCTGTCAGgcg                                                                                                                    >  5:325781/1‑149 (MQ=255)

                                     ctgcCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCgg                                                                                                                <  2:251892/148‑1 (MQ=255)

                                        ccccccGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGACAGgccgggcaa                                                                                                              <  1:258910/147‑8 (MQ=255)

                                               gCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGgc                                                                                                                     >  2:244703/1‑133 (MQ=255)

                                               gCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGgc                                                                                                                     <  1:244703/133‑1 (MQ=255)

                                                    cagcGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAgctgct                                                                                                >  1:103915/1‑149 (MQ=255)

                                                     agcGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCgg                                                                                                                               >  7:211100/1‑117 (MQ=255)

                                                     agcGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCgg                                                                                                                               <  8:211100/117‑1 (MQ=255)

                                                                     agagaCACGCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAACGATGCTGTCGGTGCCAGGCTTTc                                                                               >  2:90152/1‑149 (MQ=255)

                                                                     agagaCACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTc                                                                               <  6:118985/149‑1 (MQ=255)

                                                                     agagaCACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTc                                                                               <  3:181621/149‑1 (MQ=255)

                                                                       agaCACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCTCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGGTGTCGGTGCCAGGCTTTCGc                                                                             <  1:243575/149‑1 (MQ=255)

                                                                       agaCACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGTCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGc                                                                             <  1:90152/149‑1 (MQ=255)

                                                                       agaCACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAgg                                                                                    <  4:23029/142‑1 (MQ=255)

                                                                       agaCACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGc                                                                             >  1:198775/1‑149 (MQ=255)

                                                                       agaCACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTACTGTCGGTGCCAGGCTTTCGc                                                                             <  1:261610/149‑1 (MQ=255)

                                                                       agaCACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGATGCTGTCGGTGCCAGGcagctcg                                                                             >  3:23029/1‑143 (MQ=255)

                                                                           acCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCt                                                                           >  1:234016/1‑147 (MQ=255)

                                                                             ccATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCa                                                                                                              <  6:208299/110‑1 (MQ=255)

                                                                             ccATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCa                                                                                                              >  5:208299/1‑110 (MQ=255)

                                                                             ccATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGc                                                                       >  3:267433/1‑149 (MQ=255)

                                                                             ccATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGc                                                                       >  3:154934/1‑149 (MQ=255)

                                                                              cATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCg                                                                      >  5:12865/1‑149 (MQ=255)

                                                                                    ccTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACg                                                                >  5:109309/1‑149 (MQ=255)

                                                                                    ccTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACg                                                                >  2:251217/1‑149 (MQ=255)

                                                                                     cTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACgg                                                               >  1:351562/1‑149 (MQ=255)

                                                                                        tATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGt                                                              <  8:274516/147‑1 (MQ=255)

                                                                                        tATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGt                                                              <  1:335519/147‑1 (MQ=255)

                                                                                           gACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTa                                                         >  5:339040/1‑149 (MQ=255)

                                                                                           gACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTa                                                         >  5:210196/1‑149 (MQ=255)

                                                                                            aCTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTa                                                         <  6:12865/148‑1 (MQ=255)

                                                                                            aCTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTAc                                                        >  2:13610/1‑149 (MQ=255)

                                                                                             cTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTAc                                                        >  2:30296/1‑148 (MQ=255)

                                                                                              tATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCAGCAGAAGATCAGCGCCGCGCTGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACg                                                       >  1:263909/1‑148 (MQ=255)

                                                                                                  aaTGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAAcc                                                  >  3:375852/1‑149 (MQ=255)

                                                                                                  aaTGTGCCTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACAGTTGGTACGAAAcc                                                  <  8:184952/149‑1 (MQ=255)

                                                                                                      tgCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCTGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGAc                                               <  2:263909/148‑1 (MQ=255)

                                                                                                      tgCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACt                                              <  3:273790/149‑1 (MQ=255)

                                                                                                      tgCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACt                                              <  6:109309/149‑1 (MQ=255)

                                                                                                        cGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTc                                            <  1:251217/149‑1 (MQ=255)

                                                                                                         gTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTc                                            >  1:7138/1‑148 (MQ=255)

                                                                                                         gTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTc                                            >  2:26693/1‑148 (MQ=255)

                                                                                                          tCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCGGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTc                                            >  6:22244/1‑147 (MQ=255)

                                                                                                            cT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGa                                        >  4:58483/1‑149 (MQ=255)

                                                                                                               gCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCACCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGAcc                                      >  1:25611/1‑149 (MQ=255)

                                                                                                                cGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTAc                                                        <  5:359169/130‑1 (MQ=255)

                                                                                                                cGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCAGATGTCGGTGCCAGGCTTTCGCCTGAGCGCTAACGGTTGGTAc                                                        >  6:359169/1‑130 (MQ=255)

                                                                                                                       cGCATGGAAGTGCTGCAGAAGATCAGCACCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGaa                              <  2:25611/149‑1 (MQ=255)

                                                                                                                       cGCATGGAAGTGCTGCAGAAGAGCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTTGTACGAAACCGACTTCAAGACCGGGGCGaa                              <  3:273771/149‑1 (MQ=255)

                                                                                                                        gCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCTGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGaa                              >  6:293146/1‑148 (MQ=255)

                                                                                                                         cATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGaaa                             >  7:306612/1‑148 (MQ=255)

                                                                                                                               aGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCg                                                 >  3:321440/1‑122 (MQ=255)

                                                                                                                               aGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCg                                                 <  4:321440/122‑1 (MQ=255)

                                                                                                                                 tgctgcAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTg                    <  4:151577/149‑1 (MQ=255)

                                                                                                                                   ctgcAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGc                  <  3:58483/149‑1 (MQ=255)

                                                                                                                                    tgcAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCTCTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATc                      >  7:66946/1‑144 (MQ=255)

                                                                                                                                    tgcAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAACCCGACTTCAAGACCGGGGCGAAAAAGAATCTgg                   >  4:95624/1‑147 (MQ=255)

                                                                                                                                    tgcAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATc                      <  8:66946/144‑1 (MQ=255)

                                                                                                                                     gcAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCGGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAg                >  8:75557/1‑149 (MQ=255)

                                                                                                                                     gcAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTgg                   >  4:149738/1‑146 (MQ=255)

                                                                                                                                     gcAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTgg                   <  3:149738/146‑1 (MQ=255)

                                                                                                                                       agaagaTCAGCGCCGCGGCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAggc              <  2:54548/149‑1 (MQ=255)

                                                                                                                                       agaagaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAggc              <  6:257802/149‑1 (MQ=255)

                                                                                                                                        atcgaTCAGCGGCGCGGCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAggcg             >  5:215860/4‑149 (MQ=255)

                                                                                                                                         aagaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGaaaaa                           <  5:354675/134‑1 (MQ=255)

                                                                                                                                         aagaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGaaaaa                           >  6:354675/1‑134 (MQ=255)

                                                                                                                                         aagaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACAGGGGCGAAAAAGAATCTGGCAggcgg            >  7:349576/1‑149 (MQ=255)

                                                                                                                                           gaTCAGCGGCGCGGCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGgcggcg          >  1:361440/1‑149 (MQ=255)

                                                                                                                                           gaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAggcggc           >  2:227450/1‑148 (MQ=255)

                                                                                                                                           gaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGgcggcg          <  4:154934/149‑1 (MQ=255)

                                                                                                                                            aTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGa         <  6:325781/149‑1 (MQ=255)

                                                                                                                                            aTCAGAGGCGCGGCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGa         >  5:136570/1‑149 (MQ=255)

                                                                                                                                             tCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGa         >  3:16168/1‑148 (MQ=255)

                                                                                                                                               aGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGAc        >  2:264309/1‑147 (MQ=255)

                                                                                                                                                gcggcgcggcGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAg     >  8:257558/1‑149 (MQ=255)

                                                                                                                                                gcgccgcgccGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAg     <  6:339040/149‑1 (MQ=255)

                                                                                                                                                   ccgcgccgcTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAACCTGCTGTCGGGGCCAGGCTTGCGCCTGAGCGGTAACGGGTGGTACGAAACCGACTTCCAGCCCGGGGCGAAAAAGACTCTGGCAGGCGGCGTCAAGGcc  <  6:196173/149‑1 (MQ=255)

CGCAGTGCTTGCTAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAGGCC  >  NC_002947/1389800‑1390094

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