BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	3,725,239	2 bp→AA	100%	intergenic (‑71/‑457)	U712_19205 ← / → U712_19210	Uncharacterized protein ywcI/Minor extracellular protease vpr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	3,725,239	0	G	A	99.8%	1662.5 / ‑5.6	478	intergenic (‑71/‑458)	U712_19205/U712_19210	Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand): ref base G (0/0); major base A (239/238); minor base T (1/0); total (240/238)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*	CP006881	3,725,240	0	G	A	100.0%	1667.4 / ‑6.6	478	intergenic (‑72/‑457)	U712_19205/U712_19210	Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (239/237); total (240/238)
Rejected as polymorphism: E-value score below prediction cutoff.
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 478 total aligned reads displayed.

TTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑C‑TT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTTTGTTGAAAGCT  >  CP006881/3725097‑3725380

||

ttttCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGaaaa                                                                                                                                         >  2:5649350/1‑150 (MQ=255)

    cTGTCCGGCGTCACCCAGTCCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCACACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTCCCCAGCCGCAGCCCGCAAGAAAAGCCg                                                                                                                                     <  2:5718770/150‑1 (MQ=255)

     tgtgCGGCGTCCCCCAGTTCATCAAAAAGACCATCCACAGCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCAGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACgg                                                                                                                                    <  1:1454664/150‑1 (MQ=255)

        ccGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCg                                                                                                                                 >  2:324266/1‑150 (MQ=255)

        ccGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCg                                                                                                                                 <  2:1063915/150‑1 (MQ=255)

             gTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttc                                                                                                                            >  2:1463851/1‑150 (MQ=255)

              tCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttct                                                                                                                           >  2:1877382/1‑150 (MQ=255)

               cACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttctt                                                                                                                          >  2:428950/1‑150 (MQ=255)

                aCCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttcttt                                                                                                                         >  2:3553654/1‑150 (MQ=255)

                   cAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttcttttct                                                                                                                      >  1:2370497/1‑150 (MQ=255)

                   cAGTTCACCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAGCGTTTCAGAGAAAATCATCCCTCCGCTCCTTTTTGGCAGGCT‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttcttttct                                                                                                                      <  2:3407162/150‑1 (MQ=255)

                        caccaAGAAGACCATCCACACCCGTAAAGATCCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCCTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGTCCGATTTCTTTTCTCGCaa                                                                                                                 <  2:1292767/147‑1 (MQ=255)

                         atcaAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCaaa                                                                                                                >  1:4623750/1‑150 (MQ=255)

                             aaaaGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAg                                                                                                            >  2:3623886/1‑150 (MQ=255)

                             aaaaGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAg                                                                                                            >  1:4708094/1‑150 (MQ=255)

                              aaaGACCATACACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGc                                                                                                           <  2:5814162/150‑1 (MQ=255)

                               aaGCCCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTCCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTCCTCGCAAACTAGCg                                                                                                          <  2:425014/150‑1 (MQ=255)

                               aaGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGGCGGCCGATTTCTTTTCTCGGAAACTAGCg                                                                                                          >  1:2975605/1‑150 (MQ=255)

                                  aCCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGct                                                                                                       >  1:672231/1‑150 (MQ=255)

                                    cATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCTGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCTGAAACTAGCGGttct                                                                                                     >  1:1732236/1‑150 (MQ=255)

                                    cATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGctct                                                                                                     >  1:5721835/1‑150 (MQ=255)

                                     aTCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCTGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTa                                                                                                    >  2:3996177/1‑150 (MQ=255)

                                     aTCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTa                                                                                                    >  2:2237906/1‑150 (MQ=255)

                                              cGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGGAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTc                                                                                           >  1:1058647/1‑150 (MQ=255)

                                               gTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCa                                                                                          <  1:1373951/150‑1 (MQ=255)

                                                  aaGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGa                                                                                       >  1:3277331/1‑150 (MQ=255)

                                                   aGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAc                                                                                      >  1:29394/1‑150 (MQ=255)

                                                   aGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAc                                                                                      >  2:5479038/1‑150 (MQ=255)

                                                      tACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGActct                                                                                   <  1:1918439/150‑1 (MQ=255)

                                                           gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCGGGCA‑GC‑CTTT‑TTAGCAGCCCGGTTTCTCAGCCGCCGCCCCCAGGAAAAGACGGCCGAATTCTTTTCTCGCAAACTAGCGGGTCTAGCGAAACTCACGACCCTATGAt                                                                              >  2:2863312/1‑150 (MQ=255)

                                                           gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGAt                                                                              <  1:3417830/150‑1 (MQ=255)

                                                           gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGAt                                                                              >  1:204832/1‑150 (MQ=255)

                                                           gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGAt                                                                              >  2:5298390/1‑150 (MQ=255)

                                                              aaCGTTTCCTAGAAAAACATTCCTCCGCTTCTTTTTGGCAGGCA‑GC‑C‑TTCTTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGa                                                                           <  1:4224259/150‑1 (MQ=255)

                                                               acttttCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAgcc                                                                          <  1:3240416/147‑3 (MQ=255)

                                                                cGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGAcc                                                                         <  2:2565998/150‑1 (MQ=255)

                                                                cGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGAcc                                                                         <  2:4708094/150‑1 (MQ=255)

                                                                 gTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCg                                                                        >  1:5601241/1‑150 (MQ=255)

                                                                 gTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGCTTTCTCAGCCGCAACCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCGCTAGCAAAACTCATGACTCTATGATAggacg                                                                        >  1:4122295/1‑146 (MQ=255)

                                                                      aTAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGaa                                                                   <  2:5003953/150‑1 (MQ=255)

                                                                        aGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAag                                                                 <  2:5230271/150‑1 (MQ=255)

                                                                        aGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGAGGAAag                                                                 <  1:5655658/150‑1 (MQ=255)

                                                                         gAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAaga                                                                >  2:1704275/1‑150 (MQ=255)

                                                                           aaaaCATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAgaga                                                              >  2:892486/1‑150 (MQ=255)

                                                                             aaCATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAAt                                                            <  2:161857/150‑1 (MQ=255)

                                                                             aaCATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAAt                                                            >  2:2271256/1‑150 (MQ=255)

                                                                                atgccTCCGCTTCTTTTTGGCAGGCAGGC‑T‑TT‑TTAGGAGCCCGTTTTCTCAGCCGCACCCCGCAAGAAAAGACGACCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAAATCt                                                         >  1:4659067/4‑150 (MQ=255)

                                                                                aTCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCt                                                         <  2:4810721/150‑1 (MQ=255)

                                                                                aTCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCt                                                         <  2:3297117/150‑1 (MQ=255)

                                                                                aTCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCt                                                         >  1:2595568/1‑150 (MQ=255)

                                                                                aTCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAAACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGGCTCTATGATAGACCGTGGAAAGAGAATTCt                                                         >  1:1213909/1‑150 (MQ=255)

                                                                                 tCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCGGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTACCCCCACTAGCACAACTCATGACTCTATGATAGACCGTGCAAAGACAATTCtt                                                        >  1:3786177/1‑150 (MQ=255)

                                                                                  cccTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCttt                                                       <  2:4398608/150‑1 (MQ=255)

                                                                                   ccTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCt                                                         <  2:3329451/147‑1 (MQ=255)

                                                                                   ccTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTg                                                      <  1:5097437/150‑1 (MQ=255)

                                                                                    cTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGt                                                     >  2:3947501/1‑150 (MQ=255)

                                                                                    cTCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGt                                                     >  1:5524457/1‑150 (MQ=255)

                                                                                     tCCGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTc                                                    >  1:182465/1‑150 (MQ=255)

                                                                                      ccGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTca                                                   >  2:3842741/1‑150 (MQ=255)

                                                                                       cGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTcac                                                  <  1:3364271/150‑1 (MQ=255)

                                                                                       cGCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTcac                                                  <  2:4092719/150‑1 (MQ=255)

                                                                                        gCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTcaca                                                 >  2:2520393/1‑150 (MQ=255)

                                                                                        gCTTCTTTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTcaca                                                 <  1:5803114/150‑1 (MQ=255)

                                                                                          ttcttTTTGGCGGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAAATCATGACTCTATGATAGACCGTGGAAAGAGAATTGTTTGTCACAat                                               >  1:1527033/1‑150 (MQ=255)

                                                                                          ttcttTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAat                                               <  1:4946547/150‑1 (MQ=255)

                                                                                           tcttTTTGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAata                                              >  1:3629410/1‑150 (MQ=255)

                                                                                            cttTTTGGCAGGCA‑GCTT‑TT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCCTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAatat                                             <  1:517444/150‑1 (MQ=255)

                                                                                              ttttGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTGCTCAGCCGGAGCCCCCAAGAAAAGACGGCCGATTCCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGa                                           <  1:2863312/150‑1 (MQ=255)

                                                                                               tttGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAg                                          >  1:2494493/1‑150 (MQ=255)

                                                                                               tttGGCAGGCA‑GC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAg                                          <  2:50662/150‑1 (MQ=255)

                                                                                                   gcaggca‑gC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGGCGATTTCTTTTTTCGCAAACTAGCGGCTCTAGAAAAACTCATAACTCTATGATAGACCGTGGAAAAAGTATTTTTTGTCACAATATGAGTGaa                                      >  1:1665603/1‑150 (MQ=255)

                                                                                                   gcaggca‑gC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGaa                                      <  2:4169190/150‑1 (MQ=255)

                                                                                                    caggca‑gC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGaaa                                     <  1:4321015/150‑1 (MQ=255)

                                                                                                    caggca‑gC‑CTTT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGaaa                                     >  2:291841/1‑150 (MQ=255)

                                                                                                           agccttt‑ttAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAAcc                                <  1:4493663/145‑1 (MQ=255)

                                                                                                           agccttt‑ttAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAAcc                                >  2:2095338/6‑150 (MQ=255)

                                                                                                               ttt‑ttAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAaactaact                            >  1:2178137/2‑150 (MQ=255)

                                                                                                               ttt‑ttAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACt                            <  1:3615920/149‑1 (MQ=255)

                                                                                                               ttt‑ttAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACt                            <  2:4623750/149‑1 (MQ=255)

                                                                                                               ttt‑ttAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACt                            <  2:1827601/149‑1 (MQ=255)

                                                                                                                tt‑ttAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTa                           <  1:3069510/150‑1 (MQ=255)

                                                                                                                 t‑ttAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAg                          >  2:2419987/1‑150 (MQ=255)

                                                                                                                   ttAGCAGCCCGTTTCCTCACCCGCAGCCCGCAAGAAAATACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGt                         <  1:2649769/150‑1 (MQ=255)

                                                                                                                     agcagcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGGCCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGttt                       >  2:711200/1‑150 (MQ=255)

                                                                                                                     agcagcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGttt                       >  1:53562/1‑150 (MQ=255)

                                                                                                                          cccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGa                  <  2:2964169/150‑1 (MQ=255)

                                                                                                                           ccGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGaa                 <  1:3115422/150‑1 (MQ=255)

                                                                                                                               tttCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGttt             <  2:624019/150‑1 (MQ=255)

                                                                                                                                  ctcAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTttgt          >  2:3753902/1‑150 (MQ=255)

                                                                                                                                   tcAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTttgtt         <  1:789183/150‑1 (MQ=255)

                                                                                                                                   tcAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACGAGTTTTTAGAAGGTttgtt         >  2:2006156/1‑150 (MQ=255)

                                                                                                                                    cAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTttgttg        <  2:686452/150‑1 (MQ=255)

                                                                                                                                    cAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTttgttg        >  2:1346676/1‑150 (MQ=255)

                                                                                                                                     aGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTTTGTTGa       >  1:2270353/1‑150 (MQ=255)

                                                                                                                                      gCCGCAGCTCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTTTGTTGaa      <  2:2783807/150‑1 (MQ=255)

                                                                                                                                       ccGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTTTGTTGaaa     <  2:5386114/150‑1 (MQ=255)

                                                                                                                                        cGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTTTGTTGAAAg    <  2:5326112/150‑1 (MQ=255)

                                                                                                                                        cGCAGCCCGCAAGAAAACCCGGCCGACTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCTCAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTTTGTTGAAAg    <  2:1413872/150‑1 (MQ=255)

                                                                                                                                          cAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGTCTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTTTGTTGAAAGCt  <  1:1292032/150‑1 (MQ=255)

                                                                                                                                          cAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTTTGTTGAAAGCt  <  2:3034338/150‑1 (MQ=255)

||

TTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCA‑GC‑C‑TT‑TTAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTTTGTTGAAAGCT  >  CP006881/3725097‑3725380

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