BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,519,981	G→T	18.0%	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	18.0%	544.4 / 73.3	189	A52S (GCC→TCC)	yrfF	inner membrane protein
Reads supporting (aligned to +/- strand): ref base G (74/81); new base T (13/21); total (87/102)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 199 total aligned reads displayed.

ACCATTGTGATTTTTTTAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTATCACGCGTTACGGCATCT  >  CP009273/3519834‑3520128

aCCATTGTGATTTTTTTAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAgcg                                                                                                                                                     >  2:2069292/1‑148 (MQ=255)

 ccATTGTGATTTTTTTAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAgcgc                                                                                                                                                    <  1:2930254/148‑1 (MQ=255)

 ccATTGTGATTTTTTTAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAgcgc                                                                                                                                                    >  2:266474/1‑148 (MQ=255)

  cATTGTGATTTTTTTAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGcc                                                                                                                                                   >  2:1869833/1‑148 (MQ=255)

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

   aTTGTGATTTTTTTAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCg                                                                                                                                                  >  1:1550866/1‑148 (MQ=255)

    ttGTGATTTTTTTAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGt                                                                                                                                                 >  2:62001/1‑148 (MQ=255)

       tgATTTTTTTAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGa                                                                                                                                              >  2:3094797/1‑148 (MQ=255)

        gATTTTTTTAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCACACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGaa                                                                                                                                             >  2:1695222/1‑148 (MQ=255)

         aTTTTTTTAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGaaa                                                                                                                                            >  1:569816/1‑148 (MQ=255)

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

           ttttttAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAAt                                                                                                                                          <  1:2361379/148‑1 (MQ=255)

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

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

              tttAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCTCCGTTGAAAATTAt                                                                                                                                       >  2:1452162/1‑148 (MQ=255)

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

               ttAGCTGCTGTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCTCCGTTGAAAATTATc                                                                                                                                      <  2:979452/148‑1 (MQ=255)

                tAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCTCCGTTGAAAATTATCt                                                                                                                                     <  1:3270085/148‑1 (MQ=255)

                tAGCTGCTTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCTCCGTTGAAAATTATCt                                                                                                                                     <  1:1926410/148‑1 (MQ=255)

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

                     gctTTGCTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTgaga                                                                                                                                >  1:619609/1‑148 (MQ=255)

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

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

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

                           cTGGCCTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGa                                                                                                                          <  2:343873/148‑1 (MQ=255)

                               ccTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCTCCGTTGAAAATTATCTTGAGAGCCTGACTCa                                                                                                                      <  2:1020250/148‑1 (MQ=255)

                               ccTGCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCa                                                                                                                      <  2:823973/148‑1 (MQ=255)

                                cTGCTCACTACTTGCGGGATGGTTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAg                                                                                                                     <  2:706441/148‑1 (MQ=255)

                                  gCTCACTACTTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGt                                                                                                                   <  2:1818053/148‑1 (MQ=255)

                                      actactTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTa                                                                                                               <  2:1157498/148‑1 (MQ=255)

                                      actacgTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTa                                                                                                               <  2:1846102/148‑1 (MQ=255)

                                         actTGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAg                                                                                                            <  2:2388765/148‑1 (MQ=255)

                                           ttGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGt                                                                                                          >  2:2566325/1‑148 (MQ=255)

                                            tGCGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTg                                                                                                         >  2:1913234/1‑148 (MQ=255)

                                             gCGGGTTGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGc                                                                                                        <  2:2986573/148‑1 (MQ=255)

                                              cGGGATGGCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCTCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGcc                                                                                                       <  1:1596658/148‑1 (MQ=255)

                                                    ggCTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGccc                                                                                                 >  2:2098375/1‑148 (MQ=255)

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

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

                                                      cTGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCaa                                                                                               <  1:602495/148‑1 (MQ=255)

                                                       tGATAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCTCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAAc                                                                                              <  1:3019945/148‑1 (MQ=255)

                                                          tAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACggg                                                                                           >  2:2090626/1‑148 (MQ=255)

                                                          tAAAAGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACggg                                                                                           >  1:124032/1‑148 (MQ=255)

                                                             aaGTGCGATCCAGTCGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAACCTCACAACTGAAGAAGGTAGCTCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGc                                                                                        <  1:1015204/148‑1 (MQ=255)

                                                             aaGTGCGATCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGc                                                                                        <  1:2693288/148‑1 (MQ=255)

                                                                    aTCCAGACGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGcgcg                                                                                 <  1:3194578/148‑1 (MQ=255)

                                                                          aCGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCg                                                                           >  2:695332/1‑148 (MQ=255)

                                                                          aCGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCg                                                                           >  2:1359070/1‑148 (MQ=255)

                                                                           cGGCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGGAAATTATCTTGAGAGCCAGACGCAGGGATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGa                                                                          >  1:1457402/1‑148 (MQ=255)

                                                                            ggCGTCAGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCTCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGAt                                                                         <  2:435941/148‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                  aGCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCTCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATctctct                                                                   <  2:481147/148‑1 (MQ=255)

                                                                                   gCTGCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTg                                                                  <  1:961665/148‑1 (MQ=255)

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

                                                                                      gCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGATAGCCAGACGCACGTATTACAGGTGCCAGGCCCAACTGGAGCCAGCGCGGCACCGAACTCTCTGgcg                                                               >  1:716618/1‑148 (MQ=255)

                                                                                      gCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGgcg                                                               >  1:716927/1‑148 (MQ=255)

                                                                                      gCCCTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGgcg                                                               <  2:2310117/148‑1 (MQ=255)

                                                                                       cccTGGACCAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGgcgc                                                              >  2:1275536/1‑148 (MQ=255)

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

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

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

                                                                                               cAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGcc                                                      >  1:2298426/1‑148 (MQ=255)

                                                                                               cAACGCCTTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGcc                                                      >  1:3184610/1‑148 (MQ=255)

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

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

                                                                                                   gCCGTCGCGGATGCGCAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGaaa                                                  <  2:1544993/148‑1 (MQ=255)

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

                                                                                                               gcgcAAACGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCa                                      >  2:1435690/1‑148 (MQ=255)

                                                                                                                     aCGCGTAAACTCACACCTGAAGAACGTAGCTCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGc                                >  1:1314338/1‑148 (MQ=255)

                                                                                                                     aCGCGTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGc                                >  2:3264438/1‑148 (MQ=255)

                                                                                                                      cgcgTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCt                               <  1:1753147/148‑1 (MQ=255)

                                                                                                                        cgTAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGa                             <  1:3060044/148‑1 (MQ=255)

                                                                                                                          tAAACTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGaca                           >  2:2123052/1‑148 (MQ=255)

                                                                                                                            aaCTCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGacaca                         <  1:2017013/148‑1 (MQ=255)

                                                                                                                               tCACACCTGAAGAACGTAGCTCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGc                      <  2:3224784/148‑1 (MQ=255)

                                                                                                                               tCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGc                      >  1:85358/1‑148 (MQ=255)

                                                                                                                               tCACACCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGc                      >  1:2098857/1‑148 (MQ=255)

                                                                                                                                cacaCCTGAAGAACGTAGCTCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCt                     <  1:3118073/148‑1 (MQ=255)

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

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

                                                                                                                                 acacCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTa                    >  1:996003/1‑148 (MQ=255)

                                                                                                                                  cacCTGAAGAACGTAGCTCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTAt                   <  1:2033156/148‑1 (MQ=255)

                                                                                                                                  cacCTGAAGAACGTAGCTCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTAt                   >  2:2898468/1‑148 (MQ=255)

                                                                                                                                  cacCTGAAGAACGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTAt                   <  2:586416/148‑1 (MQ=255)

                                                                                                                                   acCTGAAGAACGTTGCGCCGTTGTAAATTATCTTGAGAGCCTGACGCTGGTATTACAGGTGCCTGGCCGAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTATc                  <  1:1488792/148‑1 (MQ=255)

                                                                                                                                          gaaCGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTATCACGCGtt           >  1:496390/1‑148 (MQ=255)

                                                                                                                                          gaaCGGAGCTCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTATCACGCGtt           <  1:1451885/148‑1 (MQ=255)

                                                                                                                                           aaCGTAGCTCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTATCACGCGTTa          >  1:230537/1‑148 (MQ=255)

                                                                                                                                           aaCGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTATCACGCGTTa          <  1:1367893/148‑1 (MQ=255)

                                                                                                                                           aaCGTAGCGCCATTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTATCACGCGTTa          >  2:2032670/1‑148 (MQ=255)

                                                                                                                                             cGTAGCGCCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTATCACGCGTTACg        >  2:1103705/1‑148 (MQ=255)

                                                                                                                                                 gcgcCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTATCACGCGTTACGGCAt    <  1:1275536/148‑1 (MQ=255)

                                                                                                                                                   gcCGTTGAAAATTATCTTGAGAGCCTGACGCAGGTATTACAGGTGCCTGGCCCAACGGGAGCCAGCGCGGCACCGATCTCTCTGGCGCTGAATGCCGAAAGCAACAACGTCATGATGCTGACACACGCTATCACGCGTTACGGCAtat  <  2:3010007/148‑3 (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: `‑`