BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009273	2,727,902	(C)_8→9	coding (492/732 nt)	yfiH ←	UPF0124 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,727,894	1	.	C	97.3%	460.2 / 4.9	110	coding (500/732 nt)	yfiH	UPF0124 family protein
Reads supporting (aligned to +/- strand): ref base . (1/2); new base C (56/51); total (57/53)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

Only 100 of 149 total aligned reads displayed.

CTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑T‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCC  >  CP009273/2727794‑2727989

   ccGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑T‑ccc                                                                                              >  2:3041405/1‑101 (MQ=255)

   ccGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑T‑ccc                                                                                              >  1:1437906/1‑101 (MQ=255)

    cGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑T‑cccc                                                                                             >  2:2459715/1‑101 (MQ=255)

    cGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑T‑cccc                                                                                             >  1:1444900/1‑101 (MQ=255)

       gCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑T‑ccccccc                                                                                          <  1:761445/101‑1 (MQ=255)

        cAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑T‑cccccccc                                                                                         <  1:2784388/101‑1 (MQ=255)

         aaGCTGATAAATNTCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTNTAGCGTCTANNGCCATNAACGCCTCGCGAA‑CC‑T‑ccccccccc                                                                                        >  2:31802/1‑100 (MQ=255)

         aaGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑T‑ccccccccc                                                                                        >  1:1709673/1‑100 (MQ=255)

           gCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑T‑cccccccccac                                                                                      <  2:1937095/101‑4 (MQ=255)

           gCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑T‑cccccccccac                                                                                      >  1:3296132/1‑98 (MQ=255)

            cTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAGACGCCTCGCGAA‑CC‑T‑cccccccccact                                                                                     <  2:1624674/101‑5 (MQ=255)

            cTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑T‑cccccccccact                                                                                     >  2:3510779/1‑97 (MQ=255)

            cTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑T‑cccccccccact                                                                                     >  1:2726479/1‑97 (MQ=255)

               aTAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCg                                                                                   <  2:3616342/101‑1 (MQ=255)

                  aaTATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCCCGCGAA‑CC‑TCCCCCCCCCACTTCGAAc                                                                                <  1:1294193/101‑1 (MQ=255)

                   atatCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAAcg                                                                               <  1:579634/101‑1 (MQ=255)

                     atCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAAcgcg                                                                             <  1:2270808/101‑1 (MQ=255)

                        cGCCAGATCCTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGt                                                                          <  2:2827750/101‑1 (MQ=255)

                           cAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAACCGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGa                                                                       <  1:624808/101‑1 (MQ=255)

                           cAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGa                                                                       <  2:616133/101‑1 (MQ=255)

                           cAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGa                                                                       >  2:3406091/1‑101 (MQ=255)

                            aGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGAc                                                                      >  2:2901664/1‑101 (MQ=255)

                            aGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGAc                                                                      >  1:3801929/1‑101 (MQ=255)

                            aGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGAc                                                                      >  1:3800400/1‑101 (MQ=255)

                               tACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCaa                                                                   <  2:2545573/101‑1 (MQ=255)

                                aCTTTTCAACATACTTAAAAAAAGCAGCACTTGCTTTAGCATTTTCTGCCCTTAACGCCTCGGTAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTAGACCAAt                                                                  <  2:3673627/101‑1 (MQ=255)

                                aCTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAAt                                                                  <  1:800681/101‑1 (MQ=255)

                                 cTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAAtt                                                                 >  2:227579/1‑101 (MQ=255)

                                 cTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAAtt                                                                 <  1:2750366/101‑1 (MQ=255)

                                  ttATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAACCGCCTCGCGAACCC‑C‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTg                                                                <  1:2550468/101‑1 (MQ=255)

                                  ttATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTg                                                                <  1:2789502/101‑1 (MQ=255)

                                  ttATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTg                                                                <  2:1775878/101‑1 (MQ=255)

                                  ttATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACCCGCGTGGACCAATTg                                                                >  1:2183176/1‑101 (MQ=255)

                                         cATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGccc                                                         <  1:1361276/101‑1 (MQ=255)

                                          aTGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCt                                                        >  1:442199/1‑101 (MQ=255)

                                           tGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTa                                                       >  1:3294746/1‑101 (MQ=255)

                                            gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGAGCAATTGCCGCCCCaaa                                                      >  2:3302195/1‑101 (MQ=255)

                                            gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTaa                                                      >  1:1621231/1‑101 (MQ=255)

                                            gCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTaa                                                      >  2:2852731/1‑101 (MQ=255)

                                             cTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAc                                                     >  1:296934/1‑101 (MQ=255)

                                             cTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAc                                                     >  2:2791355/1‑101 (MQ=255)

                                             cTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAc                                                     <  2:861165/101‑1 (MQ=255)

                                             cTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACCCGCGTGGACCAATTGCCGGCCCTAAc                                                     >  1:1114917/1‑101 (MQ=255)

                                              tgaatgaaAGCTGCATTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAcc                                                    <  2:3719052/101‑1 (MQ=255)

                                              tgaatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAcc                                                    >  1:1875218/1‑101 (MQ=255)

                                              tgaatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAAcc                                                    <  2:611553/101‑1 (MQ=255)

                                              tgaatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACAGGCTTGGGCCAATTGCCGGCCCTAAcc                                                    >  2:2457457/1‑101 (MQ=255)

                                               gaatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCa                                                   >  1:1057119/1‑101 (MQ=255)

                                                aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAg                                                  <  1:3771499/101‑1 (MQ=255)

                                                aatgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGCGGACCAATTGACGGCCCCAACCAg                                                  >  1:2267070/1‑101 (MQ=255)

                                                  tgaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGc                                                <  2:1511484/101‑1 (MQ=255)

                                                   gaaAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCg                                               <  2:3200594/101‑1 (MQ=255)

                                                       gCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGaa                                           >  1:1013650/1‑101 (MQ=255)

                                                        cTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAC‑CCGC‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAat                                          <  2:3862380/101‑1 (MQ=255)

                                                        cTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAat                                          <  2:1872226/101‑1 (MQ=255)

                                                          gCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAatat                                        <  1:352440/101‑1 (MQ=255)

                                                           cACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATAtt                                       >  2:2200682/1‑101 (MQ=255)

                                                           cACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCCAACACGCGTGGTCCAATTGCCGGCCCTAACCAGGCGAAAATAtt                                       >  1:2053964/1‑101 (MQ=255)

                                                            aCTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATAttt                                      >  2:3042376/1‑101 (MQ=255)

                                                             cTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATAtttt                                     <  1:2982854/101‑1 (MQ=255)

                                                               tGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTcc                                   >  1:3098320/1‑101 (MQ=255)

                                                               tGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTcc                                   <  1:2358702/101‑1 (MQ=255)

                                                                 cTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCgg                                 >  2:2949624/1‑101 (MQ=255)

                                                                 cTTTAGCCTCTACGGCGATAACCCCCTCGCGAACCC‑C‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCgg                                 <  2:1495231/101‑1 (MQ=255)

                                                                     aGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTa                             <  1:3143255/101‑1 (MQ=255)

                                                                        gTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCa                          >  2:884553/1‑101 (MQ=255)

                                                                        gTCTACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCa                          <  1:1338796/101‑1 (MQ=255)

                                                                        gTCTACTGCCATAAACGCCTCGCGAA‑CC‑T‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAg                         >  1:3696114/1‑101 (MQ=255)

                                                                          cTACTGCCATAAACGCCTCGCGAC‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGc                        <  2:46349/101‑1 (MQ=255)

                                                                           tACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTTTCAGCa                       >  1:2350430/1‑101 (MQ=255)

                                                                           tACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCa                       >  2:1721896/1‑101 (MQ=255)

                                                                           tACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCa                       >  1:2841186/1‑101 (MQ=255)

                                                                           tACTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCa                       <  1:473868/101‑1 (MQ=255)

                                                                            aCTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCaa                      >  1:716244/1‑101 (MQ=255)

                                                                             cTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCaaa                     >  1:1384164/1‑101 (MQ=255)

                                                                             cTGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCaaa                     <  2:680639/101‑1 (MQ=255)

                                                                              tGCCATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCaaaa                    >  1:127231/1‑101 (MQ=255)

                                                                                ccATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACa                  <  2:1520570/101‑1 (MQ=255)

                                                                                ccATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACa                  >  1:2024246/1‑101 (MQ=255)

                                                                                 cATAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAg                 <  2:3858502/101‑1 (MQ=255)

                                                                                  aTAAACGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAgg                >  1:3569132/1‑101 (MQ=255)

                                                                                    aaaCGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaa              <  2:2380498/101‑1 (MQ=255)

                                                                                    aaaCGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaa              <  2:3550314/101‑1 (MQ=255)

                                                                                    aaaCGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaa              <  2:707273/101‑1 (MQ=255)

                                                                                    aaaCGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaa              >  1:2696444/1‑101 (MQ=255)

                                                                                     aaCGCCTCGCGCC‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaaa             <  2:460791/101‑1 (MQ=255)

                                                                                     aaCGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaaa             <  2:1604405/101‑1 (MQ=255)

                                                                                     aaCGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaaa             <  1:1046099/101‑1 (MQ=255)

                                                                                     aaCGCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGGGGCCCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGaaa             >  2:1244635/1‑101 (MQ=255)

                                                                                        gCCTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCg          >  2:266002/1‑101 (MQ=255)

                                                                                          cTCGCGAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtc        >  2:1849957/1‑101 (MQ=255)

                                                                                           ccgcgAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtct       <  1:1806008/100‑1 (MQ=255)

                                                                                            cgcgAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtctc      >  2:3445624/1‑101 (MQ=255)

                                                                                            cgcgAA‑CC‑TCCCCCCCCCACTTCCAACGCGCGTGGGCCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGtctc      >  1:132171/1‑101 (MQ=255)

                                                                                             gcgAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTctct     <  1:2901664/101‑1 (MQ=255)

                                                                                              cgAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCtt    >  1:1511065/1‑101 (MQ=255)

                                                                                              cgAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCtt    <  1:1157037/101‑1 (MQ=255)

                                                                                              ccgcccc‑T‑CCCCCCCCACTGCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCtt    <  2:131703/96‑1 (MQ=255)

                                                                                               gAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTc   <  2:2875638/101‑1 (MQ=255)

                                                                                               gAA‑CC‑TCCCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTc   <  2:1437906/101‑1 (MQ=255)

CTGCCGGGCAAGCTGATAAATATCCGCCAGATACTTATCACCATGCTGAATGAAAGCTGCACTTGCTTTAGCGTCTACTGCCATAAACGCCTCGCGAA‑CC‑T‑CCCCCCCCACTTCGAACGCGCGTGGACCAATTGCCGGCCCTAACCAGGCGAGAATATTTTCCGGATTATCAGCAAAACAGGAAACCGTCTCTTCC  >  CP009273/2727794‑2727989

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 17 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`