BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,917	C→T	58.8%	intergenic (‑1000/‑168)	USA300HOU_1794 ← / → USA300HOU_1796	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,929,917	0	C	T	58.8%	82.0 / 169.3	153	intergenic (‑1000/‑168)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (32/31); new base T (51/38); total (83/70)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 154 total aligned reads displayed.

GTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAA  >  CP000730/1929777‑1930054

gTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑T‑AA‑G‑G‑TGAGGCGGGTTTTCGTTCACATCATGGCACCAGTCATCAGTGGCTGTGcc                                                                                                                                                  >  1:524609/1‑141 (MQ=255)

   aCTGTGCTTGTATATTTGCCTATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGGAGAG‑AGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAtt                                                                                                                                               <  2:315977/141‑1 (MQ=255)

     tgtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGc                                                                                                                                             <  1:60674/141‑1 (MQ=255)

     tgtgCTTGTATATTCGCCGATGATATAATAGGAGAGTCGCCTATCTCTCCGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGc                                                                                                                                             <  2:28513/141‑1 (MQ=255)

     tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGAGCTGTTGCg                                                                                                                                            <  1:567487/139‑1 (MQ=255)

     tatgCTTGTATATGCACAGATGGTAAAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                                                                                                                                            <  1:698596/139‑1 (MQ=255)

      gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGGGGCG‑GA‑C‑AA‑G‑C‑AGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATGGCg                                                                                                                                            >  1:177912/1‑141 (MQ=255)

      gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCT‑AAAGAAGCAAAAGTGAGCCG‑AGTGAACGTCGAG‑TA‑T‑AA‑G‑G‑TGATGCTTGTTTTCGTTCACATCAGAGCACCAGTCATCAGTGGCTGTACCATTGCg                                                                                                                                            >  1:370522/1‑141 (MQ=255)

      atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                                                                           <  1:105831/140‑1 (MQ=255)

        gCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtt                                                                                                                                          >  1:691625/1‑141 (MQ=255)

         cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                                                                                        >  2:642691/1‑141 (MQ=255)

           gttattTTGGCCGATGATATAATAGTAGAGTCGCTTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑T‑AA‑G‑T‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttt                                                                                                                                       <  2:395127/139‑1 (MQ=255)

            gTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattata                                                                                                                                     <  1:538346/141‑7 (MQ=255)

             tataTTCACAGATGGTATAATAATAGATTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACCAAGAG‑AGGAGGTG‑CA‑T‑TA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                                                                                                    <  1:365937/141‑8 (MQ=255)

               tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑T‑TA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACtt                                                                                                                                  >  1:1491/9‑141 (MQ=255)

               tatTTACAGATGGTGTAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatactt                                                                                                                                  <  1:198709/141‑7 (MQ=255)

                attcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCCGTCATCAGTGGCTGTGCTGTTGCGTATTATACttt                                                                                                                                 >  2:170795/8‑140 (MQ=255)

                aattaccaatgatATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCtt                                                                                                                                  <  1:333251/133‑1 (MQ=255)

                 tttGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCA‑TCAATTTAGACGCAGAG‑ACACGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCACCAGTGACTGTTCCCTTGCGtttttttactt                                                                                                                                 >  2:174548/1‑137 (MQ=255)

                 ttcacagatgatATAAAAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                                                                                <  2:23651/135‑1 (MQ=255)

                  tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                                                                               >  1:736519/6‑141 (MQ=255)

                    acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                                                                             >  1:687411/4‑141 (MQ=255)

                    acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGTATTATACTTATTgg                                                                                                                             <  2:585235/138‑1 (MQ=255)

                    acagatgggataataataGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑T‑TA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                                                                             <  1:71719/132‑1 (MQ=255)

                     cagatggtATAATAATAGAGAAGCATAACTCTCTTACG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                            >  1:767767/3‑141 (MQ=255)

                      agatggtATAATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                                                                           >  2:211307/2‑141 (MQ=255)

                         tggtataATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CC‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                        >  1:297551/4‑141 (MQ=255)

                         tggtataATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                        <  2:425958/138‑4 (MQ=255)

                         tggtataATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                        >  1:430400/4‑141 (MQ=255)

                         tggtataATAATAGAGGCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGCAGAG‑AGGAGGAG‑CATT‑AC‑T‑T‑TGCTAATCATTTTAGTTCACATAATAGCACGAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                        >  1:646422/4‑141 (MQ=255)

                          ggtataATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                       >  1:409270/3‑141 (MQ=255)

                          ggtataATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                       >  1:582952/3‑141 (MQ=255)

                          ggtataAAAATAGAGTCGCCAATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCGCATCATAGCACCAGTCATCAGTGGCTGTGCTGTTACGTATTATACTTATTGGCTTAGt                                                                                                                       >  1:488767/3‑141 (MQ=255)

                            tataATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTa                                                                                                                      >  2:697637/1‑141 (MQ=255)

                            tataATAATAGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                     <  1:538867/141‑2 (MQ=255)

                             ataataataGAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTACTCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                                    <  2:576585/141‑3 (MQ=255)

                                ataGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACg                                                                                                                  <  1:247138/141‑1 (MQ=255)

                                  aatagATTCGCCTATCTCTCAGGCG‑TCAA‑‑TGGACTAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                               <  2:770001/139‑1 (MQ=255)

                                  aatagAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                               <  2:749139/139‑1 (MQ=255)

                                   atagAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                                                                              >  2:565761/2‑140 (MQ=255)

                                    tagAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaaa                                                                                                             >  1:31180/1‑139 (MQ=255)

                                    tagAGTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                                                                             >  2:758991/1‑139 (MQ=255)

                                     agagTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑a                                                                                                           <  2:478470/141‑1 (MQ=255)

                                     agagTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑a                                                                                                           <  2:668581/141‑1 (MQ=255)

                                     agagTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                            <  1:244552/141‑4 (MQ=255)

                                     agagTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                            <  1:634147/141‑4 (MQ=255)

                                      gagTCGCCTATCCCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑Ac                                                                                                          <  2:41598/141‑1 (MQ=255)

                                       agTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                          >  2:72282/1‑136 (MQ=255)

                                       agTCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                          >  2:480994/1‑136 (MQ=255)

                                        gTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACaa                                                                                                        <  2:32723/141‑1 (MQ=255)

                                        gTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACaa                                                                                                        <  1:72687/141‑1 (MQ=255)

                                        gTCGCCGATCTCTCAGGCG‑TCAATT‑‑GAAGAAGAG‑AGGAGGTGCGT‑T‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                         >  2:196975/1‑135 (MQ=255)

                                         tCGCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                                                                                        >  2:113229/1‑134 (MQ=255)

                                           gCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                      >  1:287339/1‑132 (MQ=255)

                                           gCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                      <  2:154858/141‑10 (MQ=255)

                                           gCCTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                      <  2:155011/141‑10 (MQ=255)

                                            ccTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CC‑TCAC‑G‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCACCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataataa                                                                                                     >  1:730734/1‑131 (MQ=255)

                                            ccTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                                     >  1:491941/1‑131 (MQ=255)

                                            ccTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                                     >  2:514869/1‑131 (MQ=255)

                                            ccTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                                     >  2:616714/1‑131 (MQ=255)

                                            ccTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑T‑TACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                                     >  1:457055/1‑131 (MQ=255)

                                             cTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                    >  2:85138/1‑130 (MQ=255)

                                             cTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                    >  1:744668/1‑130 (MQ=255)

                                             cTATCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                    >  1:320741/1‑130 (MQ=255)

                                              tATCTCTCAGGCG‑TCAATT‑‑GTCGAAGAG‑AGGAGGTG‑CA‑T‑TA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                   <  2:71437/141‑13 (MQ=255)

                                               aTCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                  >  1:325699/1‑128 (MQ=255)

                                               aTCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                  >  2:277348/1‑128 (MQ=255)

                                               aTCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CATT‑AC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataagtaata                                                                                                  >  1:171442/1‑128 (MQ=255)

                                               aTCTCTCAGGCG‑TCAATT‑‑GACGAAGAG‑AGGAGGTG‑CA‑T‑TA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                  >  2:666719/1‑128 (MQ=255)

                                                tctctcAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTTTTTTTTTTTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTg                                                                                              >  2:205758/1‑141 (MQ=255)

                                                    tcAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACa‑t                                                                                         >  1:608235/1‑141 (MQ=255)

                                                    tcAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGAG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTAGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑CATGGGTGGCa‑t                                                                                         >  2:256390/1‑141 (MQ=255)

                                                      aGGCG‑TCAATTTAGACGCAGAG‑AGGAGGGG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGGGGCa‑tat                                                                                       >  2:587540/1‑141 (MQ=255)

                                                              aaTTTAGACGCAGAG‑AGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑TCAA‑‑AATAGGTGAGA‑TAT‑AGCC‑GCa                                                                              >  2:341349/1‑141 (MQ=255)

                                                                    gACGCAGAG‑AGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATa                                                                        >  1:456067/1‑141 (MQ=255)

                                                                       gctgag‑cgGAGGTG‑TA‑T‑AA‑G‑G‑GGATGCTTATTTTCGTTCCCATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATaaaa                                                                     <  2:379884/138‑1 (MQ=255)

                                                                            g‑aggagGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATccc                                                                <  1:624848/141‑1 (MQ=255)

                                                                              aggaggTG‑TG‑G‑AA‑G‑G‑TGATGCTTAGTTTCGTTCACATCATAGCACCCGTCATCAGGGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑CCAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATcctc                                                               >  1:566864/1‑139 (MQ=255)

                                                                              aggaggTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAGTAAACATcccc                                                               >  1:588547/1‑141 (MQ=255)

                                                                                gaggTG‑TA‑T‑AA‑G‑G‑TGATGATAGGGTGCTGTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTc                                                             >  2:748512/1‑141 (MQ=255)

                                                                                    tg‑tA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTa                                                         <  2:246895/141‑1 (MQ=255)

                                                                                         tt‑aA‑G‑G‑TCGTGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCg                                                      <  2:539066/140‑1 (MQ=255)

                                                                                                      tGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTATTTCTAATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAgag                                             >  1:658523/1‑139 (MQ=255)

                                                                                                      tGCTTATATTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAG‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAGTCCCCTCACTACCGCAAATAGTg                                             >  1:487649/1‑141 (MQ=255)

                                                                                                         ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAACAGTGAgg                                          >  1:49664/1‑141 (MQ=255)

                                                                                                              ttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGAtt                                     <  1:276995/141‑1 (MQ=255)

                                                                                                               ccgctCACATCATCGCTCTTGTCATGAGTGGCTGTGCCATTGCGTTTTTTTCGTCTGTGCTAAGTAGACGCAAT‑‑ATAA‑‑AACAGGTGGCA‑TAT‑AGCC‑TCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTg                                    <  2:213888/137‑1 (MQ=255)

                                                                                                                cGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTgg                                   >  1:237177/1‑141 (MQ=255)

                                                                                                                cGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACCCAAA‑‑ACAA‑‑AATCGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAAGGGATTgg                                   >  1:626273/1‑141 (MQ=255)

                                                                                                                  ttCACATCATAGCACCAGTCATCAGTGGCGGTGCGATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGGGACA‑TAC‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGCTTGgtg                                 <  2:268565/141‑1 (MQ=255)

                                                                                                                      catcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGtata                             >  2:291179/1‑141 (MQ=255)

                                                                                                                           tAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGA‑AAAAATCCCCACACTACTGC‑AATAGTGTGGGGATTGGTGcatt                             >  2:677055/1‑140 (MQ=255)

                                                                                                                             gCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑GATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAAt                      >  1:423073/1‑141 (MQ=255)

                                                                                                                              cACCGGTCATCAGTGGCTGTGCCATGGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAACTTCCCCTCACTTCCGCCAATAGTGAGGGGATTGGTGTATCAGTACATa                     <  2:377276/141‑1 (MQ=255)

                                                                                                                                  aGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTCTAAGTAAATACTTa                 >  2:118702/1‑141 (MQ=255)

                                                                                                                                    tcatcaGTGGCTGTGCCATTGCGTTTTTTTCTTCTTGGCTAAGTAGACGCAAT‑‑GCCA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCATTACCGCAAATAGTGAGGGGATTGTTGTATAAGTAAATACTTAtt               <  2:423686/141‑1 (MQ=255)

                                                                                                                                    tcatcaGTGGCTGGGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAtt               <  2:386619/141‑1 (MQ=255)

                                                                                                                                     catcaGTGGCTGTGCCATTGCGTTTTTTTTTTTTTTTCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAttt              >  2:302004/1‑141 (MQ=255)

                                                                                                                                          gTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AAAAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAACTACTTATTTTCGtt         >  2:627244/1‑141 (MQ=255)

                                                                                                                                                 tgCCATTGCGTTTTTTTCTAATTTGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCT‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATGCGTAAATCCTTATTTTGGTTGTCTTaa  >  1:409666/1‑141 (MQ=38)

GTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCG‑TCAATTTAGACGCAGAG‑AGGAGGTG‑TA‑T‑AA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAT‑‑ACAA‑‑AATAGGTGACA‑TAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAA  >  CP000730/1929777‑1930054

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