BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	658,096	C→A	70.7%	M1M (ATG→ATT) †	ybeD ←	UPF0250 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	658,096	0	C	A	70.7%	392.8 / inf	350	M1M (ATG→ATT) †	ybeD	UPF0250 family protein
Reads supporting (aligned to +/- strand): ref base C (48/55); new base A (91/156); total (139/211)
Fisher's exact test for biased strand distribution p-value = 9.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 352 total aligned reads displayed.

GCTTGGTTTTACCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAAT  >  CP009273/657949‑658245

gctTGGTTTTACCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAGAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAAt                                                                                                                                                     <  1:2944669/150‑1 (MQ=255)

gctTGGTTTTACCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAAt                                                                                                                                                     <  2:316654/150‑1 (MQ=255)

gctTGGTTTTACCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAAt                                                                                                                                                     >  2:4171939/1‑150 (MQ=255)

    ggTTTTACCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCa                                                                                                                                                 >  1:59994/1‑150 (MQ=255)

     gTTTTACCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAg                                                                                                                                                >  2:180111/1‑150 (MQ=255)

     gTTTTACCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAg                                                                                                                                                <  2:3658212/150‑1 (MQ=255)

      ttttACCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGc                                                                                                                                               <  1:3804061/150‑1 (MQ=255)

      ttttACCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGc                                                                                                                                               <  2:152416/150‑1 (MQ=255)

         tACCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTcc                                                                                                                                            <  1:319658/150‑1 (MQ=255)

            cGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGc                                                                                                                                         >  1:5310770/1‑150 (MQ=255)

             gTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGcg                                                                                                                                        >  2:1299537/1‑150 (MQ=255)

              ttGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGcgt                                                                                                                                       >  1:4868002/1‑150 (MQ=255)

                ggggTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGcgtaa                                                                                                                                     <  1:4442994/150‑1 (MQ=255)

                  ggTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGcgtaacg                                                                                                                                   >  1:3364195/1‑150 (MQ=255)

                      tAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAAt                                                                                                                               >  2:4269450/1‑150 (MQ=255)

                       aGTCACGTGGAGGATGGCGCTGTACCACTTCATCCAGCTGATCAACCTGCTCAGGTAACGCCTGCCCCATACCTTTGTAAGTAAATGGAGTAGGGAATTGAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAAtt                                                                                                                              <  2:808814/150‑1 (MQ=255)

                            ccTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATa                                                                                                                         <  1:4546462/150‑1 (MQ=255)

                             cTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAg                                                                                                                        <  1:3945338/150‑1 (MQ=255)

                                gcgcATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCaa                                                                                                                     <  2:4942139/150‑1 (MQ=255)

                                 cgcATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAAc                                                                                                                    >  1:87812/1‑150 (MQ=255)

                                  gcATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACt                                                                                                                   >  1:1614504/1‑150 (MQ=255)

                                   cATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTc                                                                                                                  >  2:144434/1‑150 (MQ=255)

                                     tGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTccc                                                                                                                <  1:4801197/150‑1 (MQ=255)

                                      ggCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCg                                                                                                               >  1:3416607/1‑150 (MQ=255)

                                        cgcTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGcc                                                                                                             >  2:1781168/1‑150 (MQ=255)

                                              aCCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAgg                                                                                                       >  2:5199132/1‑150 (MQ=255)

                                               ccACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGc                                                                                                      >  2:3067547/1‑150 (MQ=255)

                                               ccACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGc                                                                                                      <  1:1239639/150‑1 (MQ=255)

                                                  cTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAATTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCggg                                                                                                   <  1:4969571/150‑1 (MQ=255)

                                                  cTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCggg                                                                                                   <  2:3301948/150‑1 (MQ=255)

                                                    tCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAg                                                                                                 <  1:1045518/150‑1 (MQ=255)

                                                    tCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAg                                                                                                 >  2:4686798/1‑150 (MQ=255)

                                                      aaCCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCCAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGtt                                                                                               <  2:646541/150‑1 (MQ=255)

                                                      aaCCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGtt                                                                                               <  2:2449404/150‑1 (MQ=255)

                                                       accaccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGttt                                                                                              <  1:100575/150‑1 (MQ=255)

                                                       accaccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGttt                                                                                              <  1:964232/150‑1 (MQ=255)

                                                         caccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGttttt                                                                                            <  1:1562226/150‑1 (MQ=255)

                                                         caccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGttttt                                                                                            <  1:1124218/150‑1 (MQ=255)

                                                         caccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGttttt                                                                                            <  1:3000772/150‑1 (MQ=255)

                                                         caccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTGAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGttttt                                                                                            >  2:5292075/1‑150 (MQ=255)

                                                           ccTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGACTTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTa                                                                                          >  1:1470150/1‑150 (MQ=255)

                                                             tGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACt                                                                                        <  2:392447/150‑1 (MQ=255)

                                                                 cAACCAGCTCAGGTAACGCCTGCCCCATAACTGTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTGCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATg                                                                                    <  2:5080489/150‑1 (MQ=255)

                                                                  aaCCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGc                                                                                   >  2:2098184/1‑150 (MQ=255)

                                                                        cTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGta                                                                             <  2:289265/150‑1 (MQ=255)

                                                                         tCAGGTAACGCCGGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTGAGtat                                                                            <  2:92120/150‑1 (MQ=255)

                                                                            ggTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATg                                                                         <  2:1864531/150‑1 (MQ=255)

                                                                             gTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATgg                                                                        >  1:1390228/1‑150 (MQ=255)

                                                                               aaCGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATgggg                                                                      >  1:3033925/1‑150 (MQ=255)

                                                                               aaCGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATgggg                                                                      >  1:845556/1‑150 (MQ=255)

                                                                                 cGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGAc                                                                    >  1:253752/1‑150 (MQ=255)

                                                                                 cGCCTGCCCCATAACTGTGTAAGTAAAATGAGTAGGGAATTCAAGCAGGTCGTTAAGTTTGGTTTTCTTGTCAGCTCCGGCGTATCGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGAc                                                                    <  1:5144408/150‑1 (MQ=255)

                                                                                          cATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTac                                                           <  1:1114093/150‑1 (MQ=255)

                                                                                           aTAACTCTGTAACTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTGAGTATATGGGGACGGAAATTaca                                                          <  2:772105/150‑1 (MQ=255)

                                                                                                  tGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCaa                                                   >  2:2803752/1‑150 (MQ=255)

                                                                                                     aagtaaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAgtg                                                <  1:4528338/150‑1 (MQ=255)

                                                                                                        taaAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGttt                                             <  1:3505565/150‑1 (MQ=255)

                                                                                                           aaGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCTGCGTAACGTAATTAAATAGCAACTCCCGCCAGTAGGCGTGAGTTTTTTACTGATGCTTATTATATGGGGACGGAAATTACACTTTCAAGTGTTTAAt                                          <  2:4489983/150‑1 (MQ=255)

                                                                                                           aaGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAAt                                          <  1:1729729/150‑1 (MQ=255)

                                                                                                             ggAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAAttt                                        <  1:3274329/150‑1 (MQ=255)

                                                                                                              gaggagGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAAtttt                                       <  2:1589196/150‑1 (MQ=255)

                                                                                                              gAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAAtttt                                       <  2:1614485/150‑1 (MQ=255)

                                                                                                               aGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAAttttt                                      >  2:1107390/1‑150 (MQ=255)

                                                                                                                 tAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaa                                    >  1:1202099/1‑150 (MQ=255)

                                                                                                                  aGGGAAATCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaac                                   <  2:4728672/150‑1 (MQ=255)

                                                                                                                     gAATTCAAGCCGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaaccaa                                <  1:4017405/150‑1 (MQ=255)

                                                                                                                     gAATTCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaaccaa                                <  1:746775/150‑1 (MQ=255)

                                                                                                                      aaTTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaaacaaa                               >  1:3129126/1‑150 (MQ=255)

                                                                                                                        ttCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaaccaaacc                             <  2:4448177/150‑1 (MQ=255)

                                                                                                                        ttCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaaccaaacc                             <  1:5199512/150‑1 (MQ=255)

                                                                                                                        ttCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaaccaaacc                             <  1:362975/150‑1 (MQ=255)

                                                                                                                        ttCAAGCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTaaccaaacc                             <  1:4295552/150‑1 (MQ=255)

                                                                                                                            aGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTg                         <  2:4086480/150‑1 (MQ=255)

                                                                                                                             gCAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGa                        <  1:4171909/150‑1 (MQ=255)

                                                                                                                              cAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGAt                       >  1:2820231/1‑150 (MQ=255)

                                                                                                                              cAGTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTATTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGAt                       <  2:538609/150‑1 (MQ=255)

                                                                                                                               aGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATg                      <  1:4092288/150‑1 (MQ=255)

                                                                                                                               aGTTCGTTAAGTTTGGATTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGCTGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATg                      <  2:1457394/150‑1 (MQ=255)

                                                                                                                                gTTCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATgg                     <  1:180111/150‑1 (MQ=255)

                                                                                                                                 ttCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGa                    <  2:3551326/150‑1 (MQ=255)

                                                                                                                                  tCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGaa                   <  2:3882370/150‑1 (MQ=255)

                                                                                                                                  tCGTTAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGaa                   <  1:225211/150‑1 (MQ=255)

                                                                                                                                    gTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACa                 >  2:2242038/1‑150 (MQ=255)

                                                                                                                                      tAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAtt               <  1:1512419/150‑1 (MQ=255)

                                                                                                                                      tAAGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAtt               >  1:3911499/1‑150 (MQ=255)

                                                                                                                                       aaGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAtta              <  2:2575283/150‑1 (MQ=255)

                                                                                                                                       aaGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAtta              >  2:437240/1‑150 (MQ=255)

                                                                                                                                        agattGGTTAACATGTCAGCACCGGCGTAACGAAATTAAATAGCAACTCCCGCCAGAAGGAGGGAGTTATTTACTGATGCTTAGTATATGGGGAAGGAAATTACACATACAAGAGTATAATTTTTAACCAATCCACAGATCGAACAataa             >  2:2976636/4‑150 (MQ=255)

                                                                                                                                        aGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaa             <  1:4717493/150‑1 (MQ=255)

                                                                                                                                        aGTTTGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaa             >  1:1348850/1‑150 (MQ=255)

                                                                                                                                         gTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaat            >  2:2077961/1‑150 (MQ=255)

                                                                                                                                          tttGGTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaatt           <  2:1713408/150‑1 (MQ=255)

                                                                                                                                             ggTTTTCATGACAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATATTTAACCAAACCAGTGATGGAACAataatttaa        >  1:2607845/1‑150 (MQ=255)

                                                                                                                                              gTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaatttaat       <  2:620628/150‑1 (MQ=255)

                                                                                                                                              gTTTTAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACAttaatttaat       <  2:4640210/150‑1 (MQ=255)

                                                                                                                                               ttttCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATg      <  2:433434/150‑1 (MQ=255)

                                                                                                                                               ttttAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATg      <  1:1969733/150‑1 (MQ=255)

                                                                                                                                               ttttAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATg      <  1:2044869/150‑1 (MQ=255)

                                                                                                                                               ttttAATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATg      >  1:2426994/1‑150 (MQ=255)

                                                                                                                                                   cATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAAt  >  2:4837641/1‑150 (MQ=255)

GCTTGGTTTTACCGTTGGGGTGTAGTCACCTGGCGCATGGCGCTGTACCACTTCAACCACCTGATCAACCAGCTCAGGTAACGCCTGCCCCATAACTTTGTAAGTAAAAGGAGTAGGGAATTCAAGCAGTTCGTTAAGTTTGGTTTTCATGTCAGCTCCGGCGTAACGTAATTAAATAGCAACTCCCGCCAGAAGGCGGGAGTTTTTTACTGATGCTTAGTATATGGGGACGGAAATTACACTTTCAAGTGTTTAATTTTTAACCAAACCAGTGATGGAACATTAATTTAATGTAAT  >  CP009273/657949‑658245

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