BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,546,084	G→A	8.1%	H18H (CAC→CAT)	yjiC ←	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,546,084	0	G	A	8.1%	850.6 / 43.4	248	H18H (CAC→CAT)	yjiC	uncharacterized protein
Reads supporting (aligned to +/- strand): ref base G (111/117); new base A (10/10); total (121/127)
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.97e-01

Only 100 of 264 total aligned reads displayed.

CAAGTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGACATGCATGTTT  >  CP009273/4545936‑4546232

cAAGTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATgaaggaag                                                                                                                                                      <  1:1334187/149‑1 (MQ=255)

   gTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGa                                                                                                                                                   <  1:2360535/149‑1 (MQ=255)

   gTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGa                                                                                                                                                   >  1:532730/1‑149 (MQ=255)

        cTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATc                                                                                                                                              <  2:2141316/149‑1 (MQ=255)

           aTTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGt                                                                                                                                           >  1:4514150/1‑149 (MQ=255)

           aTTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGt                                                                                                                                           >  1:1732541/1‑149 (MQ=255)

           aTTTCGACAACGCTTTTTCCATCGTAATTAGGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGt                                                                                                                                           >  2:225672/1‑149 (MQ=255)

             ttCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGtaa                                                                                                                                         <  1:1656183/149‑1 (MQ=255)

                 aCAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTaataat                                                                                                                                     >  2:542365/1‑149 (MQ=255)

                   aaCGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAAttt                                                                                                                                   >  1:1614243/1‑149 (MQ=255)

                   aaCGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAAttt                                                                                                                                   <  2:4308039/149‑1 (MQ=255)

                      gCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGAc                                                                                                                                >  2:1067160/1‑149 (MQ=255)

                      gCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGAc                                                                                                                                >  2:3402158/1‑149 (MQ=255)

                      gCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGAc                                                                                                                                <  1:1101618/149‑1 (MQ=255)

                      gCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGAc                                                                                                                                >  2:2498122/1‑149 (MQ=255)

                       cTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGAct                                                                                                                               >  1:2050176/1‑149 (MQ=255)

                            tCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGaa                                                                                                                          >  2:265155/1‑149 (MQ=255)

                            tCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAATGAGCATCTGTAATAATTTGACTCTGaa                                                                                                                          >  2:4281879/1‑149 (MQ=255)

                             ccATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAg                                                                                                                         <  1:2400050/149‑1 (MQ=255)

                              cATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAATGAGCATCTGTAATAATTTGACTCTGAAgg                                                                                                                        >  2:1855065/1‑149 (MQ=255)

                               aTCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAggg                                                                                                                       <  1:1017266/149‑1 (MQ=255)

                               aTCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAggg                                                                                                                       >  2:2955379/1‑149 (MQ=255)

                                tCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGc                                                                                                                      <  1:3770832/149‑1 (MQ=255)

                                  gTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCAt                                                                                                                    <  1:2284371/149‑1 (MQ=255)

                                  gTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAATGAGCATCTGTAATAATTTGACTCTGAAGGGCAt                                                                                                                    >  1:3532724/1‑149 (MQ=255)

                                      ttATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCt                                                                                                                <  1:2554872/149‑1 (MQ=255)

                                      ttATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCt                                                                                                                >  1:2714169/1‑149 (MQ=255)

                                         tGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTaa                                                                                                             >  2:4406645/1‑149 (MQ=255)

                                         tGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTaa                                                                                                             >  2:926496/1‑149 (MQ=255)

                                         tGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTaa                                                                                                             <  1:1037974/149‑1 (MQ=255)

                                          ggTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAAt                                                                                                            <  1:3100378/149‑1 (MQ=255)

                                          ggTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAATGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAAt                                                                                                            <  2:3839342/149‑1 (MQ=255)

                                           gTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATg                                                                                                           <  2:324675/149‑1 (MQ=255)

                                           gTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATg                                                                                                           <  1:1456788/149‑1 (MQ=255)

                                            ttATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATgg                                                                                                          <  2:1561199/149‑1 (MQ=255)

                                                  ccACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCCTCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGa                                                                                                    >  2:3047872/1‑149 (MQ=255)

                                                  ccACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGa                                                                                                    <  1:1957064/149‑1 (MQ=255)

                                                  ccACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGa                                                                                                    <  1:4289852/149‑1 (MQ=255)

                                                   cACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGAc                                                                                                   >  2:3341814/1‑149 (MQ=255)

                                                    aCGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAATGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACa                                                                                                  <  2:3976248/149‑1 (MQ=255)

                                                     cGAGTAAGAGGCGGAAGATAGATGTTTTCAGAATCCATGCGTGCATATTTTAATTTAAGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTTAGCATATGTAATAATTTGACTCTGAAGTGCATAGCTTAATGGAATATACat                                                                                                 >  1:3990163/1‑149 (MQ=255)

                                                           agagGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGt                                                                                           >  2:1212889/1‑149 (MQ=255)

                                                           agagGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGt                                                                                           >  2:3549211/1‑149 (MQ=255)

                                                              ggCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTCTTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTAc                                                                                        >  2:4460684/1‑149 (MQ=255)

                                                              ggCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTAc                                                                                        >  2:601117/1‑149 (MQ=255)

                                                                  gAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCttt                                                                                    >  1:339234/1‑149 (MQ=255)

                                                                  gAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAATGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCttt                                                                                    >  2:2318728/1‑149 (MQ=255)

                                                                   aaGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCCTTGCTTAATGGCATAGACATAATAGTTACCtttt                                                                                   >  2:2058477/1‑149 (MQ=255)

                                                                    agatagatGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTa                                                                                  >  1:4122608/1‑149 (MQ=255)

                                                                        agatGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAg                                                                              >  2:4112178/1‑149 (MQ=255)

                                                                        agatGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAATGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAg                                                                              <  1:685607/149‑1 (MQ=255)

                                                                         gatGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGt                                                                             <  1:945821/149‑1 (MQ=255)

                                                                         gatGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGt                                                                             <  1:3282752/149‑1 (MQ=255)

                                                                         gatGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGt                                                                             <  1:1197176/149‑1 (MQ=255)

                                                                           tGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTAt                                                                           <  1:3467481/149‑1 (MQ=255)

                                                                            gTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTAtt                                                                          <  1:2255080/149‑1 (MQ=255)

                                                                             ttttCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTAttt                                                                         <  2:2662428/149‑1 (MQ=255)

                                                                              tttCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAATGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTc                                                                        <  2:3670521/149‑1 (MQ=255)

                                                                                   gAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGa                                                                   <  1:3454856/149‑1 (MQ=255)

                                                                                     aTCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGCAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATg                                                                 >  1:2626951/1‑149 (MQ=255)

                                                                                      tCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGt                                                                >  2:2753322/1‑149 (MQ=255)

                                                                                       ccATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGtt                                                               <  1:4442976/149‑1 (MQ=255)

                                                                                       ccATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGtt                                                               >  1:3171225/1‑149 (MQ=255)

                                                                                        cATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGttt                                                              >  2:713850/1‑149 (MQ=255)

                                                                                             gTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGaaa                                                         >  2:1372169/1‑149 (MQ=255)

                                                                                               gCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGaaaaa                                                       <  2:888987/149‑1 (MQ=255)

                                                                                               gCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGaaaaa                                                       <  2:4417633/149‑1 (MQ=255)

                                                                                               gCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGaaaaa                                                       <  1:1975405/149‑1 (MQ=255)

                                                                                                cATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAAt                                                      <  2:1614243/149‑1 (MQ=255)

                                                                                                      ttAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATaa                                                >  2:3105946/1‑149 (MQ=255)

                                                                                                        aaTTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAATGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAAt                                              <  2:2754985/149‑1 (MQ=255)

                                                                                                         aTTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAAtt                                             >  1:2737280/1‑149 (MQ=255)

                                                                                                         aTTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAAtt                                             >  1:395729/1‑149 (MQ=255)

                                                                                                          tttACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAAttt                                            >  2:371010/1‑149 (MQ=255)

                                                                                                               gtgtAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTg                                       >  2:660522/1‑149 (MQ=255)

                                                                                                                 gtAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTgta                                     <  2:3740253/149‑1 (MQ=255)

                                                                                                                 gtAAACTCACTGTCGACTATTGGGTGATGAAGGAAATGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTgta                                     <  1:3821546/149‑1 (MQ=255)

                                                                                                                  tAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTgtag                                    >  2:2920753/1‑149 (MQ=255)

                                                                                                                   aaaCTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTgtagt                                   >  1:1546345/1‑149 (MQ=255)

                                                                                                                   aaaCTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTgtagt                                   >  1:920759/1‑149 (MQ=255)

                                                                                                                       tCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAACAATAGATAAATTTTTCTGTAGTaaaa                               >  2:4355843/1‑149 (MQ=255)

                                                                                                                       tCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTaaaa                               >  2:577777/1‑149 (MQ=255)

                                                                                                                            gTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAgaga                          <  1:3731424/149‑1 (MQ=255)

                                                                                                                            gTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTATGTAGTAAAAAgaga                          >  2:3613404/1‑149 (MQ=255)

                                                                                                                            gTCGACTATTGGGTGATGAAGGAAATGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAgaga                          >  1:3008738/1‑149 (MQ=255)

                                                                                                                                 cTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTaa                     <  1:2498965/149‑1 (MQ=255)

                                                                                                                                     tGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACaa                 >  2:151426/1‑149 (MQ=255)

                                                                                                                                      gggTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACaaa                <  1:541376/149‑1 (MQ=255)

                                                                                                                                      gggTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACaaa                <  1:21125/149‑1 (MQ=255)

                                                                                                                                        gTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATg              <  2:2441766/149‑1 (MQ=255)

                                                                                                                                         tgatgaAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGa             <  1:490453/149‑1 (MQ=255)

                                                                                                                                           atgaAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGAca           >  2:2824251/1‑149 (MQ=255)

                                                                                                                                           atgaAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGAca           <  1:2867286/149‑1 (MQ=255)

                                                                                                                                            tgaAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGAcat          <  2:2649503/149‑1 (MQ=255)

                                                                                                                                             gaaggaaaTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGAcatg         >  2:174167/1‑149 (MQ=255)

                                                                                                                                              aaggaagTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGAcatgc        <  2:531963/149‑1 (MQ=255)

                                                                                                                                               aggaagTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGAcatgca       >  2:3153194/1‑149 (MQ=255)

                                                                                                                                                ggaagTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGAcatgcat      >  1:1871803/1‑149 (MQ=255)

                                                                                                                                                ggaagTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGAcatgcat      >  2:512118/1‑149 (MQ=255)

                                                                                                                                                    gTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGACATGCATGttt  >  1:4203397/1‑149 (MQ=255)

CAAGTTTTCTAATTTCGACAACGCTTTTTCCATCGTAATTATGGTTATTGCCACGAGTAAGAGGCGGAAGATAGATGTTTTCCGAATCCATGCGTGCATATTTTAATTTACGTGTAAACTCACTGTCGACTATTGGGTGATGAAGGAAGTGAGCATCTGTAATAATTTGACTCTGAAGGGCATTGCTTAATGGCATAGACATAATAGTTACCTTTTATGAGTATTTCACTGATGTTTAGAAAAATAGATAAATTTTTCTGTAGTAAAAAGAGAAGTAAACAAATGACATGCATGTTT  >  CP009273/4545936‑4546232

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: `‑`