BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,519,981	G→T	5.7%	A52S (GCC→TCC)	yrfF →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,519,981	0	G	T	5.7%	1008.6 / 21.3	263	A52S (GCC→TCC)	yrfF	inner membrane protein
Reads supporting (aligned to +/- strand): ref base G (120/128); new base T (9/6); total (129/134)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 275 total aligned reads displayed.

ACCATTGTGATTTTTTTAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTATCACGCGTTACGGCATCT  >  CP009273/3519834‑3520128

aCCATTGTGATTTTTTTAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAgcg                                                                                                                                                     <  2:2760288/148‑1 (MQ=255)

   aTTGTGATTTTTTTAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCg                                                                                                                                                  <  1:1149549/148‑1 (MQ=255)

     tgtgATTTTTTTAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGtt                                                                                                                                                >  1:3032647/1‑148 (MQ=255)

     tgtgATTTTTTTAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGtt                                                                                                                                                >  1:1190734/1‑148 (MQ=255)

     tgtgATTTTTTTAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGtt                                                                                                                                                >  2:896433/1‑148 (MQ=255)

        gATTTTTTTAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCTTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGaa                                                                                                                                             <  1:1041144/148‑1 (MQ=255)

        gATTTTTTTAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGaa                                                                                                                                             >  1:1371368/1‑148 (MQ=255)

          tttttttAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCACACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGaaaa                                                                                                                                           >  2:1991602/1‑148 (MQ=255)

          tttttttAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGaaaa                                                                                                                                           <  2:778685/148‑1 (MQ=255)

          tttttttAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGaaaa                                                                                                                                           <  1:841483/148‑1 (MQ=255)

          tttttttAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGaaaa                                                                                                                                           <  2:2703054/148‑1 (MQ=255)

           ttttttAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCTCCGTTGAAAAt                                                                                                                                          >  2:2138419/1‑148 (MQ=255)

            tttttAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAAtt                                                                                                                                         <  2:953298/148‑1 (MQ=255)

             ttttAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTa                                                                                                                                        <  2:1340015/148‑1 (MQ=255)

              tttAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTAt                                                                                                                                       <  2:2854574/148‑1 (MQ=255)

                tAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCt                                                                                                                                     <  2:284086/148‑1 (MQ=255)

                tAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCt                                                                                                                                     <  1:3742937/148‑1 (MQ=255)

                  gctgctTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTg                                                                                                                                   <  2:3606577/148‑1 (MQ=255)

                  gctgctTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTg                                                                                                                                   >  2:2019890/1‑148 (MQ=255)

                   ctgctTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTga                                                                                                                                  >  2:2772400/1‑148 (MQ=255)

                    tgctgtgctgGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTgag                                                                                                                                 <  2:1060537/148‑1 (MQ=255)

                    tgctTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTgag                                                                                                                                 <  2:318730/148‑1 (MQ=255)

                    tgctTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTgag                                                                                                                                 <  2:1667850/148‑1 (MQ=255)

                    tgctTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTgag                                                                                                                                 <  1:3674854/148‑1 (MQ=255)

                      ctTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGagag                                                                                                                               >  1:3264501/1‑148 (MQ=255)

                      ctTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGagag                                                                                                                               >  1:3264640/1‑148 (MQ=255)

                       tttGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGc                                                                                                                              <  1:2523814/148‑1 (MQ=255)

                        ttGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCTCCGTTGAAAATTATCTTGAGAGcc                                                                                                                             <  2:1729500/148‑1 (MQ=255)

                         tGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCt                                                                                                                            >  2:53193/1‑148 (MQ=255)

                          gCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTg                                                                                                                           <  1:617057/148‑1 (MQ=255)

                              gCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGc                                                                                                                       >  2:1713405/1‑148 (MQ=255)

                               ccTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCa                                                                                                                      >  1:231935/1‑148 (MQ=255)

                                 tGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAgg                                                                                                                    <  2:2209789/148‑1 (MQ=255)

                                   cTCACTACTTGCGGGATGGCTGATAAAAGTGCGCTCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTa                                                                                                                  <  2:630709/148‑1 (MQ=255)

                                   cTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTa                                                                                                                  <  1:2969981/148‑1 (MQ=255)

                                   cTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTa                                                                                                                  >  1:2329752/1‑148 (MQ=255)

                                   cTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTa                                                                                                                  <  2:630717/148‑1 (MQ=255)

                                    tCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTAt                                                                                                                 <  1:2386126/148‑1 (MQ=255)

                                     cACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTAtt                                                                                                                >  1:637444/1‑148 (MQ=255)

                                          ctTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAgg                                                                                                           <  1:681411/148‑1 (MQ=255)

                                          ctTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAgg                                                                                                           <  1:1013385/148‑1 (MQ=255)

                                                 gATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTgg                                                                                                    >  1:3015856/1‑148 (MQ=255)

                                                   tGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGcc                                                                                                  <  1:3780744/148‑1 (MQ=255)

                                                     gCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCa                                                                                                <  2:1575609/148‑1 (MQ=255)

                                                     gCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCa                                                                                                >  2:2638569/1‑148 (MQ=255)

                                                         aTAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACgg                                                                                            >  1:308618/1‑148 (MQ=255)

                                                          tAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACggg                                                                                           <  2:3712448/148‑1 (MQ=255)

                                                           aaaaGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGa                                                                                          >  2:2446331/1‑148 (MQ=255)

                                                           aaaaGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGa                                                                                          <  2:3769221/148‑1 (MQ=255)

                                                              aGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGcc                                                                                       <  1:1215500/148‑1 (MQ=255)

                                                               gTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCa                                                                                      >  1:3328930/1‑148 (MQ=255)

                                                                tGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAg                                                                                     >  2:1218717/1‑148 (MQ=255)

                                                                tGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAg                                                                                     >  1:523651/1‑148 (MQ=255)

                                                                      ccAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGc                                                                               >  1:534071/1‑148 (MQ=255)

                                                                       cAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCa                                                                              <  1:1637715/148‑1 (MQ=255)

                                                                       cAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCa                                                                              <  1:2202841/148‑1 (MQ=255)

                                                                         gACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCATTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCAcc                                                                            >  1:1920526/1‑148 (MQ=255)

                                                                             gCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGAtc                                                                        <  2:605751/148‑1 (MQ=255)

                                                                              cGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGAtct                                                                       <  2:1190734/148‑1 (MQ=255)

                                                                               gTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGAtctc                                                                      <  1:1267359/148‑1 (MQ=255)

                                                                                 cAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGAtctctc                                                                    <  1:2416619/148‑1 (MQ=255)

                                                                                 cAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGAtctctc                                                                    <  2:727389/148‑1 (MQ=255)

                                                                                   gCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTg                                                                  >  1:893578/1‑148 (MQ=255)

                                                                                    cTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTgg                                                                 >  1:2099250/1‑148 (MQ=255)

                                                                                          tGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGa                                                           <  2:848444/148‑1 (MQ=255)

                                                                                           ggACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGaa                                                          <  1:908038/148‑1 (MQ=255)

                                                                                           ggACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGaa                                                          <  2:231935/148‑1 (MQ=255)

                                                                                           ggACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGaa                                                          <  2:232730/148‑1 (MQ=255)

                                                                                              ccAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGc                                                       <  2:421073/148‑1 (MQ=255)

                                                                                                 aCGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGa                                                    >  1:2799247/1‑148 (MQ=255)

                                                                                                 aCGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGACAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGa                                                    >  2:2067127/1‑148 (MQ=255)

                                                                                                          cgGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGcaacaa                                           <  1:968322/148‑1 (MQ=255)

                                                                                                            gATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCTCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACg                                         >  1:642463/1‑148 (MQ=255)

                                                                                                            gATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACg                                         >  2:1172089/1‑148 (MQ=255)

                                                                                                              tGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTc                                       <  2:3305210/148‑1 (MQ=255)

                                                                                                               gcgcAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCa                                      <  2:1371804/148‑1 (MQ=255)

                                                                                                                cgcAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCat                                     <  1:2772400/148‑1 (MQ=255)

                                                                                                                 gcAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCatg                                    <  1:2161868/148‑1 (MQ=255)

                                                                                                                 gcAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCatg                                    <  1:1855716/148‑1 (MQ=255)

                                                                                                                  cAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCatga                                   <  1:2138524/148‑1 (MQ=255)

                                                                                                                  cAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCatga                                   <  2:1427453/148‑1 (MQ=255)

                                                                                                                  cAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCatga                                   <  1:1289413/148‑1 (MQ=255)

                                                                                                                      cgcgTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCt                               >  1:2909469/1‑148 (MQ=255)

                                                                                                                      cgcgTAAACTCACACCTGAAGAACGTAGCGCCGTGGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCt                               >  1:40403/1‑148 (MQ=255)

                                                                                                                              cTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACg                       <  2:1391986/148‑1 (MQ=255)

                                                                                                                              cTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACg                       <  1:1694454/148‑1 (MQ=255)

                                                                                                                                cacaCCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCt                     <  1:1415532/148‑1 (MQ=255)

                                                                                                                                cacaCCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCt                     >  1:3347977/1‑148 (MQ=255)

                                                                                                                                  cacCTGAAGAACGTAGCGCCGTTGTAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTAt                   >  1:1004777/1‑148 (MQ=255)

                                                                                                                                  cacCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTAt                   >  2:205446/1‑148 (MQ=255)

                                                                                                                                   acCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTATc                  <  1:2019890/148‑1 (MQ=255)

                                                                                                                                    ccTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTATCa                 >  1:3679321/1‑148 (MQ=255)

                                                                                                                                    ccTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTATCa                 <  2:2893430/148‑1 (MQ=255)

                                                                                                                                       gaagaaCGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTATCAcgc              >  1:3062851/1‑148 (MQ=255)

                                                                                                                                       gaagaaCGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTATCAcgc              >  2:1887829/1‑148 (MQ=255)

                                                                                                                                         agaaCGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTATCACGCGt            >  2:65055/1‑148 (MQ=255)

                                                                                                                                           aaCGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTATCACGCGTTa          >  2:3781881/1‑148 (MQ=255)

                                                                                                                                                 gcgcCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTATCACGCGTTACGGCAt    >  2:3042664/1‑148 (MQ=255)

                                                                                                                                                 gcgcCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTATCACGCGTTACGGCAt    >  2:2014171/1‑148 (MQ=255)

                                                                                                                                                   gcCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTATCACGCGTTACGGCAtct  <  2:1949871/148‑1 (MQ=255)

ACCATTGTGATTTTTTTAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTATCACGCGTTACGGCATCT  >  CP009273/3519834‑3520128

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