BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,027,520	G→T	7.9%	C339F (TGT→TTT)	trkH →	potassium transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,027,520	0	G	T	7.9%	1116.5 / 16.9	348	C339F (TGT→TTT)	trkH	potassium transporter
Reads supporting (aligned to +/- strand): ref base G (165/144); new base T (35/4); total (200/148)
Fisher's exact test for biased strand distribution p-value = 8.45e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 356 total aligned reads displayed.

CCTCGTACTGTGGTTTCATAATGTCTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATCGCGCAC  >  CP009273/4027371‑4027669

ccTCGTACTGTGGTTTCATAATGTCTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATg                                                                                                                                                       <  2:1364383/150‑1 (MQ=255)

 cTCGTACTGTGGTTTCATAATGTCTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGt                                                                                                                                                      <  1:2428121/150‑1 (MQ=255)

  tCGTACTGTGGTTTCATAATGTCTACAGTTCGGCGCTGATGACAAATAACCAGGCGTTTTTCCAGGTGGTGTCGAAGGCGAAAACCGCCGGGATTAAAAATGACAACAAAGCCCACAGGCCGCAC‑ATTTTTCAGGTACTGCTTTTAtttt                                                                                                                                                     >  1:1687175/1‑150 (MQ=255)

    gTACTGTGGTTTCATAATGTCTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCa                                                                                                                                                   <  2:2166446/150‑1 (MQ=255)

      aCTGTGGTTTCATAATGTCTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGc                                                                                                                                                 <  2:5147207/150‑1 (MQ=255)

       cTGTGGTTTCATAATGTCTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGcc                                                                                                                                                >  2:2888917/1‑150 (MQ=255)

        tgtgGTTTCATAATGTCTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCt                                                                                                                                               >  2:2238202/1‑150 (MQ=255)

         gtgGTTTCATAATGTCTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTTCCGGGTCTGCTTTTTTTTTGAGccgt                                                                                                                                              >  1:2252401/1‑148 (MQ=255)

         gtgGTTTCATAATGTCTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCtt                                                                                                                                              >  2:2749884/1‑150 (MQ=255)

         gtgGTTTCATAATGTCTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCtt                                                                                                                                              >  2:1065596/1‑150 (MQ=255)

          tgGTTTCATAATGTCTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCttt                                                                                                                                             <  2:2119976/150‑1 (MQ=255)

            gTTTCATAATGTCTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTTCCGGTACTGCTTTTATTTTCAGCCTTTAt                                                                                                                                           >  1:430184/1‑150 (MQ=255)

            gTTTCATAATGTCTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTAt                                                                                                                                           <  1:4329792/150‑1 (MQ=255)

            gTTTCATAATGTCTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTAt                                                                                                                                           <  2:2316729/150‑1 (MQ=255)

              ttCATAATGTCTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATcg                                                                                                                                         <  1:3693065/150‑1 (MQ=255)

                cATAATGTCTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATcggc                                                                                                                                       >  2:1547839/1‑150 (MQ=255)

                    aTGTCTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGtt                                                                                                                                   <  1:1001089/150‑1 (MQ=255)

                    aTGTCTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGtt                                                                                                                                   <  2:4087999/150‑1 (MQ=255)

                     tGTCTACTGTTCGGCGCTTATGACAATTAACAAGTGTTTTTTCCAGGTGGTGTCTTTGTCGAAAACCGCCGGGTTTACATCTGACAGCATTTCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTg                                                                                                                                  <  1:1357267/150‑1 (MQ=255)

                      gTCTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTgc                                                                                                                                 <  1:4172934/150‑1 (MQ=255)

                        cTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTgcgc                                                                                                                               >  2:4727480/1‑150 (MQ=255)

                          aCAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCCCCCGGTTTAAAACAGACAGCATTTCCCCCAGGCCGCGC‑TTTTTTGCGGGACGGCTTTTTTTTTAAACCCTTATAGGCGGTTGGGCCg                                                                                                                             >  1:1480362/1‑150 (MQ=255)

                              ttCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTc                                                                                                                         <  1:474432/150‑1 (MQ=255)

                                 ggCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCCCCGGGTTTAAAACAGACAGCAATTCCCGCTGGCCGCTCTTTTTTGCCGGT‑CAGCATTTTTTTTAAGCCATTATAAGCGGGTGCGCCCGGGCaaa                                                                                                                      >  1:4533802/1‑149 (MQ=255)

                                  gcgcTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACg                                                                                                                     <  2:702155/150‑1 (MQ=255)

                                   cgcTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTAAAACTGACAGCATTTCCCGCTGGCCGCTC‑TTTTTTCCGGGACTGCTTTTATGTTAAGCCTTTATAGGCGGGTTCGCCGGGTCaaacg                                                                                                                    >  1:764443/1‑147 (MQ=255)

                                    gcTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACggg                                                                                                                   <  1:3868205/150‑1 (MQ=255)

                                      tgatgaCAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCg                                                                                                                 >  2:1071939/1‑150 (MQ=255)

                                      tgatgaCAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCAATGCCCGCTGGCCGCTC‑TTTTTTCCCGTACAGCTTTTAAGTTAAGCCTTTATAGGCGGTTGCGCCGGGTCAACGGGCg                                                                                                                 >  2:4925654/1‑150 (MQ=255)

                                        atgaCAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGt                                                                                                               <  1:82783/150‑1 (MQ=255)

                                            cAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACACCTGACAGCAGTGCCCGCTGGCCGCTG‑TGTTTGCCGGGCCGGCGGGTAGGTTCACCCGTGATCGGCGGGTGCGGCGGGTCAACGGGCGGGGGcc                                                                                                           >  2:3736393/1‑150 (MQ=255)

                                            cAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGcc                                                                                                           <  1:4782961/150‑1 (MQ=255)

                                            cAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCCCCGGGGTTACAAAAGACAGCAATGCCCGCTGGCCGCGCATTTTTGCCGG‑ACAGCTTTTATTTTAAGCCATTATAGGCGGGTGCGCCCGGGCAAACCGCGGGGGcc                                                                                                           >  1:718149/1‑150 (MQ=255)

                                              aTTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTg                                                                                                         >  2:3890704/1‑150 (MQ=255)

                                              aTTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAAATGACAGCATTGCCCGCTGGCCGCGC‑TTTTTTCCGGGACATCTTTTATTTTAAGCCTTTATAGGCGGTTGCGCCGGGTCAAAGGGCGGTGGccag                                                                                                         >  1:5206170/1‑148 (MQ=255)

                                                tAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCGGaa                                                                                                       >  2:4671787/1‑150 (MQ=255)

                                                 aaCCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAg                                                                                                      <  2:1374917/150‑1 (MQ=255)

                                                 aaCCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAg                                                                                                      >  2:4313787/1‑150 (MQ=255)

                                                   ccAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGt                                                                                                    <  1:3802198/150‑1 (MQ=255)

                                                   ccAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGt                                                                                                    <  2:4565978/150‑1 (MQ=255)

                                                   ccAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGt                                                                                                    <  1:3033821/150‑1 (MQ=255)

                                                   ccAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCCCCGGGTTTACAACTGAAAGCATTTCCCGCTGGCCGCTC‑TTTTTTCCGGTAATTCTTTTATTTTAAGCCTTTATAGGCGGTTTCGCCGGGTCAACCGGCGGTGGCCAGAAGGt                                                                                                    >  1:5094989/1‑150 (MQ=255)

                                                    cAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTg                                                                                                   >  2:2445480/1‑150 (MQ=255)

                                                    cAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTg                                                                                                   >  2:177114/1‑150 (MQ=255)

                                                     aGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTTCCGGTAATGCTTTTATGTTCAGCCTTTAACGGCGGTTGCGCCGGGTCAAAGGGCGGTGGCCTGAAGGTGa                                                                                                  >  2:4040736/1‑150 (MQ=255)

                                                     aGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGa                                                                                                  <  2:4527635/150‑1 (MQ=255)

                                                     aGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGa                                                                                                  <  2:436671/150‑1 (MQ=255)

                                                     aGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCAGAAGGTGa                                                                                                  >  2:2384773/1‑150 (MQ=255)

                                                     aGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑GTTTTTGCGGTACTGCGTTTATTTTCAGCCGTTATAGGCGGTTTCGCCCGGTCAACCGGCGGTGGCCGGGAGGTGa                                                                                                  >  1:3272093/1‑150 (MQ=255)

                                                       gCGTTTTTCCTGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATc                                                                                                <  1:2188474/150‑1 (MQ=255)

                                                       gCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATc                                                                                                <  2:324141/150‑1 (MQ=255)

                                                           ttttCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCa                                                                                            >  2:1170772/1‑150 (MQ=255)

                                                            tttCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTAAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCAt                                                                                           >  1:3342559/1‑150 (MQ=255)

                                                             ttCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACATCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATc                                                                                          <  1:1514276/150‑1 (MQ=255)

                                                             ttCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATc                                                                                          <  2:3300042/150‑1 (MQ=255)

                                                               cctggtggtGTCGATGGCGACAACCGCCGGGTTTACAACTGACATCATTGCCCGCTGGCCTCTC‑TTTTTGCCGTTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCct                                                                                        <  1:982886/147‑1 (MQ=255)

                                                               ccAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCct                                                                                        <  1:5127690/150‑1 (MQ=255)

                                                                 aGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCctgc                                                                                      >  2:2650514/1‑150 (MQ=255)

                                                                  ggtggtGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCctgct                                                                                     <  1:3786350/150‑1 (MQ=255)

                                                                           gATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTaa                                                                            >  2:2242692/1‑150 (MQ=255)

                                                                                  aCAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAgggg                                                                     >  2:2819471/1‑150 (MQ=255)

                                                                                   cAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGa                                                                    <  2:314808/150‑1 (MQ=255)

                                                                                    aaCCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTTCCGGGACTGCTTTTATGTTCAGCCATTATCCGCGGTTGCGCCGCGTCAACGGGCGGGGGCCAGAAGGTGAGCCCCAACCTGCAGCTGTTTTAGCACGGGaa                                                                   >  1:3608650/1‑150 (MQ=255)

                                                                                    aaCCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGaa                                                                   >  2:3482382/1‑150 (MQ=255)

                                                                                     aCCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATAGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAAc                                                                  >  1:2737633/1‑150 (MQ=255)

                                                                                      ccgccgGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAAcc                                                                 <  2:1831797/150‑1 (MQ=255)

                                                                                           gggTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAg                                                            >  2:617369/1‑150 (MQ=255)

                                                                                             gTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCGGCTGTTTAAGCAGGGGAACCGTGAGCt                                                          >  1:2020953/1‑150 (MQ=255)

                                                                                               ttACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGa                                                        <  2:4506889/150‑1 (MQ=255)

                                                                                                tACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTTCCGGTACTGGTTTTATGTTCAGCCTTTATAGGCGGTGGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCGGCTGCTGTTTAAGCAGGGGGACCGTGAGCGGaa                                                       >  1:10650/1‑150 (MQ=255)

                                                                                                  cAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCGGAAGGTGATCCGCATCCTGCGGCTGTTTAAGCAGGGGAACCGTGAGCGGAAAc                                                     >  1:1789036/1‑150 (MQ=255)

                                                                                                     cTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGGTTTTATGTTAAGCCTTTATCGGCGGGTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGGGAGCTGAAACGAc                                                  >  1:1395279/1‑150 (MQ=255)

                                                                                                     cTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGAc                                                  <  2:3285675/150‑1 (MQ=255)

                                                                                                      tGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACt                                                 <  1:3482938/150‑1 (MQ=255)

                                                                                                       gACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTg                                                <  1:1241495/150‑1 (MQ=255)

                                                                                                            cATTGCCCGCTGGCCGCTG‑TTTTTTCCGGTACTGCTTTTATGTTAAGCCGTTATAGGCGGTTTGGCCGGGTCAACCGGCGGTGGCCTGAAGGTGATCCGCATCCAGCGGCTGTTTAATCAGGGGAACCGTGAGCTGAAACGACTGGTGCa                                           >  1:922591/1‑150 (MQ=255)

                                                                                                             aTTGCCCGCTGGCCGCTC‑TTTTTTCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATACGCAACCTGCTGCTGTTTAATCAGGGGAACCCTGAGCTGAAACGACAGGTGCAt                                          >  1:537913/1‑150 (MQ=255)

                                                                                                               tGCCCGGTGGCCGCTA‑TTTTTGCGGGTACTGATTTTTTTTTCCGCCTTTATCGGCGTTTGCGCCGGGTCAACGGGCGTTGGCCTGATGGTTATCCGCATCCTGCTGCTGTTTAAGCTGGGTAACCGTTAGCTGAAACGACTGGTGCATcc                                        <  2:1993301/150‑1 (MQ=255)

                                                                                                                gCCCGCTGGCCGCTC‑TTTTTTCCGGTACAGCTTTTATTTTAAGCCTTTATAGGCGGTTTCGCCGGGTCAACGGGCGGTGGCCAGAATGTGATACGCATCCAGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAAAGACAGGTGCATccc                                       >  1:116968/1‑149 (MQ=255)

                                                                                                                   cgcgGGCCGCTC‑TTTTTGCCGGTCCTGCTTTTATGTTCAGCCTTTTTCGGCGGTTGCGCGGGGACACCGGGCGGTGGCCTGAAGGTTATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAAc                                    <  2:307953/150‑1 (MQ=255)

                                                                                                                     cTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGc                                  >  2:4435174/1‑150 (MQ=255)

                                                                                                                      tGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGcc                                 >  2:829796/1‑150 (MQ=255)

                                                                                                                      tGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTAAGCCTTTATCCGCGGTAGCGCCGGGTCAACGGGCGGTTGCCTGAAGGTGATCCGCAAACTGCTGCAGTTTAAGCAGGGGAACCGTGAGATGAAACGAATGGTGCATCCGAACGcc                                 >  1:4982251/1‑150 (MQ=255)

                                                                                                                           gCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCta                            >  1:2851662/1‑150 (MQ=255)

                                                                                                                              c‑tTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAg                         >  1:2925657/1‑150 (MQ=255)

                                                                                                                              c‑tTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGACGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAg                         >  2:1441025/1‑150 (MQ=255)

                                                                                                                                   ttGCCGGTCCGGCTTTTATGTTCCGCCTTTTTGGGCGGTTGCGCGGGGTCAAGGGGGGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTTTTTAAGCTGGGGAACCGTGTGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCAtt                     <  1:2992802/150‑1 (MQ=255)

                                                                                                                                   ttGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCGGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGTAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCAtt                     <  2:882544/150‑1 (MQ=255)

                                                                                                                                   ttGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCAtt                     >  1:3960012/1‑150 (MQ=255)

                                                                                                                                    tGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTa                    >  1:3200169/1‑150 (MQ=255)

                                                                                                                                    tGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTa                    >  2:3610270/1‑150 (MQ=255)

                                                                                                                                          tACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTgg              >  1:4654412/1‑150 (MQ=255)

                                                                                                                                             tGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGa           >  1:2562386/1‑150 (MQ=255)

                                                                                                                                             tGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGa           >  1:2097562/1‑150 (MQ=255)

                                                                                                                                             tGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGa           >  2:5048059/1‑150 (MQ=255)

                                                                                                                                               cTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAAGGGGCGGTGGCCTGAAGGTGAACCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAAt         <  2:1276995/150‑1 (MQ=255)

                                                                                                                                                ttttATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATc        >  2:4578672/1‑150 (MQ=255)

                                                                                                                                                 tttATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATcg       <  2:764443/150‑1 (MQ=255)

                                                                                                                                                   tATGTTCAGCCTTTATGGGCGGTTGCGCGGGGTCACGGGGCGGTGGCCTGAAGGTTATCCGCATCCTGCTGCTGTTTAAGCTGGGGAACCGTGAGCTGACACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATcgcg     <  1:3418511/150‑1 (MQ=255)

                                                                                                                                                   tATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATcgcg     >  1:2600613/1‑150 (MQ=255)

CCTCGTACTGTGGTTTCATAATGTCTACAGTTCGGCGCTGATGACAATTAACCAGGCGTTTTTCCAGGTGGTGTCGATGGCGACAACCGCCGGGTTTACAACTGACAGCATTGCCCGCTGGCCGCTC‑TTTTTGCCGGTACTGCTTTTATGTTCAGCCTTTATCGGCGGTTGCGCCGGGTCAACGGGCGGTGGCCTGAAGGTGATCCGCATCCTGCTGCTGTTTAAGCAGGGGAACCGTGAGCTGAAACGACTGGTGCATCCGAACGCCGTCTATAGCATTAAGCTGGGGAATCGCGCAC  >  CP009273/4027371‑4027669

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