BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	647,422	G→T	100%	intergenic (‑110/‑269)	citC ← / → citA	[citrate [pro‑3S]‑lyase] ligase/sensory histidine kinase in two‑component regulatory system with CitB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	647,422	0	G	T	98.0%	822.4 / 6.5	253	intergenic (‑110/‑269)	citC/citA	[citrate [pro‑3S]‑lyase] ligase/sensory histidine kinase in two‑component regulatory system with CitB
Reads supporting (aligned to +/- strand): ref base G (1/4); new base T (129/118); total (131/122)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 253 total aligned reads displayed.

AACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGT  >  CP009273/647276‑647567

 aCGTTTTACGCGGGTGAAAATATCATTGCCGAACCTAATAAATAGGATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGta                                                                                                                                                <  2:568838/149‑1 (MQ=255)

 aCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGta                                                                                                                                                <  1:1386054/149‑1 (MQ=255)

 aCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGta                                                                                                                                                <  1:147195/149‑1 (MQ=255)

   gTTGTACGCGGGTGAAAATATCAGTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGtata                                                                                                                                              <  2:4436685/149‑1 (MQ=255)

    ttttACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTatat                                                                                                                                             <  1:783148/149‑1 (MQ=255)

     tttACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATc                                                                                                                                            <  1:3947787/149‑1 (MQ=255)

      ttACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCt                                                                                                                                           <  1:3202074/149‑1 (MQ=255)

        aCGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCttt                                                                                                                                         <  1:3012228/149‑1 (MQ=255)

          gcgGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCttttt                                                                                                                                       <  1:2219193/149‑1 (MQ=255)

              gTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTcc                                                                                                                                   <  2:3270201/149‑1 (MQ=255)

               tGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCt                                                                                                                                  >  2:670247/1‑149 (MQ=255)

               tGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTACTCCt                                                                                                                                  <  2:2317167/149‑1 (MQ=255)

               tGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTACTCCt                                                                                                                                  <  2:2643574/149‑1 (MQ=255)

                gAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCtt                                                                                                                                 <  1:85656/149‑1 (MQ=255)

                    atatCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGc                                                                                                                             <  1:3816146/149‑1 (MQ=255)

                    atatCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGc                                                                                                                             <  1:2631686/149‑1 (MQ=255)

                    atatCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGc                                                                                                                             <  1:2710841/149‑1 (MQ=255)

                      atCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGTGGGTGATCGATGTATATCTTTTTCTCCTTTCGCta                                                                                                                           <  1:937876/149‑1 (MQ=255)

                      atCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCta                                                                                                                           <  2:3688911/149‑1 (MQ=255)

                        cATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCtata                                                                                                                         >  1:3945793/1‑149 (MQ=255)

                        cATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGGTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCtata                                                                                                                         >  1:128955/1‑149 (MQ=255)

                          ttGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACa                                                                                                                       >  2:237876/1‑149 (MQ=255)

                             ccGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATcc                                                                                                                    >  1:1751035/1‑149 (MQ=255)

                                aaCATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTaa                                                                                                                 >  1:3336106/1‑149 (MQ=255)

                                 aCATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAg                                                                                                                <  2:4608217/149‑1 (MQ=255)

                                  cATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAgg                                                                                                               <  2:4249585/149‑1 (MQ=255)

                                   ataataAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGa                                                                                                              >  1:3090943/1‑149 (MQ=255)

                                    taataaATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGAAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAg                                                                                                             >  1:2428375/1‑149 (MQ=255)

                                     aataaataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGt                                                                                                            >  2:1284411/1‑149 (MQ=255)

                                         aataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTAttt                                                                                                        <  2:2468899/149‑1 (MQ=255)

                                         aataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTAttt                                                                                                        >  2:4109158/1‑149 (MQ=255)

                                           taGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCg                                                                                                      >  1:4184963/1‑149 (MQ=255)

                                            aGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCgg                                                                                                     <  1:2343913/149‑1 (MQ=255)

                                             gTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTTGCTATACATCCTAAGGAGTATTTCGGc                                                                                                    >  2:3101411/1‑149 (MQ=255)

                                             gTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGc                                                                                                    >  1:192292/1‑149 (MQ=255)

                                                    gAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAAt                                                                                             >  2:3057460/1‑149 (MQ=255)

                                                     aaGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAAtt                                                                                            <  2:661124/149‑1 (MQ=255)

                                                       ggTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAAtttt                                                                                          >  2:1038273/1‑149 (MQ=255)

                                                       ggTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAAtttt                                                                                          >  2:1172207/1‑149 (MQ=255)

                                                       ggTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAAtttt                                                                                          >  1:3034108/1‑149 (MQ=255)

                                                          gCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGat                                                                                       <  2:1214683/149‑1 (MQ=255)

                                                           cATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGatt                                                                                      <  1:1924600/149‑1 (MQ=255)

                                                            aTGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGattt                                                                                     >  1:3957999/1‑149 (MQ=255)

                                                                 gttATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCCAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTc                                                                                >  1:4099852/1‑149 (MQ=255)

                                                                  ttATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTca                                                                               >  1:4710879/1‑149 (MQ=255)

                                                                   tATCGATTTGCAACGAATGTTGTTCAATGGTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTcac                                                                              <  1:4028906/149‑1 (MQ=255)

                                                                        aTTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATaga                                                                         >  1:2944498/1‑149 (MQ=255)

                                                                        aTTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATaga                                                                         <  2:3578457/149‑1 (MQ=255)

                                                                        aTTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATaga                                                                         <  2:4380658/149‑1 (MQ=255)

                                                                        aTTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATaga                                                                         <  2:614999/149‑1 (MQ=255)

                                                                        aTTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATaga                                                                         <  2:4380927/149‑1 (MQ=255)

                                                                         tttGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATagag                                                                        >  2:441343/1‑149 (MQ=255)

                                                                            gCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTa                                                                     >  2:2606175/1‑149 (MQ=255)

                                                                            gCAACGAATGTTGTTCAATGGTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTa                                                                     <  2:3742828/149‑1 (MQ=255)

                                                                             cAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTGATGTTCACTTGGGAGAAAGGGGGTGATCGATGTCTATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATGTCACATAGAGTTAg                                                                    <  2:727161/149‑1 (MQ=255)

                                                                              aaCGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGt                                                                   <  2:3925142/149‑1 (MQ=255)

                                                                              aaCGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTAACATAGAGTTAGt                                                                   >  2:4176290/1‑149 (MQ=255)

                                                                               aCGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTg                                                                  >  1:2507566/1‑149 (MQ=255)

                                                                                 gAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGt                                                                <  1:4186354/149‑1 (MQ=255)

                                                                                 gAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGt                                                                <  1:3496838/149‑1 (MQ=255)

                                                                                    tgttgtTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGtttt                                                             >  1:241714/1‑149 (MQ=255)

                                                                                    tgttgtTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGtttt                                                             <  2:1921488/149‑1 (MQ=255)

                                                                                     gttgttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGttttc                                                            <  1:4777793/149‑2 (MQ=255)

                                                                                       tgttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtcta                                                          <  1:1138087/149‑1 (MQ=255)

                                                                                        gttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATATAGTTAGTGGTTTtctat                                                         >  1:4365441/1‑149 (MQ=255)

                                                                                        gttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctat                                                         >  1:3789288/1‑149 (MQ=255)

                                                                                         ttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatt                                                        <  2:2075086/149‑1 (MQ=255)

                                                                                          tCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctattt                                                       <  1:1725265/149‑1 (MQ=255)

                                                                                             aTGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatttatt                                                    <  2:1621767/149‑1 (MQ=255)

                                                                                             aTGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatttatt                                                    <  1:4612862/149‑1 (MQ=255)

                                                                                                ttGCAAACTGATAACCTTTTATTTTCACTTAGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTaa                                                 >  2:1773002/1‑149 (MQ=255)

                                                                                                      aCTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGAttt                                           >  2:908458/1‑149 (MQ=255)

                                                                                                      aCTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGAttt                                           >  2:3361002/1‑149 (MQ=255)

                                                                                                       cTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGAtttt                                          >  2:3517221/1‑149 (MQ=255)

                                                                                                        tGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTa                                         >  1:254300/1‑149 (MQ=255)

                                                                                                        tGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTa                                         >  1:411857/1‑149 (MQ=255)

                                                                                                          aTAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAg                                       <  2:4747191/149‑1 (MQ=255)

                                                                                                              ccTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGtttt                                   <  2:324696/149‑1 (MQ=255)

                                                                                                              ccTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGtttt                                   <  1:1728293/149‑1 (MQ=255)

                                                                                                              ccTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCCTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGtttt                                   >  2:83084/1‑149 (MQ=255)

                                                                                                               cTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGGATATCTTTTTCTCCTTTAGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTACACATAGAGATAGTGGTTTTCTATATATTTAATGATTTTAAGtttat                                  >  1:3455421/1‑147 (MQ=255)

                                                                                                                ttttATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGtttttt                                 >  1:2799247/1‑149 (MQ=255)

                                                                                                                   tATTTTCACTTAGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGTTTTttaat                              <  1:2088098/149‑1 (MQ=255)

                                                                                                                   tATTTTCACTTAGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGTTTTttaat                              <  1:4421653/149‑1 (MQ=255)

                                                                                                                        tCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATg                         <  1:2100546/149‑1 (MQ=255)

                                                                                                                          aCTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTa                       >  2:1606317/1‑149 (MQ=255)

                                                                                                                            ttGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAAt                     >  1:1463054/1‑149 (MQ=255)

                                                                                                                            ttGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAAt                     >  1:1463266/1‑149 (MQ=255)

                                                                                                                             tGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAAtt                    >  2:2246076/1‑149 (MQ=255)

                                                                                                                              gggAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTa                   <  1:479310/149‑1 (MQ=255)

                                                                                                                              gggAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTa                   >  2:1448678/1‑149 (MQ=255)

                                                                                                                              gggAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTa                   <  1:877079/149‑1 (MQ=255)

                                                                                                                               ggAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTAc                  >  1:1352272/1‑149 (MQ=255)

                                                                                                                                  gaAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTAAGaa               >  1:3404669/1‑149 (MQ=255)

                                                                                                                                    aaGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAAt             <  1:1081798/149‑1 (MQ=255)

                                                                                                                                     aGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATg            >  2:3373295/1‑149 (MQ=255)

                                                                                                                                      gggggTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGa           >  2:2855565/1‑149 (MQ=255)

                                                                                                                                      gggggTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGa           >  1:1603041/1‑149 (MQ=255)

                                                                                                                                           tGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAACTTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGc      >  2:387535/1‑149 (MQ=255)

                                                                                                                                           tGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAGATGACTCGc      >  2:2779647/1‑149 (MQ=255)

AACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGT  >  CP009273/647276‑647567

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