BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,873,695	C→A	56.9%	R296L (CGG→CTG)	recF ←	gap repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,873,695	0	C	A	56.9%	15.1 / inf	323	R296L (CGG→CTG)	recF	gap repair protein
Reads supporting (aligned to +/- strand): ref base C (75/64); new base A (99/85); total (174/149)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

Only 100 of 340 total aligned reads displayed.

TCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTG  >  CP009273/3873546‑3873844

tCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACAcc                                                                                                                                                       <  2:2327189/150‑1 (MQ=255)

  gAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCAGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACAccgc                                                                                                                                                     <  2:5050218/150‑1 (MQ=255)

   aaTTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACAcagcc                                                                                                                                                    >  1:335137/1‑150 (MQ=255)

       ttCGTCCGACATGTCTATACCGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGccgccc                                                                                                                                                <  2:4014373/150‑1 (MQ=255)

        tCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGccgcccg                                                                                                                                               <  1:2064264/150‑1 (MQ=255)

        tCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGccgcccg                                                                                                                                               <  2:3494986/150‑1 (MQ=255)

          gTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCt                                                                                                                                             >  1:1578321/1‑150 (MQ=255)

            ccGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCttt                                                                                                                                           >  2:4142700/1‑150 (MQ=255)

             cGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTc                                                                                                                                          <  2:2578427/150‑1 (MQ=255)

             agACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTc                                                                                                                                          <  2:809085/149‑1 (MQ=255)

              gACATGTCTATAACGTGTTCAGCACTGATCGCTCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCa                                                                                                                                         >  1:203029/1‑150 (MQ=255)

              gACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCa                                                                                                                                         >  2:661709/1‑150 (MQ=255)

               aCATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCGATCAGGTAGAGACACCGCCGCCCGCTTTCAc                                                                                                                                        <  2:309426/150‑1 (MQ=255)

                cATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACg                                                                                                                                       >  2:496380/1‑150 (MQ=255)

                    tCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTg                                                                                                                                   >  1:128061/1‑150 (MQ=255)

                     cTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGa                                                                                                                                  >  1:1906536/1‑150 (MQ=255)

                     cTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGa                                                                                                                                  <  2:2782142/150‑1 (MQ=255)

                     cTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGa                                                                                                                                  >  2:894937/1‑150 (MQ=255)

                       ataACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAgg                                                                                                                                <  1:2157290/150‑1 (MQ=255)

                         aaCGTGTTCAGCACTGATCGCGCAGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGaa                                                                                                                              >  2:5245065/1‑150 (MQ=255)

                               ttCAGCACAGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTAAAGAGGCTGGCAAGCAGCCCGCGACGCACATCATCAAGCTCAGAGGCAAAAACAACAATAAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTcc                                                                                                                        >  2:1533265/1‑150 (MQ=255)

                                 cAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCtt                                                                                                                      <  1:1428248/150‑1 (MQ=255)

                                  aGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTg                                                                                                                     <  1:1140597/150‑1 (MQ=255)

                                      cTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGcc                                                                                                                 <  1:2631083/150‑1 (MQ=255)

                                       tGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCa                                                                                                                <  1:4400673/150‑1 (MQ=255)

                                       tGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCa                                                                                                                <  2:1068953/150‑1 (MQ=255)

                                       tGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCa                                                                                                                <  1:2224417/150‑1 (MQ=255)

                                        gATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAg                                                                                                               <  1:1051603/150‑1 (MQ=255)

                                        gATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAg                                                                                                               >  2:4248709/1‑150 (MQ=255)

                                         aTCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGa                                                                                                              >  1:311741/1‑150 (MQ=255)

                                           cgcgCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACg                                                                                                            >  1:4000829/1‑150 (MQ=255)

                                              gcTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTaa                                                                                                         >  2:3933002/1‑150 (MQ=255)

                                              gcTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCACATCATCAAGCAAAGAGGCAAAATCAACTATCAGGTATAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACCTaa                                                                                                         >  1:4647125/1‑150 (MQ=255)

                                                tGACAAAGTCCTGTGATTGCGTCGCTTTTAAGCGTCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCCGAGGCAAAATCATCAATCAGGTTTAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACggcagg                                                                                                       >  2:2268713/1‑145 (MQ=255)

                                                tGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAgg                                                                                                       >  1:3803097/1‑150 (MQ=255)

                                                 gACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGgc                                                                                                      >  1:880283/1‑150 (MQ=255)

                                                  aCAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGgcg                                                                                                     <  1:1490806/150‑1 (MQ=255)

                                                  aCAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGgcg                                                                                                     >  1:4149617/1‑150 (MQ=255)

                                                    aaaGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGca                                                                                                   >  1:2818128/1‑150 (MQ=255)

                                                      aGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGcaca                                                                                                 >  2:4368864/1‑150 (MQ=255)

                                                         ccTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCa                                                                                              >  2:2799026/1‑150 (MQ=255)

                                                         ccTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCa                                                                                              <  2:3618659/150‑1 (MQ=255)

                                                            gtgATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACa                                                                                           <  1:1268531/150‑1 (MQ=255)

                                                             tgATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAg                                                                                          <  2:5158504/150‑1 (MQ=255)

                                                             tgATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGGAAGGCGCACATCAACAg                                                                                          >  1:957697/1‑150 (MQ=255)

                                                               aTTGCGTCGCTTTTAATCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCt                                                                                        >  1:3939796/1‑150 (MQ=255)

                                                               aTTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCt                                                                                        >  2:1512520/1‑150 (MQ=255)

                                                               aTTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCt                                                                                        >  2:4656660/1‑150 (MQ=255)

                                                                    gTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGc                                                                                   <  1:3233499/150‑1 (MQ=255)

                                                                     tCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCt                                                                                  <  2:2676411/150‑1 (MQ=255)

                                                                      cGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTg                                                                                 >  1:73450/1‑150 (MQ=255)

                                                                        cTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGcc                                                                               <  1:180449/150‑1 (MQ=255)

                                                                        cTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGcc                                                                               <  2:3345133/150‑1 (MQ=255)

                                                                        cTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGcc                                                                               <  1:4732088/150‑1 (MQ=255)

                                                                          tttAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCa                                                                             <  1:3609990/150‑1 (MQ=255)

                                                                          tttAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCa                                                                             <  2:2082205/150‑1 (MQ=255)

                                                                          ttgaaGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCa                                                                             <  1:4359951/147‑1 (MQ=255)

                                                                           ttAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCAc                                                                            >  1:450307/1‑150 (MQ=255)

                                                                                  gCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATaa                                                                     >  1:3175830/1‑150 (MQ=255)

                                                                                   cTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAg                                                                    <  2:4929008/150‑1 (MQ=255)

                                                                                   cTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAg                                                                    >  2:1339430/1‑150 (MQ=255)

                                                                                   cTGGCAAGCAGCCCGAGAAGATCATCATCAAGCTCATAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCAGCTTTCACGGGTGAGGAACTCTACTTGCGCCAGACGTAAGGCGCACATCAACAGATTAAGCTGCCCACGCGATAAg                                                                    <  1:1169029/150‑1 (MQ=255)

                                                                                     ggCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGt                                                                  <  2:1875008/150‑1 (MQ=255)

                                                                                       cAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTAt                                                                <  1:4656660/150‑1 (MQ=255)

                                                                                        aaGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATc                                                               >  2:2847085/1‑150 (MQ=255)

                                                                                        aaGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATc                                                               >  2:3715890/1‑150 (MQ=255)

                                                                                        aaGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATc                                                               <  2:3821605/150‑1 (MQ=255)

                                                                                              cccGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCAc                                                         <  1:3476102/150‑1 (MQ=255)

                                                                                              cccGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCAc                                                         <  1:2649245/150‑1 (MQ=255)

                                                                                               ccGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCAcc                                                        <  2:4640777/150‑1 (MQ=255)

                                                                                                  cgACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGgc                                                     >  1:1658639/1‑150 (MQ=255)

                                                                                                    aCGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGgcgc                                                   >  1:4817296/1‑150 (MQ=255)

                                                                                                       cTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTCTCTTCCACCGGCGCAcc                                                >  2:3658755/1‑150 (MQ=255)

                                                                                                            catcaAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCgg                                           >  1:4591978/1‑150 (MQ=255)

                                                                                                            catcaAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCgg                                           >  1:4982475/1‑150 (MQ=255)

                                                                                                             atcaAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGgc                                          >  2:4338810/1‑150 (MQ=255)

                                                                                                             atcaAGCACAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGgc                                          >  2:3866361/1‑150 (MQ=255)

                                                                                                                 aGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACCTCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGa                                      <  1:826753/150‑1 (MQ=255)

                                                                                                                   cTCAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAAt                                    <  2:3207202/150‑1 (MQ=255)

                                                                                                                     cAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGc                                  >  1:4111421/1‑150 (MQ=255)

                                                                                                                     cAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGc                                  >  1:2434506/1‑150 (MQ=255)

                                                                                                                     cAGAGGCAAAATCATCTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGc                                  >  2:2659184/1‑150 (MQ=255)

                                                                                                                       gagGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTCTCTTCCACCGGCGCACCGTCGGCGCGAATGCGt                                >  2:1760691/1‑150 (MQ=255)

                                                                                                                       gagGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGt                                >  2:3811981/1‑150 (MQ=255)

                                                                                                                         ggCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTaa                              <  2:1393687/150‑1 (MQ=255)

                                                                                                                               atcatcTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGc                        >  1:5178234/1‑150 (MQ=255)

                                                                                                                                tcatcTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCgat                       <  1:2871571/150‑3 (MQ=255)

                                                                                                                                tcatcTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCt                       >  1:3718141/1‑150 (MQ=255)

                                                                                                                                tcatcTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCGTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCt                       <  2:4268770/150‑1 (MQ=255)

                                                                                                                                 catcTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCtt                      <  1:469832/150‑1 (MQ=255)

                                                                                                                                    cTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTtgt                   >  1:5003936/1‑150 (MQ=255)

                                                                                                                                    cTATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTtgt                   <  2:2939949/150‑1 (MQ=255)

                                                                                                                                     tATCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTtgtg                  >  1:2724740/1‑150 (MQ=255)

                                                                                                                                      aTCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGc                 <  2:1927857/150‑1 (MQ=255)

                                                                                                                                      aTCAGGTAGAGACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGc                 >  1:3048257/1‑150 (MQ=255)

                                                                                                                                             agagACACAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGt          >  1:484017/1‑150 (MQ=255)

                                                                                                                                              gagaCACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGACTGCGTAAGTCCGCTTTGTGCGGGCCGTg         >  2:3054571/1‑150 (MQ=255)

                                                                                                                                                gaCACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTgcg       <  1:5295222/150‑1 (MQ=255)

                                                                                                                                                 acacAGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTgcgc      >  1:4263442/1‑150 (MQ=255)

                                                                                                                                                     cgccgcCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGgtg  >  1:5037838/1‑150 (MQ=255)

TCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGCAGCCCGCGACGCTCATCATCAAGCTCAGAGGCAAAATCATCTATCAGGTAGAGACACCGCCGCCCGCTTTCACGGGTGAGGAACTCTCCTTGCGCCAGACGTAAGGCGCACATCAACAGCTTAAGCTGCCCACGCGATAAGGTATCTTCCACCGGCGCACCGTCGGCGCGAATGCGTAAGTCCGCTTTGTGCGGGCCGTGCGCGGTG  >  CP009273/3873546‑3873844

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