BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006881	3,725,237	A→G	intergenic (‑69/‑460)	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,237	0	A	G	99.5%	1289.0 / NA	366	intergenic (‑69/‑460)	U712_19205/U712_19210	Uncharacterized protein ywcI/Minor extracellular protease vpr
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (168/196); minor base C (0/2); total (168/198)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 366 total aligned reads displayed.

TTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑C‑TTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTTTGTTGAAAGCT  >  CP006881/3725097‑3725380

 tttCTGTCCGGCGCCCCCCAGTTCATCAAAAAGACCATCCCCACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTCCTCAGCCGCAGCCCGCAAGAAAAg                                                                                                                                        <  1:1604922/150‑1 (MQ=255)

  tgcTGTCGGGCGCCCCCCAGTTCATCAAACACACCATCCACACCCGTGAGGATCCAAGCCAACGTGTCATAGAGAGACAGCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGa                                                                                                                                       <  2:3847702/148‑1 (MQ=255)

     tGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGCTTTCTCAGCCGCAGGCCGCAAGAAAAGACgg                                                                                                                                    >  1:3058112/1‑150 (MQ=255)

      gTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGc                                                                                                                                   >  2:39088/1‑150 (MQ=255)

      gTCCGGCGTCACCCAGTCCACCAAAAAGACCATCCACCCCCGTATAGATACAAGCAAACGTTTCTTAGAAAACCATCCCTCCCCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGc                                                                                                                                   <  2:3588867/150‑1 (MQ=255)

       tCCGGCGTCACCCAGTTCAACAAAAAGACCACCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGcc                                                                                                                                  <  2:528826/150‑1 (MQ=255)

        ccGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCg                                                                                                                                 >  1:1832992/1‑150 (MQ=255)

         cGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGa                                                                                                                                <  1:749968/150‑1 (MQ=255)

                aCCCAGCTCATCAAAAAGCCCATCCACACCCGTAAAGATACAGGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttcttt                                                                                                                         <  2:1007401/150‑1 (MQ=255)

                   cAGTCCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGAtttcttttct                                                                                                                      <  2:2132833/150‑1 (MQ=255)

                    aGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTtctc                                                                                                                     <  1:3372843/150‑1 (MQ=255)

                    aGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTtctc                                                                                                                     >  2:2259883/1‑150 (MQ=255)

                     gTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCg                                                                                                                    <  1:2845032/150‑1 (MQ=255)

                      ttCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGc                                                                                                                   >  2:3022551/1‑150 (MQ=255)

                       tcatcaAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCa                                                                                                                  <  1:2449526/150‑1 (MQ=255)

                         atcaAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCCAGAAAAGACGGCCGATTTCTTTTCTCGCaaa                                                                                                                >  2:955736/1‑150 (MQ=255)

                         atcaAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCaaa                                                                                                                <  2:1370670/150‑1 (MQ=255)

                          tcaAAAAGACCATCCACCCCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAAc                                                                                                               <  2:2781186/150‑1 (MQ=255)

                          tcaAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAAc                                                                                                               <  2:1728775/150‑1 (MQ=255)

                           caAAAAGACCATCCACCCCGGAAAAGATACAAGCAAACTTTTCTTAGAAAAACATCCCTCCGCTACTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTCTACTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACt                                                                                                              <  2:1610445/150‑1 (MQ=255)

                             aaaaGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAg                                                                                                            <  2:3481720/150‑1 (MQ=255)

                              aaaGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGGTTTCCTTTCTCGCAAACTAGc                                                                                                           >  1:319230/1‑150 (MQ=255)

                              aaaGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGc                                                                                                           <  1:957849/150‑1 (MQ=255)

                                aGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTCCTCGCAAACTAGCgg                                                                                                         >  2:1997291/1‑150 (MQ=255)

                                 gACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCTGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCGGATTTCTTTTCTCGCAAACTAGCGGc                                                                                                        >  1:1163940/1‑150 (MQ=255)

                                     aTCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTa                                                                                                    >  2:231702/1‑150 (MQ=255)

                                        cacaCCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCa                                                                                                 <  1:1576132/150‑1 (MQ=255)

                                         acacCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaa                                                                                                >  1:1710856/1‑150 (MQ=255)

                                          cacCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaaa                                                                                               >  2:1453973/1‑150 (MQ=255)

                                          cacCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaaa                                                                                               >  2:2908318/1‑150 (MQ=255)

                                           acCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCaaaa                                                                                              >  2:434946/1‑150 (MQ=255)

                                               gTAAAGATCCAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCa                                                                                          <  1:93216/150‑1 (MQ=255)

                                               gTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCa                                                                                          >  2:2503480/1‑150 (MQ=255)

                                                 aaaGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATg                                                                                        <  1:2503480/150‑1 (MQ=255)

                                                 aaaGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATg                                                                                        <  1:3660884/150‑1 (MQ=255)

                                                 aaaGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTACTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATg                                                                                        >  2:3714639/1‑150 (MQ=255)

                                                   aGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGGCGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGAc                                                                                      >  1:100601/1‑150 (MQ=255)

                                                    gATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCTGGCAG‑C‑CTTTT‑TAGCAGCCCGGTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGGTCTCTTTTCTCCCACACTAGCGGCCCCGGCAAAACCCATGAct                                                                                     >  1:2618360/1‑150 (MQ=255)

                                                      tACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCTGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAGGAAAAAACGGCCGACTTCTTTTCTCCCAAACTAGCGGCCCTAGCAAAACTCATGActct                                                                                   >  1:2095385/1‑150 (MQ=255)

                                                      tACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCCCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTTCCAAAACTCATGActct                                                                                   >  1:3121885/1‑150 (MQ=255)

                                                       ccAAGCAAACGTTTCATAGAAAACCATCCCTCCCCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGGTTTCTCAGCCGCAGCCCGCGAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTa                                                                                  <  2:3391865/149‑1 (MQ=255)

                                                           gCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGGAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGGCTCTATGAt                                                                              >  1:502209/1‑150 (MQ=255)

                                                            cAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATa                                                                             >  2:640328/1‑150 (MQ=255)

                                                              aaCGTTTCATAGCAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGa                                                                           <  1:2249432/150‑1 (MQ=255)

                                                               aCGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGAc                                                                          >  2:2392310/1‑150 (MQ=255)

                                                                cGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGAcc                                                                         <  1:1379336/150‑1 (MQ=255)

                                                                cGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGAcc                                                                         <  2:48612/150‑1 (MQ=255)

                                                                 gTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCg                                                                        <  2:1120933/150‑1 (MQ=255)

                                                                    tCATAGAAATACAACCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TCGCAGCCCGTTTTCCCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTgg                                                                     <  2:192875/150‑1 (MQ=255)

                                                                      aTAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGaa                                                                   <  1:851768/150‑1 (MQ=255)

                                                                       tAGAGAAGCATCCCTCCGCTTCGTCTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGaaa                                                                  <  2:3314605/150‑1 (MQ=255)

                                                                        aGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCCAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAag                                                                 >  1:1211818/1‑150 (MQ=255)

                                                                           aaaaCATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAgaga                                                              >  1:1993212/1‑150 (MQ=255)

                                                                                aTCCCTCCGCTTCTTTTGGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCt                                                         <  2:1559155/150‑1 (MQ=255)

                                                                                  cccTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGGTAGACCGTGGAAAGAGAATTCttt                                                       >  2:2553356/1‑150 (MQ=255)

                                                                                    cTCCGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAACCTAGCGGCTCTAGCAAAACCCATGACCCTATGATAGACCGTGGAAAGAGAATTCTTTGt                                                     <  2:3644050/150‑1 (MQ=255)

                                                                                      ccGCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTca                                                   <  2:1956899/150‑1 (MQ=255)

                                                                                        gCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTcaca                                                 >  2:1657889/1‑150 (MQ=255)

                                                                                        gCTTCTTTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTcaca                                                 <  2:2435101/150‑1 (MQ=255)

                                                                                           tcttTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGGTAGACCGTGGAAAGAGAATTCTTTGTCACaaac                                              >  1:2476201/1‑148 (MQ=255)

                                                                                            cttTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGGAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAAGTCTTTGGCACaacaa                                             >  1:1669344/1‑147 (MQ=255)

                                                                                            cttTTTGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCCCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAatat                                             <  2:3454865/150‑1 (MQ=255)

                                                                                             tttttGGCAGGCAGCC‑C‑TTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATg                                            <  2:753948/150‑1 (MQ=255)

                                                                                                ttGGCAGGCAG‑CTT‑TTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGGAAGAGAATTCTTTGTCACAATAGGAgg                                         >  1:1306426/1‑149 (MQ=255)

                                                                                                 tGGCAGGCAG‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGAAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTg                                        <  1:626489/150‑1 (MQ=255)

                                                                                                  ggcaggcaG‑C‑C‑TTTCTAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGTTTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACCCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGa                                       >  1:2164333/1‑150 (MQ=255)

                                                                                                   gcaggcag‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGaa                                      >  2:14221/1‑150 (MQ=255)

                                                                                                     aggcag‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGaaaa                                    >  2:2086950/1‑150 (MQ=255)

                                                                                                     aggcag‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGaaaa                                    <  2:247719/150‑1 (MQ=255)

                                                                                                        cag‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAGAAc                                 >  1:2215232/1‑150 (MQ=255)

                                                                                                        cag‑C‑CTTTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAAc                                 <  1:788009/150‑1 (MQ=255)

                                                                                                           agtctttt‑tAGCAGCCCGTTCTCCCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAAcc                                <  1:3152027/145‑1 (MQ=255)

                                                                                                            gcctttt‑tAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGGGAAAAAACCa                               >  1:1518366/5‑150 (MQ=255)

                                                                                                             cctttt‑tAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCaa                              <  1:1390383/147‑1 (MQ=255)

                                                                                                             cctttt‑tAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCCATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCaa                              >  1:3587574/4‑150 (MQ=255)

                                                                                                             cctttt‑tAGCAGCCCGTTTCCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCACTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCaa                              <  1:2315846/147‑1 (MQ=255)

                                                                                                              ctttt‑tAGCAGCCCGTTTTCGCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCCAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAAc                             <  1:1083987/148‑1 (MQ=255)

                                                                                                                ttt‑tAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTa                           >  2:1078407/1‑150 (MQ=255)

                                                                                                                ttt‑tAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTa                           >  1:3715509/1‑150 (MQ=255)

                                                                                                                 tt‑tAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAg                          <  2:1753882/150‑1 (MQ=255)

                                                                                                                  t‑tAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGt                         <  2:3120555/150‑1 (MQ=255)

                                                                                                                  t‑tAGCAGCCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCCCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGt                         <  1:2937926/150‑1 (MQ=255)

                                                                                                                  t‑tAGCAGCCCGTTTTCCCAGCCGCCGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGt                         <  1:3883100/150‑1 (MQ=255)

                                                                                                                  t‑tAGCAGCCCGTTTCCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGt                         <  2:1860909/150‑1 (MQ=255)

                                                                                                                    tAGCAGCCCGTTTTCCCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGGAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGtt                        <  2:2614587/150‑1 (MQ=255)

                                                                                                                    tAGCAGCCCGTTTTCCCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGtt                        <  2:884345/150‑1 (MQ=255)

                                                                                                                      gcagcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGGAAGGGGATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGtttt                      >  1:147867/1‑150 (MQ=255)

                                                                                                                       cagcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAAATCTTTGTCACAATATGAGTGAAAAAACCAACTAGttttt                     >  1:409037/1‑150 (MQ=255)

                                                                                                                       caccccGTTTTCCCAGCCGCAGCCCCCAAGAAAAGACGGCCGATTTCTTTTCACGCAACCTAGCGGCTCTAGCAAAACCCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGttttt                     <  2:2182210/147‑1 (MQ=255)

                                                                                                                        agcCCGTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAACGGAATTCTTTGTCACAATATGAGTGGAAAAACCAACTAGTTTTTa                    >  1:54018/1‑150 (MQ=255)

                                                                                                                            cGTTTTCTCAGCCGCAGCCCGCAAGAAAAGCCGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAg                <  1:2051631/150‑1 (MQ=255)

                                                                                                                             gTTTTCTCAGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGt               <  1:2738128/150‑1 (MQ=255)

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

                                                                                                                                    cagccccagcccGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTttgttg        <  1:2553356/150‑1 (MQ=255)

                                                                                                                                     aGCCGCAGCCCGCAAGAAAAGACGGCGGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTTTGTTGa       <  1:3816491/150‑1 (MQ=255)

                                                                                                                                     aGCCGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTTTGTTGa       <  1:3022551/150‑1 (MQ=255)

                                                                                                                                       ccGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAGACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTCGCCACAATATGAATGAGAAAACCAACTAGTTTTTAGAAGGTTTGTTGaaa     >  1:3314575/1‑150 (MQ=255)

                                                                                                                                       ccGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGGATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTTTGTTGaaa     >  1:1439321/1‑150 (MQ=255)

                                                                                                                                       ccGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTGTTTAGAAGGTTTGTTGaaa     >  1:1926471/1‑150 (MQ=255)

                                                                                                                                        cGCAGCCCGCAAGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGGTTTGTTGAAAg    >  2:3146374/1‑150 (MQ=255)

TTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAACATCCCTCCGCTTCTTTTTGGCAGGCAG‑C‑C‑TTT‑TAGCAGCCCGTTTTCTCAGCCGCAGCCCACGGGAAAAGACGGCCGATTTCTTTTCTCGCAAACTAGCGGCTCTAGCAAAACTCATGACTCTATGATAGACCGTGGAAAGAGAATTCTTTGTCACAATATGAGTGAAAAAACCAACTAGTTTTTAGAAGTTTTGTTGAAAGCT  >  CP006881/3725097‑3725380

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